Georesistivity observational experiment based on encoded source
-
摘要: 地电阻率的高精度观测是实现地震预测预报的前提之一.面对日益严重的电磁干扰，提出了基于编码源循环互相关辨识技术的地电阻率观测方法，其实质是将待测地质体视为待辨识系统，利用编码源信号激励供电电极A和B产生的电流信号作为系统输入信号，测量电极M与N之间的电压信号作为系统响应输出，将输入和输出信号严格同步采样为时间序列，分别与参考信号进行循环互相关运算并转换至频率域计算获得待探测地电阻率谱(幅度和相位).由于系统环境的干扰和随机噪声与编码源信号不相关，通过循环互相关运算可以达到抑制环境随机噪声和干扰的目的.这种地电阻率观测体系在环境干扰较大的甘肃省兰州观象台和陇南汉王地震台站利用现有的观测场地和线路进行了观测试验，测量结果显示，数据的一致性好、均方差小，说明该方法在强干扰环境下具有较好的抗电磁干扰能力，观测频带较以往直流电法测量有较大的扩展.该方法为现有地面地电阻率台站持续发展提供了技术保障，可为地震预报与科学研究提供高质量的地电阻率观测数据.Abstract: The accurate observation of earth resistivity is an effective method for earthquake prediction. In the face of increasingly severe electromagnetic interference, a method of ground resistivity observation based on coded source circular cross-correlation identification is proposed. The measured geoelectrical structure can be treated as the system to be identified, using the coded current signal excited by power supply electrodes A and B as the input signal of the system, and the voltage signal between the electrodes M and N is the system output response, and the input and output signals are synchronously sampled into the time series, then a reference signal is generated according to sampled coded current time series. Next, the circular cross-correlation operation between input signal and the reference signal, and between output signal and the reference signal were performed separately, and the cross-correlation time series were converted into the frequency domain by fast Fourier transform and the apparent resistivity spectrum (amplitude and phase) of the earth to be detected can be deduced. A short experiment using this method was carried out at two stations of Lanzhou observatory and Hanwang station, the results showed that this method has a strong ability to restrain noise, and can achieve high frequency resolution and measurement accuracy. The multi-frequency apparent resistivity spectrum provides the conditions for the removal of electromagnetic coupling. The feasible technical scheme is put forward in this project, want to develop instrument, to optimize the main experimental parameters, and get a higher precision resistivity monitoring and the phase monitoring ability. The result can provide the data support for the study of the characteristics of time, location and intensity of seismic anomaly information and its relationship with strong earthquakes, and it also has practical significance and application value in the field of earthquake prevention and disaster reduction.
-
引言
地电阻率观测是开展地震预测预报的有效方法(赵玉林, 钱复业, 1978, 1979；钱复业，赵玉林，1980; 钱家栋，1993；钱家栋等，2013). 1966年以来，我国的地震地电阻率观测经历了从地表测量到井下观测的发展历程，在观测技术、理论方法、数据应用等方面均取得了较大进展(赵和云，1994; 赵家骝等，1996；汤吉等, 1998；Ma, 2002；Lu et al, 2004, 2016；黄清华，2005；Huang, 2008；张学民等，2009；杜学彬, 2010；高曙德等, 2010, 2013, 2016；王兰炜等, 2011；Zhao et al，2011；赵国泽等，2012；杜学彬等, 2015).目前, 地震台站受到日益严重的电磁干扰，对地电阻率观测造成了巨大的影响，干扰引起的视电阻率变化淹没了震前孕育的地电阻率变化(高曙德，2016).地震台站亟需引入新的观测技术和方法来提高地电阻率观测精度.
相关辨识法是一种可以有效去除随机噪声和干扰的系统辨识方法. Duncan等(1980)最早尝试了基于伪随机二进制序列(pseudo random binary sequence，简写为PRBS，简称m序列)为信号源的宽频带地电阻率测量方法，证明采用编码源进行大深度地电阻率观测的技术路线是可行的. Ziolkowski等(2007)基于PRBS编码源技术开发出了多道瞬变电磁法，对大埋深高阻薄层有较高的探测能力.国内近几年也掀起对基于PRBS编码源电磁勘探法理论和应用的研究热情. 赵璧如等(2006)介绍了基于CDMA(code division multiple access)技术的地电阻率高精度测量方法；汤井田和罗维斌(2008)以及罗维斌等(2012, 2016)提出了基于逆重复M序列伪随机频率域电磁勘探方法原理，并讨论了对各种干扰的压制方法以及提高测量精度的措施.逆重复M序列由2个PRBS序列(m序列)隔位取反产生，周期是m序列的2倍，自相关峰正反向交替，无直流分量，所以逆重复M序号是系统辨识一种比较理想的信号源.基于PRBS伪随机电法及观测系统研究也有新的进展(淳少恒等，2014；罗先中等，2014；罗延钟等, 2015, 2016)；齐彦福等(2015)利用方波响应移位叠加和电流导数与阶跃响应褶积两种方法实现了理论m序列和实际发射波形的全时正演模拟，武欣等(2015)研究了m序列伪随机编码源电磁响应的精细辨识问题；王若等(2016)利用伪随机编码源激发下的时域电磁信号合成与理论值对比，可用于检验后续数据处理方法的可靠性和有效性；张世中等(2001)采用有限长的伪随机码来控制供电周期(变周期供电测量法)以减少测量的随机误差, 并将其应用于ZD8B地电仪提高观测系统的测量精度；张宇等(2014, 2016)利用相关检测技术研制了交流地电阻率观测仪器，并在江宁台进行观测试验，其结果显示抗干扰能力得到了改善，但频段相对较窄.
近年来编码源电阻率测量方法在地学勘探领域得到了快速发展且应用广泛，主要采用m序列为激发信号源.该方法在地震系统目前尚处于试验阶段, 主要原因是没有配套的观测仪器，其次是现阶段地震系统原有的观测线路和观测场地无法改造，该方法能否观测到可靠连续的数据，还需要作进一步试验研究.本文基于2n系列伪随机电磁法理论(何继善，2010)，采用循环互相关法获取系统频率响应的方法(相关辨识方法是用伪随机序列(m序列或逆重复M序列)作为激励信号输入系统，由观测电场响应和输入电流信号辨识得到的地电系统频率响应，反映待探测大地系统的阻抗信息)，所获取的阻抗谱可直接用于反演解释，也可反算出系统的传递函数，从而获得地电系统的阶跃响应以及冲激响应，可得到复电阻率、阻抗及Cole-Cole模型参数(罗维斌, 李梅，2015).
本文基于相关辨识方法开展仪器的研制和实用化创新，在环境干扰较大的甘肃省兰州观象台和陇南汉王地震台利用现有的观测场地和线路进行试验观测，获取相关数据，以检验该方法在强干扰环境下的抗干扰能力，以期拓宽直流电法测量范围，提高地电阻率的观测精度.
1. 编码源理论和循环互相关法
1.1 基本原理
采用PRBS或逆重复M序列作为激励信号波形，利用互相关法辨识未知系统的频率响应，系统输入、输出的互相关为(刘立超等，2014；罗维斌等，2016)：
(1) 式中，hs(τ)为观测系统的冲激响应，he(τ)为待辨识未知系统的冲激响应，Ru(τ)为输入信号的自相关，Ryu(τ)为系统输入与输出信号的互相关函数.在复数域，Hs(s)和He(s)分别为观测系统自身hs(τ)和待辨识系统he(τ)的拉普拉斯变换, 式(1) 在频率域可变换为
(2) 式中，Wyu(jω)为系统输入与输出信号的互相关函数Ryu(τ)的拉普拉斯变换，Wu(jω)为系统输入信号自相关函数Ru(τ)的拉普拉斯变换，He(jω)为系统的频率响应.
基于循环互辨识法的工作流程如图 1所示.在计算数据中引入了一个参考信号S(按照输入激励信号的时间序列周期和采样率触发软件程序生成的无干扰信号)，与系统输入信号U(t)和系统输出信号yout(t)作循环互相关，由循环互相关时间序列的快速傅里叶变换得到互功率谱，由式(3) 和(4) 恢复出系统的频率响应H(ω)和相位φ(ω)，这样处理的优点是即使输入和输出两端都存在工频噪声, 也能得到有效抑制.
(3) 
(4) ![]()
图 1 基于逆重复M序列系统辨识技术流程图Figure 1. Workflow of system identification based on inverse repeated M-sequence式中，H(ω)为频率响应，PCys为参考信号S与输出信号yout的互功率谱，PCus为参考信号S与输入信号U(t)的互功率谱，φ(ω)为相位，mod为模运算符.
对于地电阻率谱而言，当观测值是接收电极间的电压值UMN时间序列时(MN为接收电场的测量间距)，以上计算过程可改写为
(5) 式中，PCUMNs(ω)为参考信号S与接收电压信号UMN的互功率谱，PCIABs(ω)为参考信号S与供电电流信号IAB(t)的互功率谱，
(6) 式中K为观测装置系数
(7) 式中，A和B为供电电极，M和N为测量电极.由于采用互相关技术，本测量方法对观测场地干扰有较强的压制能力.
1.2 循环互相关法
对于两个时间序列的互相关计算，采用非循环互相关法计算时，在序列移位过程中需要在序列头部和尾部补零来实现两个序列长度相等，计算结果会在相关脉冲峰两侧产生误差干扰；循环互相关法计算时，将等长的两个时间序列分别首尾环状相连，这样在计算中序列循环移位，无需补零(Roinila et al, 2009)，且互相关峰值相等，这个特性在多周期信号功率谱叠加中优势显著.
2. 试验结果与分析
2.1 观测方案
利用甘肃省地震局地电阻率台站的观测装置开展试验, 仪器选用由罗维斌等(2016)试制的测量样机，数据采样率为1000 Hz，对发送电流和接收电压多周期同步采集为时间序列，研究同步采样和非同步采样对于辨识结果的影响，通过计算观测结果来调整和优化仪器参数，从而确定最佳的观测参数.
2.2 兰州观象台试验
2016年7月30日选用兰州观象台进行试验.该台站位于黄河Ⅲ级阶地(次生黄土)，第四纪覆盖层厚10—30 m，岩性为亚砂土、粉砂、沙砾石，下伏砂岩.地下水位埋深5—30 m，含水层厚度大于5 m，含水层岩性为粉砂和沙砾石; 台站的电测曲线为Q型(李乐进，1985). 1970年5月开始观测, 布设有3道测线，其方位分别为N5°W, N40°E，N50°W，近年受站外蔬菜大棚建设和工业及城镇化建设的干扰影响增大, 测量数据波动较大(图 2)，其中2015年初至2016年初，N5°W连续10个月出现脉冲状干扰，最大干扰值可达150 Ω·m，测值曲线无法清楚地看到工程焊接等干扰，所以在绘图时对图 2作了限幅来显示工程建设和焊接的干扰. 2016年初N40°E测道供电极B附近电焊作业，资料出现振荡变化；2016年9月开始建设和工程焊接使N5°W测值振荡变化持续了8个月，给资料分析带来了不确定的因素, 资料的可用性大大降低，该台目前已规划选址搬迁.本次试验利用N40°E现有观测线路(AB=900 m，MN=300 m), 用8号镀锌铁丝悬挂电缆线作为测量及供电导线，电极为铅板，埋深2.5 m.
![]()
图 2 2015年1月—2017年6月兰州观象台直流电法观测的地电阻率曲线Figure 2. The georesistivity curves recorded at the Lanzhou observatory by direct current method from January 2015 to June 20172.2.1 输入信号和输出电压
为了验证信号序列长度压制干扰的能力，设计编码源信号产生参数：寄存器位数m=7，序列周期254位；晶振频率5 MHz，分频后获得5.24 Hz时钟频率激励移位寄存器，产生的逆重复M序列周期为63.5 s.直流电压110 V，供电电流1.5 A. 图 3为试验振荡器产生的逆M序列信号波形、功率谱和循环自相关时间序列.
![]()
图 3 逆M序列信号波形、功率谱和循环自相关时间序列Figure 3. Signal waveform, power spectrum and circular autocorrelation of inverse repeated M sequence由于本项工作是试验性的，在台站未敷设新的观测线路，只有借用工作的线路在正常观测的间隙进行试验，且由于试验样机的参数是通过模拟得到，需在实际观测中经不断地调整和改正来逐步完善，所以观测进行了一天. 2016年7月30日9时开始试验，由于整点台站的其它仪器工作，所以从每个整点后30分钟开始工作，记录时长为20分钟，试验工作至17时50分结束；用1000 Hz采样频率同时记录激发电流信号和接收电压信号.从记录的电压曲线(图 4)看出，测量信号峰峰电压可达3500 mV(各种干扰及开关瞬时电磁感应、介质响应等综合值)，观测的编码信号地电介质响应电压为500 mV (图 4中黑点)，这两个信号的强度相差7倍，所以利用传统的地电观测方法，当地下介质的电性发生微弱变化时，基本上被噪声信号所淹没，想检测出目标层电阻率的真实变化则比较困难，而利用编码源的信号互相关理论，可以有效压制不相关噪声，分辨出有用信息，分辨信号能力较直流电法有一定的改善.
![]()
图 4 兰州观象台记录的测量电压(空心圆)和加载的编码源信号的电压(黑点)Figure 4. The value of ambient noise (open circles) and the encoding pseudorandom signals (dark dots) recorded at Lanzhou observatory2.2.2 地电阻率计算
对图 4的电压信号按图 1的流程使用上述数据处理方法，产生参考信号，然后通过循环互相关法计算输入、输出信号与参考信号的互相关函数，再转换至频率域获得阻抗频率响应，乘以装置系数K即可得到不同频率的地电阻率，结果如图 5所示.可见：在不同时段、同一地点测量的地电阻率变化形态相同；以频率1 Hz为分界线，频率小于1 Hz的地电阻率变化相对稳定，地电阻率在6—10 Ω·m范围内变化，且随着频率的减小，地电阻率也在变小，反映了地下介质的电性结构在深部是低阻层，这与该台站电性分层一致(李乐进，1985)：表层为第四纪覆盖层(亚砂土、粉砂、沙砾石)，随着深度增加下伏为砂岩. 20世纪70年代该地区无电磁干扰测量电测曲线为Q型(钱复业等，1987)，根据该台站电测深曲线反演获得上部3层地电断面的参数：h1=3.5 m, ρ1= 80 Ω·m；h2= 15 m, ρ2=30 Ω·m；h3=300 m，ρ3=9 Ω·m.随频率增大，基于1 Hz以上频率点计算的地电阻率数值变大，分析认为有两个原因：其一，随着频率的增大，反映表层的信息更加丰富；其二，供电线路与接收线路在一起平行布线产生了电磁耦合. 图 6给出了2016年7月30日兰州观象台在不同频率点(0.05—1 Hz)对应时段记录的地电阻率曲线，可以看出低频点(0.05—0.2 Hz)的地电阻率与直流观测数值基本接近(图 2)，且受电磁耦合影响小，反映出地电阻率随时间的变化，且随着频率的变化，电阻率的变化层次也十分清晰，勾勒出不同频率电阻率的分布和变化，为研究目标层的电性精细变化提供了基础.
![]()
图 5 兰州观象台在不同时段记录0.01—25 Hz频率点的地电阻率曲线Figure 5. Georesistivity spectrum during 0.01-25 Hz at different time observed at Lanzhou observatory![]()
图 6 2016年7月30日兰州观象台在不同频率点(0.01—1.03 Hz)对应时段记录的地电阻率曲线Figure 6. The georesistivity histogram at different time moment recorded at Lanzhou obser-vatory on July 30, 2016 for different frequencies (0.01-1.03 Hz)2.3 陇南汉王台试验
陇南汉王地电台地处南北地震带中段，武都山字型弧状顶端转折部位.测区地形北高南低, 西高东低，北靠山、南临白龙江，覆盖厚度由东至西、由南向北从约22 m增至100 m以上，表层为碎砾石沙土，下伏中上志留纪碳硅质板岩、含黄铁矿片岩和磁黄铁矿砂岩互层，测区位于磁异常区附近，电测深曲线为Q型①. 1974年建台布设3个地电阻率测道: N85°E道和N54°W道(AB=1 km，MN=0.333 km), N73°W道(AB=0.924 km，MN=0.308 km).近年来该地电测区由于受到铁路和高速公路、电站及工程建设的影响，观测资料质量大大降低. 图 7为2013—2014年N85°E测道的地电阻整点值曲线. 2013年7月开始由于电站的工作附加高频干扰，8月底A供电极附近填土建设电阻率曲线上升一个台阶；由于兰渝铁路汉王段通过汉王地电台的测区，其中铁路路基横跨N85°E道和N73°W道，共用供电电极(B1和B3)需要搬迁，故在2013年将该电极移至铁路的同侧路基下方；2014年2月初开始，铁路建设开始筑土碾压路基，在地面加高了3—4 m，在此过程中电阻率出现了持续振荡变化；由于测区环境严重干扰，2015年地面观测被迫停测.本次试验就选用2016年10月3—4日N85°E测道旧线路在原场地进行，以检验仪器的抗干扰性和数据观测质量.
① 国家地震局兰州地震研究所台群处编. 1983.甘肃省地震台志(中).
![]()
图 7 2013—2014年陇南汉王台N85°E测道地电阻率观测曲线Figure 7. The georesistivity curve recorded by the channel of N85°E at Hanwang observatory in 2013-20142.3.1 输入信号和输出电压
为了更利于压制干扰，通过仿真优选，使产生的信号序列周期是工频的整数倍，来设置编码源信号产生参数：寄存器位数m=6，序列周期126位；晶振频率6.3 MHz，分频后获得15 Hz时钟频率激励移位寄存器，产生的逆重复M序列周期为16.8 s.直流电压110 V，供电电流1.8 A.
2016年10月3日9时开始试验，从每个整点后20分开始工作，记录时间30分钟，工作至晚上8时结束；使用2个采集器按1000 Hz采样频率记录输入电流和输出电压信号.
2.3.2 地电阻率计算
地电阻率计算按照式(1)—(6) 计算.从2016年10月3日和4日两天的地电阻率观测(图 8)可看出，除10月3日10时的地电阻率测值变化较大外，其它时段测值比较稳定，10月3日0.06 Hz地电阻率在21.3 Ω·m上下波动，较平稳，10月4日0.06 Hz地电阻率在21.35 Ω·m上下波动，更为平稳.其它频率显示出相似的特征，高频点由于电磁耦合阻值升高，1.07 Hz的地电阻率在24.68 Ω·m上下波动.对10月3—4日0.06—1.07 Hz频率范围内电阻率变化的分析结果(表 1)表明，10月3日地电阻率的均值、变化幅度、变化率(日变化幅值与日观测均值的比值乘以100%)均大于4日，其原因在于3日数据采集没有同步，4日同步采样后变幅明显减小，各频率点的日变化率均小于1%.随着频率减小，地电阻率也减小，反映了陇南汉王地电测区地层的电性变化，表层为碎砾石沙土覆盖层，下伏含黄铁矿片岩和磁黄铁矿砂岩互层碳硅质板岩低阻层，与20世纪70年代测区无电测干扰测量电测深曲线Q型变化形态相吻合.
![]()
图 8 2016年10月3日(a)和4日(b)陇南汉王台在不同时段记录不同频率点的地电阻率曲线Figure 8. The georesistivity curves at different frequency at Hanwang observatory on October 3 (a) and 4 (b), 2016表 1 2016年10月3日和4日陇南汉王台测试地电阻率变化幅值及变化率Table 1. The amplitude and variation of georesistivity ratio on October 3 and 4, 2016 at Hanwang observatory, Longnan city频率/Hz 2016年10月3日 2016年10月4日 ρs/(Ω·m) Δρs/(Ω·m) Δρs /ρs ρs/(Ω·m) Δρs/(Ω·m) Δρs /ρs 0.06 21.43 0.65 3.03 21.38 0.15 0.70 0.19 21.50 0.25 1.16 21.47 0.21 0.98 0.31 21.69 0.42 1.92 21.69 0.09 0.41 0.44 21.98 0.36 1.63 22.02 0.17 0.77 0.57 22.40 0.55 2.45 22.42 0.15 0.67 0.69 22.79 1.09 4.76 22.91 0.12 0.52 0.82 23.32 1.26 5.39 23.44 0.19 0.81 0.94 23.86 1.56 6.54 24.06 0.21 0.87 1.07 24.44 2.13 8.72 24.72 0.18 0.73 3. 讨论与结论
甘肃省兰州观象台和陇南汉王地震台站的试验验证了本文所述的编码源电阻率测量原理是可行的，利用编码源地电仪在原有线路和场地可观测到高质量的地电阻率数据.尽管试验时间短，不能刻画出地电阻率的长期变化，但是通过试验得到的数据说明因干扰停测的汉王地电观测台，用编码源地电仪可以在原测区、原线路下工作，这说明改进的观测仪器方案是能够推广的.试验工作取得的成果如下：
1) 将相关辨识技术应用于地电阻率观测，可极大地压制观测场地的电磁干扰，提高地电阻率观测精度；
2) 由于采用编码源宽频信号激发，可得到一系列频率的地电阻率谱，利用这些地电阻率谱以期获得其它反映地震孕育电性变化的指标参数；
3) 通过两个地电台站试验，形成了地电阻率高精度观测方案，由辨识结果实时处理得到地电阻率随时间的变化形态，且随着频率的变化勾勒出不同频率电阻率的分布，对微小的变化具有较强的分辨力, 为研究目标层的电性变化提供了基础；
4) 获得了一些优化参数，运用这些参数改进编码源仪器可观测到高质量数据.两次试验结果表明，编码源产生参数如序列周期长度、叠加周期数、采样率等对观测结果均有明显影响.若序列周期是工频信号周期的偶数倍，将有助于提升抗干扰能力，故选序列长度以6阶寄存器产生为宜，时钟频率使产生的序列周期须在9—36 s范围内选取.
后续工作将通过仪器研发和观测试验完善观测系统.由于待监测地电系统非常复杂，地质情况也复杂，观测装置及布极的深度位置和极距大小均会影响到电流场在地电系统中的扩散路径，因此利用地震台站的现有观测装置，进一步研究不同深度电极装置对辨识结果的影响，对于地电阻率高精度观测有现实意义.
进一步研究去电磁耦合问题.电磁耦合是地电观测普遍存在的问题，电磁耦合去耦好坏会对地电阻率的观测精度产生直接的影响，利用编码源宽域广谱的特点，一次供电可获得一系列频率的地电阻率谱，为电磁去耦提供了条件(可用多频去耦，该方法后续撰文介绍)；相位信息也能较灵敏地反映地电阻率变化，所以本次试验对相位的观测进行了尝试，需要收发装置的采样频率严格同步，才能满足观测要求.目前通过方法技术创新研发的编码源地电阻率观测仪器已在装配，所形成的技术方法可以为利用现有观测装置进行地震台站地电阻率观测提供技术支持和实践基础.
本研究在实施过程中得到了中国地震局地质研究所赵国泽研究员、预测研究所钱家栋研究员、北京大学黄清华教授、长江大学胡文宝教授的关注, 得到了甘肃省地震局兰州观象台和陇南地震台的大力协助, 两位匿名审稿专家对本文提出了建设性的修改意见和建议，作者在此一并表示感谢.
-
![]()
图 1 基于逆重复M序列系统辨识技术流程图
Figure 1. Workflow of system identification based on inverse repeated M-sequence
![]()
图 2 2015年1月—2017年6月兰州观象台直流电法观测的地电阻率曲线
Figure 2. The georesistivity curves recorded at the Lanzhou observatory by direct current method from January 2015 to June 2017
![]()
图 3 逆M序列信号波形、功率谱和循环自相关时间序列
Figure 3. Signal waveform, power spectrum and circular autocorrelation of inverse repeated M sequence
![]()
图 4 兰州观象台记录的测量电压(空心圆)和加载的编码源信号的电压(黑点)
Figure 4. The value of ambient noise (open circles) and the encoding pseudorandom signals (dark dots) recorded at Lanzhou observatory
![]()
图 5 兰州观象台在不同时段记录0.01—25 Hz频率点的地电阻率曲线
Figure 5. Georesistivity spectrum during 0.01-25 Hz at different time observed at Lanzhou observatory
![]()
图 6 2016年7月30日兰州观象台在不同频率点(0.01—1.03 Hz)对应时段记录的地电阻率曲线
Figure 6. The georesistivity histogram at different time moment recorded at Lanzhou obser-vatory on July 30, 2016 for different frequencies (0.01-1.03 Hz)
![]()
图 7 2013—2014年陇南汉王台N85°E测道地电阻率观测曲线
Figure 7. The georesistivity curve recorded by the channel of N85°E at Hanwang observatory in 2013-2014
![]()
图 8 2016年10月3日(a)和4日(b)陇南汉王台在不同时段记录不同频率点的地电阻率曲线
Figure 8. The georesistivity curves at different frequency at Hanwang observatory on October 3 (a) and 4 (b), 2016
表 1 2016年10月3日和4日陇南汉王台测试地电阻率变化幅值及变化率
Table 1 The amplitude and variation of georesistivity ratio on October 3 and 4, 2016 at Hanwang observatory, Longnan city
频率/Hz 2016年10月3日 2016年10月4日 ρs/(Ω·m) Δρs/(Ω·m) Δρs /ρs ρs/(Ω·m) Δρs/(Ω·m) Δρs /ρs 0.06 21.43 0.65 3.03 21.38 0.15 0.70 0.19 21.50 0.25 1.16 21.47 0.21 0.98 0.31 21.69 0.42 1.92 21.69 0.09 0.41 0.44 21.98 0.36 1.63 22.02 0.17 0.77 0.57 22.40 0.55 2.45 22.42 0.15 0.67 0.69 22.79 1.09 4.76 22.91 0.12 0.52 0.82 23.32 1.26 5.39 23.44 0.19 0.81 0.94 23.86 1.56 6.54 24.06 0.21 0.87 1.07 24.44 2.13 8.72 24.72 0.18 0.73 
下载: 导出CSV
-
淳少恒, 陈儒军, 耿明会. 2014.伪随机m序列及其在电法勘探中的应用进展[J].地球物理学进展, 29(1): 439-446. doi: 10.6038/pg20140164. Chun S H, Chen R J, Geng M H. 2014. Review of the pseudo-random m sequence and its application in electrical prospecting of exploration geophysics[J]. Progress in Geophysics, 29(1): 439-446. doi: 10.6038/pg20140164 (in Chinese).
杜学彬. 2010.在地震预报中的两类视电阻率变化[J].中国科学:地球科学, 40(10): 1321-1330. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201010004.htm Du X B. 2011. Two types of changes in apparent resistivity in earthquake prediction[J]. Science China Earth Sciences, 54(1): 145-156. doi: 10.1007/s11430-010-4031-y
杜学彬, 刘君, 崔腾发, 范莹莹, 安张辉, 闫睿, 王丽. 2015.两次近距离大震前成都台视电阻率重现性、相似性和各向异性变化[J].地球物理学报, 58(2): 576-588. doi:10.6038/cjg2015022. Du X B, Liu J, Cui T F, Fan Y Y, An Z H, Yan R, Wang L. 2015. Repeatability, similarity and anisotropy changes in apparent resistivity recorded by station Chengdu at near distances before two great earthquakes[J]. Chinese Journal of Geophysics, 58(2): 576-588. doi: 10.6038/cjg20150220 (in Chinese).
高曙德, 汤吉, 杜学彬, 刘小凤, 苏永刚, 陈彦平, 狄国荣, 梅东林, 詹艳, 王立凤. 2010.汶川8.0级地震前后电磁场的变化特征[J].地球物理学报, 53(3): 512-525. doi: 10.3969/j.issn.0001-5733.2010.03.005. Gao S D, Tang J, Du X B, Liu X F, Su Y G, Chen Y P, Di G R, Mei D L, Zhan Y, Wang L F. 2010. The change characteristics of electromagnetic field before to after Wenchuan MS8.0 earthquake[J]. Chinese Journal of Geophysics, 53(3): 512-525. doi: 10.3969/j.issn.0001-5733.2010.03.005 (in Chinese).
高曙德, 汤吉, 孙维怀. 2013.盈江5.8级和缅甸7.2级地震前电磁异常[J].地球物理学报, 56(5): 1538-1548. doi:10.6038/cjg20130512. Gao S D, Tang J, Sun W H. 2013. Electromagnetic anomaly before the Yingjiang MS5.8 and Myanmar MS7.2 earthquakes[J]. Chinese Journal of Geophysics, 56(5): 1538-1548. doi:10.6038/cjg20130512 (in Chinese).
高曙德. 2016.深井地电观测技术在地震监测中的应用探讨[J].地球物理学进展, 31(5): 2078-2088. doi: 10.6038/pg20160526. Gao S D. 2016. Discussion on the deep well geoelectric observation technique applied in earthquake monitoring[J]. Progress in Geophysics, 31(5): 2078-2088. doi: 10.6038/pg2016052 (in Chinese).
高曙德, 赵国泽, 汤吉, 苏永刚, 詹艳, 王立凤. 2016. 2008年中国大陆6.0级以上地震前后电磁脉冲异常现象[J].地震地质, 38(4): 987-1004. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ201604015.htm Gao S D, Zhao G Z, Tang J, Su Y G, Zhan Y, Wang L F. 2016. Electromagnetic anomalous pulses before and after earthquakes (M > 6.0) in China continent in 2008[J]. Seismology and Geology, 38(4): 987-1004 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ201604015.htm
何继善. 2010.广域电磁法和伪随机信号电法[M].北京:高等教育出版社: 1-214. He J S. 2010. Wide Field Electromagnetic Sounding Methods and Pseudo-Random Signal Coding Electrical Method[M]. Beijing: Higher Education Press: 1-214 (in Chinese).
黄清华. 2005.地震电磁信号传播的控制模拟实验[J].科学通报, 50(16): 1774-1778. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB200516016.htm Huang Q H. 2005. Controlled analogue experiments on propagation of seismic electromagnetic signals[J]. Chinese Science Bulletin, 50(17): 1957-1961. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB200516016.htm
李乐进. 1985.甘肃省地电台址条件研究[J].西北地震学报, 7(2): 64-76. http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ198502008.htm Li L J. 1985. Study on the conditions of the station location of telluric current in Gansu Province[J]. Northwestern Seismological Journal, 7(2): 64-76 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ198502008.htm
刘立超, 丁凯来, 林君, 周逢道, 刘长胜. 2014.基于伪随机系统辨识的电磁法仪器标定[J].仪器仪表学报, 35(8): 1721-1728. http://www.cnki.com.cn/Article/CJFDTOTAL-YQXB201408006.htm Liu L C, Ding K L, Lin J, Zhou F D, Liu C S. 2014. Instrument calibration with electromagnetic method based on pseudorandom system identification[J]. Chinese Journal of Scientific Instrument, 35(8): 1721-1728 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-YQXB201408006.htm
罗维斌, 李庆春, 汤井田. 2012.编码电磁测深[J].地球物理学报, 55(1): 341-349. doi: 10.6038/j.issn.0001-5733.2012.01.035. Luo W B, Li Q C, Tang J T. 2012. Coded source electromagnetic sounding method[J]. Chinese Journal of Geophysics, 55(1): 341-349. doi: 10.6038/j.issn.0001-5733.2012.01.035 (in Chinese).
罗维斌, 李梅. 2015. 一种获取待辨识系统频率响应函数的方法及装置: 中国, CN104794313A[P]. [2015-07-22]. http://www.pss-system.gov.cn/sipopublicsearch/patentsearch/portal2HomeSearch-portalSearch.shtml. Luo W B, Li M. 2015. Method and Device for Obtaining Frequency Response Function of System to be Identified: China, CN104794313A[P]. [2015-07-22]. http://www.pss-system.gov.cn/sipopublicsearch/patentsearch/portal2HomeSearch-portalSearch.shtml (in Chinese).
罗维斌, 李梅, 韩海涛, 高曙德. 2016.一种物探仪器系统特性高精度标定方法[J].物探与化探, 40(4): 816-819. http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH201604030.htm Luo W B, Li M, Han H T, Gao S D. 2016. A high accuracy calibration method for geophysical prospecting instrument system characteristics[J]. Geophysical and Geochemical Exploration, 40(4): 816-819 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH201604030.htm
罗先中, 李达为, 彭芳苹, 罗林涛, 赵金水. 2014.抗干扰编码电法仪的实现及应用[J].地球物理学进展, 29(2): 944-951. doi: 10.6038/pg20140263. Luo X Z, Li D W, Peng F P, Luo L T, Zhao J S. 2014. Implementation and applications of an coded electrical instrument with anti-interference ability[J]. Progress in Geophysics, 29(2): 944-951. doi: 10.6038/pg20140263 (in Chinese).
罗延钟, 陆占国, 孙国良, 黄伟才, 高隆钦. 2015.新一代电法勘查仪器:伪随机信号电法仪[J].地球物理学进展, 30(1): 411-415, doi: 10.6038/pg20150159. Luo Y Z, Lu Z G, Sun G L, Huang W C, Gao L Q. 2015. New generation of instruments for electrical and electromagnetic prospecting: Instruments using pseudo random signal[J]. Progress in Geophysics, 30(1): 411-415. doi: 10.6038/pg20150159 (in Chinese).
罗延钟, 陆占国, 孙国良, 黄伟才, 高隆钦. 2016.伪随机信号电法仪器的"第二方案"算法[J].地球物理学进展, 31(6): 2604-2608. doi: 10.6038/pg20160634. Luo Y Z, Lu Z G, Sun G L, Huang W C, Gao L Q. 2016. Algorithms of second solution for instruments of electrical and electromagnetic prospecting using pseudo random signal[J]. Progress in Geophysics, 31(6): 2604-2608. doi: 10.6038/pg20160634 (in Chinese).
齐彦福, 殷长春, 王若, 蔡晶. 2015.多通道瞬变电磁m序列全时正演模拟与反演[J].地球物理学报, 58(7): 2566-2577. doi: 10.6038/cjg20150731. Qi Y F, Yin C C, Wang R, Cai J. 2015. Multi-transient EM full-time forward modeling and inversion of m-sequences[J]. Chinese Journal of Geophysics, 58(7): 2566-2577. doi: 10.6038/cjg20150731 (in Chinese).
钱复业, 赵玉林. 1980.地震前地电阻率变化十例[J].地震学报, 2(2): 186-197. http://www.dzxb.org/Magazine/Show?id=28655 Qian F Y, Zhao Y L. 1980. Ten examples of changes in earth-resistivity prior to strong earthquakes[J]. Acta Seismo-logica Sinica, 2(2): 186-197 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28655
钱复业, 赵玉林, 许同春. 1987.地电阻率季节干扰变化分析[J].地震学报, 9(3): 289-302. http://www.dzxb.org/Magazine/Show?id=28209 Qian F Y, Zhao Y L, Xu T C. 1987. Analyses of the seasonally variation of disturbance in georesistivity[J]. Acta Seismologica Sinica, 9(3): 289-302 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28209
钱家栋. 1993.与大震孕育过程有关的地电阻率变化研究[J].中国地震, 9(4): 341-350. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD199304006.htm Qian J D. 1993. A study on the changes in geoelectrical resistivity associated with preparatory process of great earthquakes in China[J]. Earthquake Research in China, 9(4): 55-64. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD199304006.htm
钱家栋, 马钦忠, 李劭秾. 2013.汶川MS8.0地震前成都台NE测线地电阻率异常的进一步研究[J].地震学报, 35(1): 4-17. http://www.dzxb.org/Magazine/Show?id=28812 Qian J D, Ma Q Z, Li S N. 2013. Further study on the anomalies in apparent resistivity in the NE configuration at Chengdu station associated with Wenchuan MS8.0 earthquake[J]. Acta Seismologica Sinica, 35(1): 4-17 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28812
汤吉, 赵国泽, 王继军, 李文军, 詹艳. 1998.张北—尚义地震前后电阻率的变化及分析[J].地震地质, 20(2): 164-171. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ802.011.htm Tang J, Zhao G Z, Wang J J, Li W J, Zhan Y. 1998. Variation and analysis of resistivity before and after the Zhangbei-Shangyi earthquake[J]. Seismology and Geology, 20(2): 164-171 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ802.011.htm
汤井田, 罗维斌. 2008.基于相关辨识的逆重复m序列伪随机电磁法[J].地球物理学报, 51(4): 1226-1233. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200804036.htm Tang J T, Luo W B. 2008. Pseudo-random electromagnetic exploration based on invert-repeated m-sequence correlation identification[J]. Chinese Journal of Geophysics, 51(4): 1226-1233 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200804036.htm
王兰炜, 朱旭, 朱涛, 张世中, 刘大鹏, 胡哲, 张宇. 2011.地电阻率多极距观测系统及试验研究[J].地震, 31(1): 20-31. http://youxian.cnki.com.cn/yxdetail.aspx?filename=NMHB201706095&dbname=CJFDPREP Wang L W, Zhu X, Zhu T, Zhang S Z, Liu D P, Hu Z, Zhang Y. 2011. Multi-separation array geo-electrical resistivity observation system and its experimental observation[J]. Earthquake, 31(1): 20-31 (in Chinese). http://youxian.cnki.com.cn/yxdetail.aspx?filename=NMHB201706095&dbname=CJFDPREP
王若, 王妙月, 底青云, 薛国强, 殷长春, 雷达. 2016.伪随机编码源激发下的时域电磁信号合成[J].地球物理学报, 59(12): 4414-4423. doi: 10.6038/cjg20161203. Wang R, Wang M Y, Di Q Y, Xue G Q, Yin C C, Lei D. 2016. Signal synthesis of time-domain electromagnetic response excited by a pseudo random binary sequence source[J]. Chinese Journal of Geophysics, 59(12): 4414-4423. doi: 10.6038/cjg20161203 (in Chinese).
武欣, 薛国强, 底青云, 张一鸣, 方广有. 2015.伪随机编码源电磁响应的精细辨识[J].地球物理学报, 58(8): 2792-2802. doi: 10.6038/cjg20150814. Wu X, Xue G Q, Di Q Y, Zhang Y M, Fang G Y. 2015. Accurate identification for the electromagnetic impulse response of the earth with pseudo random coded waveforms[J]. Chinese Journal of Geophysics, 58(8): 2792-2802. doi: 10.6038/cjg20150814 (in Chinese).
张世中, 赵家骝, 王燕琼. 2001.一种消除地电阻率观测中随机误差的新方法[J].西北地震学报, 23(1): 88-91. http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ200101016.htm Zhang S Z, Zhao J L, Wang Y Q. 2001. A new method to remove stochastic error in measuring apparent resistivity[J]. Northwestern Seismological Journal, 23(1): 88-91 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ200101016.htm
张学民, 李美, 关华平. 2009.汶川8.0级地震前的地电阻率异常分析[J].地震, 29(1): 108-115. doi: 10.3969/j.issn.1000-3274.2009.01.014. Zhang X M, Li M, Guan H P. 2009. Anomaly analysis of earth resistivity observations before the Wenchuan earthquake[J]. Earthquake, 29(1): 108-115. doi: 10.3969/j.issn.1000-3274.2009.01.014 (in Chinese).
张宇, 王兰炜, 张兴国, 朱旭, 刘大鹏, 颜蕊. 2014.相关检测技术在低频交流地电阻率观测中的应用[J].地球物理学进展, 29(4): 1973-1979. doi: 10.6038/pg20140468. Zhang Y, Wang L W, Zhang X G, Zhu X, Liu D P, Yan R. 2014. Application of correlation detection technology in low-frequency AC geo-resistivity observation[J]. Progress in Geophysics, 29(4): 1973-1979. doi: 10.6038/pg20140468 (in Chinese).
张宇, 张兴国, 王兰炜, 马小溪, 赵庆福, 袁慎杰, 王子影. 2016.新型地电阻率交流观测系统研究及江宁台观测试验[J].地震学报, 38(5): 807-810. doi: 10.11939/jass.2016.05.015. Zhang Y, Zhang X G, Wang L W, Ma X X, Zhao Q F, Yuan S J, Wang Z Y. 2016. A new AC geo-electrical resistivity observation system and experimental observation in Jiangning seismic station[J]. Acta Seismologica Sinica, 38(5): 807-810. doi: 10.11939/jass.2016.05.015 (in Chinese).
赵璧如, 赵健, 张洪魁, 钱卫, 赵玉林, 钱复业. 2006. PS100型IP到端可控源高精度大地电测仪系统: CDMA技术首次在地电阻率测量中的应用[J].地球物理学进展, 21(2): 675-682. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ200602052.htm Zhao B R, Zhao J, Zhang H K, Qian W, Zhao Y L, Qian F Y. 2006. The PS100 high precision earth-electricity instrument system (IP to IP) with controllable source: Application of CDMA technology to the measurement of earth-resistivity for the first time[J]. Progress in Geophysics, 21(2): 675-682 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ200602052.htm
赵国泽, 王立凤, 詹艳, 汤吉, 肖骑彬, 陈小斌, 王继军, 蔡军涛, 汪晓, 杨静. 2012.地震预测人工源极低频电磁新技术(CSELF)和第一个观测台网[J].地震地质, 34(4): 576-585. doi: 10.3969/j.issn.0253-4967.2012.04.004. Zhao G Z, Wang L F, Zhan Y, Tang J, Xiao Q B, Chen X B, Wang J J, Cai J T, Wang X, Yang J. 2012. A new electromagnetic technique for earthquake monitoringCSELF and the first observational network[J]. Seismology and Geology, 34(4): 576-585. doi: 10.3969/j.issn.0253-4967.2012.04.004 (in Chinese).
赵和云. 1994.地电阻率趋势变化的形态特征与地震[J].地震学报, 16(3): 368-375. http://www.dzxb.org/Magazine/Show?id=28627 Zhao H Y. 1994. The pattern characteristics of the tendency variations of earth resistivity and its relation to earthquakes[J]. Acta Seismologica Sinica, 16(3): 368-375 (in Chinese). http://www.dzxb.org/Magazine/Show?id=28627
赵家骝, 李海亮, 李建, 王燕琼. 1996.用正向供电技术提高视电阻率观测精度[J].西北地震学报, 18(2): 31-36. http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ602.005.htm Zhao J L, Li H Y, Li J, Wang Y Q. 1996. Improving the observation accuracy of apparent resistivity by using up-down transmission[J]. Northwestern Seismological Journal, 18(2): 31-36 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ602.005.htm
赵玉林, 钱复业. 1978.唐山7.8级强震前震中周围形变电阻率的下降异常[J].地球物理学报, 21(3): 181-190. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX197803000.htm Zhao Y L, Qian F Y. 1978. Electrical resistivity anomaly observed in and around the epicentral area prior to the Tang-shan earthquake of 1976[J]. Acta Geophysica Sinica, 21(3): 181-190 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX197803000.htm
赵玉林, 钱复业. 1979.唐山地震临震前后形变电阻率的脉冲[J].地球物理学报, 22(1): 25-31. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX197901003.htm Zhao Y L, Qian F Y. 1979. Short duration resistivity variations before and after the M=7.8 Tangshan earthquake[J]. Acta Geophysica Sinica, 22(1): 25-31 (in Chinese). http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX197901003.htm
Duncan P M, Hwang A, Edwards R N, Bailey R C, Garland G D. 1980. The development and applications of a wide band electromagnetic sounding system using a pseudo-noise source[J]. Geophysics, 45(8): 1276-1296. doi: 10.1190/1.1441124
Huang Q H. 2008. Seismicity changes prior to the MS8.0 Wenchuan earthquake in Sichuan, China[J]. Geophys Res Lett, 35(23): L23308. doi: 10.1029/2008GL036270.
Lu J, Xue S Z, Qian F Y, Zhao Y L, Guan H P, Mao X Q, Ruan A G, Yu S R, Xiao W J. 2004. Unexpected changes in resistivity monitoring for earthquakes of the Longmen Shan in Sichuan, China, with a fixed Schlumberger sounding array[J]. Phys Earth Planet Inter, 145(14): 87-97. http://www.sciencedirect.com/science/article/pii/S0031920104001621
Lu J, Xie T, Li M, Wang Y L, Ren Y X, Gao S D, Wang L W, Zhao J L. 2016. Monitoring shallow resistivity changes prior to the 12 May 2008 M8.0 Wenchuan earthquake on the Longmen Shan tectonic zone, China[J]. Tectonophysics, 675: 244-257. doi: 10.1016/j.tecto.2016.03.006
Ma Q Z. 2002. The boundary element method for 3-D DC resistivity modeling in layered earth[J]. Geophysics, 67(2): 610-617. doi: 10.1190/1.1468622
Roinila T, Helin T, Vilkko M, Suntio T, Koivisto H. 2009. Circular correlation based identification of switching power converter with uncertainty analysis using fuzzy density approach[J]. Simul Modell Pract Theory, 17(6): 1043-1058. doi: 10.1016/j.simpat.2009.03.006
Zhao J L, Wang L W, Qian J D. 2011. Research on geo-electrical resistivity observation system specially used for earthquake monitoring in China[J]. Earthquake Science, 24(6): 497-511. doi: 10.1007/s11589-011-0812-4
Ziolkowski A, Hobbs B A, Wright D. 2007. Multitransient electromagnetic demonstration survey in France[J]. Geophysics, 72(4): F197-F207. doi: 10.1190/1.2735802
计量
- 文章访问数: 757
- HTML全文浏览量: 373
- PDF下载量: 23
