Characteristics of subway stray current in geoelectrical resistivity observation
-
摘要: 在对地铁杂散电流产生机理的讨论的基础上，定量计算了地铁运行时杂散电流在地电阻率观测中所呈现的影响幅度，结果表明地铁杂散电流传播范围可以达到几十甚至上百千米。通过对城市周边的北京通州、天津青光、宝坻、塘沽、江苏江宁和辽宁新城子等六个地电阻率观测台站观测到的该类信号的研究，分析了其幅值、主要频率范围以及空间分布特征，结果表明该类信号的幅值从几mV至几十mV不等，与源距关系密切，周期主要集中在50—200 s范围内，在分析地震异常前兆信号时该类信号可使信噪比降低10—30 dB左右，其相对方差最大超出标准20倍左右。根据地铁运行时杂散电流传播的特征，本文提出了几种压制该类信号的措施，为识别地震前地电阻率异常信息及排除噪声提供依据。Abstract: The geoelectrical resistivity observation is one of the most important methods in the study of earthquake precursor, and the characteristics of the stray current during subway operation in the geoelectrical resistivity observation are useful for analyzing the geoelectric observation data and the anomaly variation before earthquakes. Based on the study on the generation mechanism, the quantitatively results of the influence of stray current on geoelectrical resistivity observation were given in this paper. The result shows that the effect distance can reach dozens or even a hundred kilometers. Through monitoring the stray current in some geoelectrical resistivity observation stations around cities, such as Tongzhou station in Beijing, Qingguang, Baodi and Tanggu stations in Tianjin, Jiangning station in Jiangsu and Xinchengzi station in Liaoning, the amplitude and frequency range of the stray current signal as well as its spatial distribution characteristics were analyzed. It shows that the amplitude of the signal ranges from several to tens millivolt, and the period range is mainly from 50 s to 200 s. Due to the influence of stray current, the signal-to-noise ratio is decreased by 10 to 30 dB, and the accuracy of the observation data is more than 20 times worse than the specified requirement. According to the characteristics of the stray current propagation, several methods which can be used to suppress the signal were proposed in this paper, and it will provide a foundation for the identification of anomaly information from the observation data and anti-interference technology study.
-
引言
地震波传播特性与地震动输入是重大工程地震响应机制与灾变机理和抗震设计理论研究的基础，由于涉及地震学与工程学的交叉，这一基础研究滞后于地震响应分析，因此一直是重大工程抗震设计的瓶颈。对于局部不规则场地，其地震波传播特性极其复杂，表现为地形效应十分突出（郭明珠等，2013；高玉峰等，2021）。地震波传播至局部不规则场地后由于散射而导致场地地震动发生放大、衰减以及空间变化，而这种地形效应往往加剧桥梁、隧道和大坝等工程的震害。以2008年汶川地震为例，该地震由于发生在我国西部山区，相较于唐山地震这类平原地区的地震，其震害具有非常显著的差异，该地震造成6 140座桥梁、156条隧道和1 996座水坝受损（杜修力等，2008；景立平等，2009），可见局部地区地形条件的影响不容忽视。
地震波地形效应最早始于对实际地震动记录的分析，在1971年美国圣弗尔南多（San Fernando）地震中，Trifunac和Hudson （1971）发现帕克伊马（Pacoima）大坝上的台站观测数据出现异常的地震动放大现象。为了进一步研究地形效应，世界上很多地震台阵得以布设。1991年我国台湾省翡翠河谷上布置了六台强震仪，1992年9月花莲ML5.5地震触发了此台阵，获得了沿地形分布的地面运动加速度记录，结果显示河谷两侧的地震动差异明显，表现出了明显的地形效应（Huang，Chiu，1995）。Spudich等（1996）在研究加州北岭（Northridge）地震记录时发现，相对于山脚处，山顶出现了明显的地震动放大效应。Parolai等（2004）通过对比分析1999年土耳其伊兹米特（Izmit）地震的17条观测数据，发现复杂的地形条件下会存在高达5倍的地震动放大效应。基于2008年汶川地震期间自贡西山地形影响台阵的观测数据，王海云和谢礼立（2010）详细分析了地形及地质条件对地震动的影响，其结果显示山体地形放大效应随高程的增加而愈加明显。此后，杨宇等（2011）和唐晖等（2012）也针对自贡台阵观测数据进行了深入研究。王文才等（2020）针对2018年陕西宁强MS5.3地震的地震动数据进行分析，其结果显示了山体对地震中频段的地震动放大效应。
为了揭示地形效应的机理，关于地形对地震波的散射和衍射的数值方法和解析研究得以大量开展，数值方法主要包括有限元法（廖振鹏，刘晶波，1992；章小龙等，2017；孙纬宇等，2019）、有限差分法（Boore，1972）、边界元法（林皋，关飞，1990；巴振宁等，2017；Liu et al，2018；何颖等，2019）、谱元法（贺春晖等，2017；于彦彦等，2017；刘中宪等，2020）等。虽然数值方法应用灵活，但是其结果的准确性需要解析研究的验证。解析研究主要是波函数展开法，该方法不仅可以揭示波散射问题的物理本质，还可以检验各种数值方法的精度，近年来得到广泛关注（Trifunac，1973；Yuan，Liao，1994；Tsaur，Chang，2009；Gao et al，2012；Zhang et al，2017，2019；张宁等，2021）。
V形河谷作为常见的河谷形态，由于地震作用下河谷底部会产生应力奇异问题，其地形效应的理论解一直未得到有效解决，直到最近才由Tsaur教授带领的研究小组取得突破，分别得到了浅河谷（Tsaur，Chang，2008）和深河谷（Tsaur et al，2010）的地震SH波散射问题解答。Zhang等（2012a）为了应用方便，给出了不同深度河谷地震波散射的统一解答，随后Zhang等（2012b）研究了非对称V形河谷的地震SH波放大效应，并将此问题拓展至近源地震动的情况（Gao，Zhang，2013）。针对河谷场地上的地震动地形观测台阵，目前仅我国台湾省翡翠河谷上建立的六台强震仪在1992年花莲地震中记录到了相应的数据（Huang，Chiu，1995），且由于此河谷呈现典型的V形形态，本文将运用V形河谷场地地震SH波散射问题的频域解答，模拟其监测点位置处的地震动，进而从理论上揭示地形对地震动响应的影响规律。
1. 翡翠河谷地形及台阵记录
翡翠河谷位于我国台湾省花莲县，其上建有翡翠大坝，是为台北地区450万人供水的水源工程。在1977年到1990年间，科研人员在翡翠大坝附近分别布设了三个SMA-1型强震仪，主要用于观测河谷地形对翡翠大坝附近地震动在空间分布上的影响规律。但是由于台站距离翡翠大坝较近，其观测结果会受到大坝的影响，Huang和Chiu （1995）又于1991年在距离大坝300 m处沿翡翠河谷横断面布置了六台强震仪（SC1—SC6，高程分别为170，120，70，70，120，160 m），台站的具体位置如图1所示。
![]()
图 1 翡翠河谷地震动观测台阵剖面及台站分布示意图（Huang，Chiu，1995）Figure 1. Definition sketch for the cross section of Feitsui canyon and the location of six stations （Huang，Chiu，1995）1992年9月的花莲ML5.5地震触发了此台阵，由此获得了沿地形分布的地面运动加速度记录，此次地震的震源深度为50 km，震中位于监测台阵南偏西方向距台阵130 km处，震中至翡翠大坝的方向基本平行于台阵观测面。Huang和Chiu （1995）根据实测地震动分离出了垂直于监测剖面方向的地震动分量，即二维SH波地震动记录，为后续地震动模拟提供了便利。台阵观测数据显示出显著的地震动地形效应，各点峰值加速度（peak ground acceleration，简写为PGA）记录见表1 （Huang，Chiu，1995）。由表1可以看出，河谷迎波面相对于背波面表现出较强的地震动放大现象，迎波侧PGA最大可达9.7 cm/s2，而背波侧最大只有6.0 cm/s2；同时可以看出越靠近河谷底部，记录到的地震动PGA越小。为了分析这一现象产生的原因，本文将从V形河谷周围地震波传播解析理论出发得到各点的地震动模拟结果，以揭示其地震响应规律。
表 1 1992年9月花莲地震中SC1—SC6台站记录到的地震动PGA （Huang，Chiu，1995）Table 1. PGA values at stations SC1−SC6 during Hualien earthquake in September 1992 （Huang，Chiu，1995）台站名 PGA/（cm·s−2） 台站名 PGA/（cm·s−2） 台站名 PGA/（cm·s−2） SC1 9.7 SC2 6.5 SC3 5.2 SC4 3.6 SC5 4.5 SC6 6.0 2. 理论模拟
2.1 频域内地震动场推导
由翡翠河谷剖面图可以看到，河谷形状大致呈V形，地震波自剖面左侧斜入射。本文将首先推导线源SH波入射下非对称V形河谷地震波传播理论，得到频域内考虑地形效应的翡翠河谷各点的地震动传递函数。由于翡翠河谷属于深度较浅的河谷，因此本文将针对浅的非对称V形河谷地震波传播理论进行推导。河谷的二维简化模型如图2所示，此图表示半空间中一个非对称V形河谷，深度为d，半宽为b1 （左侧）和b2 （右侧）；模型的介质假设为弹性、各向同性、均质；介质的剪切模量为G和剪切波速为vS。入射波假设为简谐的线源柱面SH波，圆频率为ω，位移在z方向（出平面）。在河谷底部建立局部坐标系（r1，θ1），在河谷中点建立整体坐标系（r，θ），采用半径为（b1＋b2）/2的半圆形辅助边界把整个空间分成两个子区域①和②。震源位置在整体极坐标系（r，θ）中为（r0，θ0），对于本研究的问题，需使用镜像法来考虑相对于水平地面对称的两个波源，对于波源及其镜像分别使用两个极坐标系（rf，θf）和（
${{r'}_{\rm{f}}},\;\;{{\theta '} _{\rm{f}}}$ ）。通过四个坐标系的建立，引入半圆形辅助边界进行区域分解，得到区域①和区域②的运动方程（Sanchez-Sesma，1985；Gao，Zhang，2013）：
$$ \frac{{{\partial ^2}{u_1}}}{{\partial r_1^2}} {\text{＋}} \frac{1}{{{r_1}}}\frac{{\partial {u_1}}}{{\partial {r_1}}} {\text{＋}} \frac{1}{{r_1^2}}\frac{{{\partial ^2}{u_1}}}{{\partial \theta _1^2}} {\text{＋}} {k^2}{u_1} {\text{＝}} 0 {\text{，}} $$ (1) $$ \frac{{{\partial ^2}{u_2}}}{{\partial {r^2}}} {\text{＋}} \frac{1}{r}\frac{{\partial {u_2}}}{{\partial r}} {\text{＋}} \frac{1}{{{r^2}}}\frac{{{\partial ^2}{u_2}}}{{\partial {\theta ^2}}} {\text{＋}} {k^2}{u_2} {\text{＝}} - \frac{{{\text{δ}} （r {\text{－}} {r_0}）{\text{δ}} （\theta {\text{－}} {\theta _0}）}}{r}{\text{，}} $$ (2) 式中，
${\text{δ}} （ \cdot ）$ 为狄拉克δ函数，k＝ω/vS为剪切波数。对于稳态出平面问题，u1和u2分别表示区域①和区域②在频域内的位移。引入镜像法将入射波源镜像至（
${{r}'}_{\rm{f}}{\text{，}}{{\theta }'} _{\rm{f}}$ ）坐标系，通过波函数展开法可以得到区域①和区域②内含未知系数的位移表达式（Zhang et al，2012b）：$$ {u_1}({r_1}{\text{，}}{\theta _1}) {\text{＝}} \sum\limits_{n {\text{＝}} 0}^\infty {{A_n}{{\text{J}}_{n\upsilon }}(k{r_1})\cos \left[ {n\upsilon （ {{\theta _1} {\text{＋}} {\beta _1}} ）} \right]}{\text{，}} $$ (3) $$ {u_2}（r{\text{，}}\theta ） {\text{＝}} {u^{\rm{f}}} {\text{＋}} \sum\limits_{n {\text{＝}} 0}^\infty {{B_n}{\text{H}}_{2n}^{（1）}（kr）\cos （2n\theta ）} {\text{＋}} \sum\limits_{n {\text{＝}} 0}^\infty {{C_n}{\text{H}}_{2n {\text{＋}} 1}^{（1）}（kr）\sin [（2n + 1）\theta ]} {\text{，}} $$ (4) 式中：
${{\rm{J}}_n}（\cdot ）$ 为n阶第一类贝塞尔函数，${\rm{H}}_n^{（1）}（ \cdot ）$ 为n阶第一类汉克尔函数；$\;{\beta _1} {\text{＝}} \pi {\text{－}}\arctan （{b_1}/d）$ ，$\;{\beta _2} {\text{＝}} \pi {\text{－}}\arctan （{b_2}/d ）$ ，$ \upsilon $ ＝ π/（β1＋β2）；${u^{\rm{f}}}{\text{＝}}[ {{\rm{H}}_0^{(1)}(k{r_{\rm{f}}})＋{\rm{H}}_0^{(1)}(k{{r'}_{\rm{f}}})} ]/{\rm{H}}_0^{（1）}（k{r_0}）$ ，表示自由场。上述位移表达式需要满足河谷表面的应力自由边界条件和水平地表的应力自由边界条件：
$$ \tau _{{\theta _1}{{\textit{z}}_1}}^{（1）} {\text{＝}} \frac{G}{{{r_1}}}\frac{{\partial {u_1}}}{{\partial {\theta _1}}} {\text{＝}} 0 \quad\quad {\theta _1} {\text{＝}} - {\beta _1}{\text{，}}{\beta _2} {\text{，}} $$ (5) $$ \tau _{\theta {\textit{z}}}^{（2）} {\text{＝}} \frac{G}{r}\frac{{\partial {u_2}}}{{\partial \theta }} {\text{＝}} 0 \quad\quad \theta {\text{＝}} {\text{±}} \frac{\pi }{2} {\text{，}} r {\text{≥}} \frac{{{b_{\text{1}}}{\text{ + }}{b_{\text{2}}}}}{2} {\text{，}} $$ (6) 同时满足两个子区域之间的位移连续和应力连续条件：
$$ {u_1}（r{\text{，}}\theta ） {\text{＝}} {u_2}（r{\text{，}}\theta） \quad\quad r {\text{＝}} \frac{{{b_{\text{1}}}{\text{＋}}{b_{\text{2}}}}}{2} {\text{，}} - \frac{\pi }{2} {\text{≤}} \theta {\text{≤}} \frac{\pi }{2} {\text{，}} $$ (7) $$ \tau _{r{\textit{z}}}^{（1）}（r{\text{，}}\theta ） {\text{＝}} \tau _{r{\textit{z}}}^{（2）}（r{\text{，}}\theta ） \quad\quad r {\text{＝}} \frac{{{b_{\text{1}}}{\text{＋}}{b_{\text{2}}}}}{2} {\text{，}} - \frac{\pi }{2} {\text{≤}} \theta {\text{≤}} \frac{\pi }{2} {\text{.}} $$ (8) 借助区域匹配策略和Graf加法公式，通过恰当的求解方案，最终可以得到位移场各未知系数An，Bn和Cn的解答，进而各区域的位移即可根据方程（3）和（4）精确得到。具体求解方案可参考Gao和Zhang （2013）。
2.2 地震动时程模拟
地震动记录多数是在地表面实测的加速度时程，测点一般位于不受任何结构影响的自由场地。理论上，对于地震SH波，自由场地加速度幅值是地下基岩地震动幅值的两倍，因此，为了得到结构抗震所需的输入地震波，对于平坦场地，将地表地震动记录的幅值减半作为输入的基岩地震动；而对于河谷场地，需要利用我们提出的地震动传播理论得到地形放大因子（相对于平坦场地）来获取场地地震动。
为了在给定地震动加速度时程输入的情况下得到河谷场地的空间变化地震动，需建立考虑地形放大效应的河谷场地地震动时程合成方法。这一方法建立的关键在于理解并应用河谷引起的地形放大因子，说明如下：V形河谷的波函数级数解以单位幅值稳态SH波作为入射波，如果河谷不存在，不同位置的水平地表位移幅值恒相等，为自由场|uf|；而河谷地形效应导致不同位置的地表位移幅值在自由场|uf|上下波动，也就是说，在某一频率下，如果地表某一点的位移幅值大于自由场|uf|，那么这一点的地面运动相对于地表自由场放大，反之，如果小于自由场|uf|，地面运动减小。根据这一原理，将这些波函数级数解给出的河谷地表位移除以自由场|uf|，即可得到河谷地形的地震动放大因子，该放大因子为复数（包含幅值放大因子和相位调整因子）。因此如果得到了上述河谷地形的地震动放大因子随频率变化的函数，则可根据基岩地震动得到场地某点的实际地震动，具体步骤参见代登辉（2019）。
2.3 翡翠河谷测点地震动模拟
由于监测台阵距震中130 km，震源深度为50 km。因此震源相对于河谷的位置为（x0，y0）＝（−130 km，50 km）。根据Huang和Chiu （1995）的数据，翡翠河谷的模型参数为：深度d＝350 m，宽度b1 ≈ 800 m，b2 ≈ 700 m，介质密度ρ＝2.67 g/cm3，剪切波速vS＝1500 m/s。
运用频域内地震动场的推导公式，可以得到河谷场地目标位置的频域地震动放大因子，地震动放大因子为复数，包含幅值和相位的信息。图3给出了台站SC1—SC6所在位置的地震动放大因子随频率变化的曲线，计算的频率范围为0—10 Hz。从图中可以明显地看出，在0—10 Hz内，SC1—SC3台站处的地震动放大因子最大可达1.5左右，而SC4—SC6台站处的地震动放大因子均小于1。尽管在某些频率（如5 Hz左右）下，迎波侧的地震动放大因子会小于背波侧，但是迎波侧整体表现出相对于背波侧的地形放大现象，这也与实测地震动结果相吻合。
![]()
图 3 翡翠河谷上SC1—SC6台站所在位置的地震动放大因子随频率f的变化Figure 3. Variation of ground motion amplification factors as a function of frequency at the stations SC1−SC6 in Feitsui canyon对于翡翠河谷台阵记录到的1992年9月的花莲地震，其主要频率在5 Hz以内，且主要集中在2—4 Hz，持时约30 s。由于缺乏相应的地震记录数据，为了更可靠地模拟该场地的地震响应，本文通过对Huang和Chiu （1995）一文中的数据进行读点，并依据实测数据的频谱结构适当调整，构造出符合上述频谱特征的人工地震动时程，并以SC3点作为地震动输入参考点，计算其它各点的地震响应。地震动输入参考点SC3的地震动时程如图4a所示，持时为30 s，PGA为5.2 cm/s2，相应的傅里叶谱如图4b所示，可见所构造的地震波与实测地震动波形基本一致，主要频率集中在2—4 Hz。
![]()
图 4 输入参考点SC3台站的地震动时程（a）和相应的傅里叶谱（b）Figure 4. Ground motion history （a） and corresponding Fourier spectrum （b） at the input reference point SC3基于输入地震动和频域地震动放大因子，可以进一步得到河谷各点的地震动响应时程，关键步骤解释如下：
1） 将输入地震动加速度时程记录进行两次数值积分，得到位移时程uinp（t），然后将位移时程进行傅里叶变换，得到其傅里叶谱
$$ {u_{{\rm{inp}}}}（ \omega ） {\text{＝}} \int_{ - \infty }^{ {\text{＋}} \infty } {{u_{{\rm{inp}}}}（t）{{\rm{e}}^{ - {\rm{i}}\omega t}}{\rm{d}}t} {\text{；}} $$ (9) 2） 根据频域内地震动场的推导，求解不同频率简谐振动下河谷场地各个位置的地形放大因子，得到河谷场地目标位置的频域放大传递函数A（ω）；
3） 根据频域放大因子调整目标位置地震动的傅里叶谱
$$ {u_{{\rm{outp}}}}（\omega ） {\text{＝}} A（\omega ）{u_{{\rm{inp}}}}（\omega ）{\text{；}} $$ (10) 4） 根据调整后的傅里叶谱，通过傅里叶逆变换来求解考虑地形效应的地震动时程，即
$$ {u_{{\rm{outp}}}}（t） {\text{＝}} \frac{1}{{2\pi }}\int_{ - \infty }^{ {\text{＋}} \infty } {{u_{{\rm{outp}}}}（\omega ）{{\rm{e}}^{{\rm{i}}\omega t}}{\rm{d}}t} {\text{，}} $$ (11) 将位移时程uoutp（t）进行一次数值微分后得到速度时程，进行两次数值微分后得到加速度时程。
通过计算得到的地震动放大因子和地震动时程模拟方法得到各测点位置处的地震动加速度时程，如图5所示。表2给出了模拟得到的地震动PGA与监测结果的对比。从图5和表2可以看出，SC1—SC3点的地震动PGA明显大于SC4—SC6点，这说明模拟得到的各测点地震动地形分布规律与观测结果一致，即翡翠河谷迎波侧的地震动明显大于背波侧。这是因为，河谷的存在阻挡了地震波的传播，这使得地震波能量大部分集中在河谷迎波面，河谷背波面的地震运动主要是源于河谷引起的散射波的作用。同时可以看出：在河谷迎波面，SC1点的地震动PGA大于SC2点，SC3点最小；在河谷背波面，SC6点的地震动PGA大于SC5点，SC4点最小。这与观测结果相一致，同时可以得出：对于河谷同一坡面，随着所处高程的增加，其地震动也相应增大。从表2中可以看出，以SC3点作为地震动输入参考点，模拟得到的SC2和SC4点的地震动PGA与实测结果最为接近，SC1，SC5和SC6点的模拟结果与实测结果有些许差别。这是因为地震动不只受地形因素的影响，还会受很多其它因素的干扰，如测点所处位置的地质差异、河谷底部沉积层的影响等，所以离地震动输入参考点SC3越远，地震动受其它因素的干扰越大，其结果差别越明显。但是，这种差异不影响结果的指导性，从整体上看模拟结果能够较好地反映河谷场地表现出的地震动地形效应。
![]()
图 5 翡翠河谷上SC1—SC6台站的模拟加速度时程Figure 5. Response acceleration time histories at the stations SC1−SC6 in Feitsui canyon表 2 SC1−SC6台站的模拟地震动PGA与监测结果对比Table 2. Comparisons of PGA values at the stations SC1−SC6 between the records and simulation台站名 监测PGA
/（cm·s−2）模拟PGA
/（cm·s−2）台站名 监测PGA
/（cm·s−2）模拟PGA
/（cm·s−2）台站名 监测PGA
/（cm·s−2）模拟PGA
/（cm·s−2）SC1 9.7 7.8 SC3 5.2 5.2 SC5 4.5 3.6 SC2 6.5 6.3 SC4 3.6 3.1 SC6 6.0 3.9 图6给出了计算得到的各监测点位置的地震动傅里叶谱幅值，可以看出地形效应对地震波的主要频率范围影响较小，但对其傅里叶谱幅值的影响较大。SC1—SC3点的傅里叶谱幅值基本能反映其时域情况，SC1点的傅里叶谱幅值整体大于SC2点和SC3点，SC4—SC6点的傅里叶谱幅值基本相同，因此其在时域内的加速度时程也相差较小，三者的地震动PGA较为接近。由频谱特征分析可知：SC1—SC3点的地震波主要频率集中在2—4 Hz，在该区间内傅里叶频谱幅值明显大于背波侧SC4—SC6点的傅里叶频谱幅值；而SC4—SC6点的地震波主要频率集中在3—5 Hz；在4—5 Hz的频率范围内，六个观测点的傅里叶谱幅值几乎保持不变。这是因为在2—4 Hz内，SC1—SC3点相对于SC4—SC6点产生较大的地形放大，而在4—5 Hz内，六个观测点的地形放大因子差别不大，这一规律同样可以在图3中看到。
![]()
图 6 模拟得到的翡翠河谷上SC1—SC6各点加速度傅里叶谱幅值Figure 6. Response Fourier amplitudes of acceleration at the stations SC1−SC6 in Feitsui canyon3. 讨论与结论
翡翠河谷地形影响地震动观测台阵在1992年花莲地震中获得了相应的地形影响地震动记录，通过分析河谷各监测点的地震动可知，河谷迎波面相对于背波面表现出明显的地震动放大现象。为了分析此现象产生的原因，本文利用线源SH波入射下非对称V形河谷地震波传播解析理论，通过构造频段接近实际观测数据的地震波，并以SC3点作为地震动输入参考点，模拟得到了河谷各监测点位置的地震动。模拟结果进一步揭示了河谷地形效应对地震动影响的规律：河谷地形使得地震波在传播过程中产生散射，导致河谷场地表现出了与平坦场地不同的地震动规律。在河谷迎波侧，散射波与入射波相长干涉，表现出了地震动放大现象；在河谷背波侧，河谷对剪切波产生了过滤作用，表现出了地震动衰减现象。结果进一步表明，河谷同侧地震动表现出了地震动PGA随高程增加而增加的趋势。通过对比实测地震动与模拟结果可知，本文提出的河谷场地地震波传播解析理论能够较好地模拟地形效应引起的地震动差异情况，河谷两侧的地震动差异以及所处高程对地震动的影响均能够通过本解析理论得到较好的模拟。本解析理论不同于其它数值方法，能够实现河谷场地全域全过程的地震动模拟，同时本文提出的研究思路实现了解析理论的工程应用，本文中考虑地形效应的地震动模拟方法可为V形河谷场地重大工程的抗震设计提供精确的地震动输入，所得结果对地震区划的精细划分有一定的指导意义。
-
![]()
图 1 地铁单车运行时电流示意图（a）及简化等效电路图（b）
Figure 1. The schematic diagram of current （a） and equivalent circuit （b） for subway single-train operation
![]()
图 3 各观测台南北（左）、东西（右）测道记录的地铁杂散电流信号
（a） 北京通州台 （2014-04-21 13:00开始）；（b） 天津青光台 （2014-05-16 16:00开始）；（c） 天津塘沽台（2014-05-13 17:00开始）；（d） 天津宝坻台 （2014-05-12 11:00开始）
Figure 3. The stray current signals recorded from north-south （left） and east-west （right） channel at different observations
（a） Tongzhou station in Beijing （from 13:00，April 21，2014）；（b） Qingguang station in Tianjin （from 16:00，May 16，2014）；（c） Tanggu station in Tianjin （from 17:00，May 13，2014）；（d） Baodi station in Tianjin （from 11:00，May 12，2014）
![]()
图 3 各观测台南北（左）、东西（右）测道记录的地铁杂散电流信号
（e） 江苏江宁台 （2014-10-13 15:00开始）；（f） 辽宁新城子台 （2014-07-30 17:00开始）
Figure 3. The stray current signals recorded from north-south （left） and east-west （right） channe at different observations
（e） Jiangning station in Jiangsu （from 15:00，October 13，2014）；（f） Xinchengzi station in Liaoning （from 17:00，July 30，2014）
![]()
图 4 地铁运行时段和停运时段功率谱密度PSD对比
（a） 北京通州台；（b） 天津青光台；（c） 天津塘沽台；（d） 天津宝坻台；（e） 江苏江宁台；（f） 辽宁新城子台
Figure 4. Comparison of PSDs in subway operation with those in non-operation time period
（a） Tongzhou station in Beijing；（b） Qingguang station in Tianjin；（c） Tanggu station in Tianjin；（d） Baodi station in Tianjin；（e） Jiangning station in Jiangsu；（f） Xinchengzi station in Liaoning
![]()
图 5 天津青光台地电阻率小时观测数据曲线（2019年12月30日至31日）
（a） 南北方向地电阻率ρNS；（b） 南北方向测量均方差σNS
Figure 5. The curve of hourly geoelectrical resistivity observation data of Qingguang （30th to 31th，Dec 2019）
（a） Geoelectrical resistivity ρNS of NS direction；（b） RMS value σNS of NS direction
![]()
图 6 江苏江宁台地电阻率小时观测数据曲线（2017年1月1日至2日）
（a） 东西方向地电阻率ρEW；（b） 东西方向测量均方差σEW
Figure 6. The curve of hourly geoelectrical resistivity observation data of Jiangning （1st to 2nd，Jan 2017）
（a） Geoelectrical resistivity ρEW of EW direction；（b） RMS value σEW of EW direction
![]()
图 7 天津青光台日均值数据曲线（2017年1月1日至30日）
（a） 东西方向地电阻率ρEW；（b） 东西方向相对均方差δEW
Figure 7. The daily mean value of geoelectrical resistivity observation data （1st to 30th，January 2017）
（a） Geoelectrical resistivity ρEW of EW direction；（b） Relative RMS value δEW of EW direction
![]()
图 8 江苏江宁台日均值观测数据曲线（2017年1月1日至30日）
（a） 东西方向地电阻率ρEW；（b） 东西方向相对均方差δEW
Figure 8. The daily mean value of geoelectrical resistivity observation data （1st to 30th ，Jan 2017）
（a） Geoelectrical resistivity ρEW of EW direction；（b） Relative RMS value δEW of EW direction
![]()
图 9 江苏江宁台地电阻率交、直流小时观测数据曲线（2016年9月29日至30日）
（a） 南北方向地电阻率；（b） 东西方向地电阻率；（c） 南北方向测量均方差；（d） 东西方向测量均方差
Figure 9. Comparison of AC & DC hourly geoelectrical resistivity observation data of Jiangning station
（a） Geoelectrical resistivity of NS direction；（b） Geoelectrical resistivity of EW direction；（c） RMS value of NS direction；（d） RMS value of EW direction
表 1 地铁杂散电流导致的地电阻率测量相对误差δ随距离的变化（P1＝75，P2＝1，C＝3.33）
Table 1 The influence of stray current on geoelectrical resistivity varies with distance （P1＝75，P2＝1，C＝3.33）
测区至地铁
距离r1/kmδ dAB＝1000 m dAB＝500 m dAB＝200 m 10 8.530% 2.132% 0.341% 20 2.132% 0.533% 0.085% 30 0.948% 0.237% 0.038% 40 0.533% 0.133% 0.021% 50 0.341% 0.085% 0.014% 60 0.237% 0.059% 0.009% 70 0.174% 0.044% 0.007% 80 0.133% 0.033% 0.005% 90 0.105% 0.026% 0.004% 100 0.085% 0.021% 0.003% 
下载: 导出CSV
表 2 各个台站观测装置及与地铁最近距离
Table 2 The configuration in geoelectrical resistivity observation station and its distance from the subway
序号 台站名称 与地铁最近距离/km 观测装置 测道方向 供电极距/m 测量极距/m 装置系数/m 1 北京通州台 11.0 南北、东西 1 760 320 7 351 2 天津青光台 7.6 南北、东西 1 000 316 2 237 3 天津塘沽台 城区线：40.0
九号线：7.0南北、东西
南北、东西南北：1 500
东西：1 000南北：500
东西：3003 141
2 3824 天津宝坻台 52.0 南北、东西 1 000 200 3 769 5 江苏江宁台 城区线：30.0
机场线：3.0南北、东西
南北、东西1 000
1 000300
3002 382
2 3826 辽宁新城子台 19.4 南北、东西 1 000 300 2 382
下载: 导出CSV
表 3 不同台站地铁运行时杂散电流不同周期功率占比
Table 3 The distribution of power in different periods at different stations during subway operation
周期范围/s 北京通州台 天津青光台 天津塘沽台 天津宝坻台 江苏江宁台 辽宁新城子台 50—100 10.17% 19.02% 8.74% 14.78% 33.34% 24.26% 100—150 29.56% 60.50% 11.66% 54.55% 6.77% 57.57% 150—200 56.87% 18.13% 76.35% 27.00% 9.13% 15.67% 200—250 2.38% 1.68% 0.85% 1.12% 39.28% 2.10% 250—300 0.17% 0.20% 1.81% 0.57% 0.77% 0.00% 300—350 0.83% 0.47% 0.59% 1.97% 10.71% 0.40%
下载: 导出CSV
表 4 六个台站记录的地铁杂散电流信号幅度及测量信噪比
Table 4 The stray current signal amplitude and the SNR for the six stations
台站名称 测道
方向人工电位差/mV 噪声幅度/mV 测量信噪比/dB 信噪比降低
/dB停运时段 运行时段 停运时段 运行时段 北京通州 南北 7 0.14 1.33 34.0 14.4 19.6 东西 7 0.17 1.16 32.3 15.6 16.7 天津青光 南北 9 0.09 4.96 40.0 5.2 34.8 东西 9 0.36 2.25 28.0 12.0 15.9 天津塘沽 南北 7 0.06 1.94 41.3 11.1 30.2 东西 10 0.40 0.79 28.0 22.0 5.9 天津宝坻 南北 41 0.07 0.48 55.4 38.6 16.7 东西 41 0.06 0.25 56.7 44.3 12.4 江苏江宁 南北 93 0.19 3.60 53.8 28.2 25.6 东西 73 0.56 13.86 42.3 14.4 27.9 辽宁新城子 南北 26 0.13 1.02 34.0 14.4 19.6 东西 26 0.10 0.46 32.3 15.6 16.7
下载: 导出CSV
表 5 六个台站的地电阻率变幅与相对方差
Table 5 The maximum variation and relative RMS value of geoelectrical resistivity for six stations
台站 电阻率最大变幅/Ω·m 相对均方差 南北 东西 南北 东西 北京通州 1.83 3.99 0.65% 1.75% 天津青光 1.29 0.48 3.73% 1.72% 天津塘沽 1.72 0.59 7.57% 2.55% 天津宝坻 0.32 0.22 0.36% 0.20% 江苏江宁 1.50 5.22 1.13% 5.37% 辽宁新城子 0.09 0.05 0.15% 0.08%
下载: 导出CSV
表 6 六个台站的理论避让距离（P1＝75，P2＝1，δ＝0.3%）
Table 6 The minimum distance between observation station and subway for the six stations （P1＝75，P2＝1，δ＝0.3%）
台站名称 供电极距/m 测量极距/m C＝dAB/dMN 最小避让距离/km 北京通州台 1 760 320 5.50 97 天津青光台 1 000 316 3.16 53 天津塘沽台 1 500 500 3.00 79 天津宝坻台 1 000 200 5.00 55 江苏江宁台 1 000 300 3.33 53 辽宁新城子 1 000 300 3.33 53 注：P1为地电阻率测量时噪声与信号之比，P2为地铁到台站电阻率与台站电阻率之比，δ为相对均方差。
下载: 导出CSV
-
董亮,姚知林,葛彩刚,石超杰,陈金泽. 2021. 地铁杂散电流干扰下管地电位波动特征的傅里叶分析[J]. 表面技术,50(2):294–303. doi: 10.16490/j.cnki.issn.1001-3660.2021.02.031 Dong L,Yao Z L,Ge C G,Shi C J,Chen J Z. 2021. Fourier analysis of the fluctuation characteristics of pipe-to-soil potential under metro stray current interference[J]. Surface Technology,50(2):294–303 (in Chinese).
杜学彬,马占虎,叶青,谭大诚,陈军营. 2006. 与强地震有关的视电阻率各向异性变化[J]. 地球物理学进展,21(1):93–100. doi: 10.3969/j.issn.1004-2903.2006.01.015 Du X B,Ma Z H,Ye Q,Tan D C,Chen J Y. 2006. Anisotropic changes in apparent resistivity associated with strong earthquakes[J]. Progress in Geophysics,21(1):93–100 (in Chinese).
桂燮泰,关华平,戴经安. 1989. 唐山、松潘地震前视电阻率短临异常图象重现性[J]. 西北地震学报,11(4):71–75. Gui X T,Guan H P,Dai J A. 1989. The short-term and immediate anomalous pattern recurrences of the apparent resistivity before the Tangshan and Songpan earthquakes of 1976[J]. Northwestern Seismological Journal,11(4):71–75 (in Chinese).
胡先茂. 2011. 城市轨道交通对埋地管线影响范围的分析[J]. 现代城市轨道交通,(3):74–76. doi: 10.3969/j.issn.1672-7533.2011.03.024 Hu X M. 2011. Analysis on influence scope of urban rail transit on buried pipelines[J]. Modern Urban Transit,(3):74–76 (in Chinese).
李金铭. 2005. 地电场与电法勘探[M]. 北京: 地质出版社: 60–67. Li J M. 2005. Geoelectric Field and Electrical Exploration[M]. Beijing: Geological Publishing House: 60–67 （in Chinese）.
李雷,刘楠康,曾文,邵华锋. 2019. 轨道交通杂散电流在复杂地质条件下的分布[J]. 广东电力,32(8):133–140. doi: 10.3969/j.issn.1007-290X.2019.008.017 Li L,Liu N K,Zeng W,Shao H F. 2019. Distribution of stray current in rail transit under complex geological conditions[J]. Guangdong Electric Power,32(8):133–140 (in Chinese).
林江,唐华,于海学. 2002. 地铁迷流腐蚀及其防护技术[J]. 建筑材料学报,5(1):72–76. doi: 10.3969/j.issn.1007-9629.2002.01.014 Lin J,Tang H,Yu H X. 2002. Protection of stray current corrosion in metro[J]. Journal of Building Materials,5(1):72–76 (in Chinese).
刘国兴. 2005. 电法勘探原理与方法[M]. 北京: 地质出版社: 7–13. Liu G X. 2005. The Principle and Method of Electrical Exploration[M]. Beijing: Geological Publishing House: 7–13 （in Chinese）.
刘文权. 2014. 城市地铁杂散电流对埋地输油管道的危害[J]. 全面腐蚀控制,28(11):29–32. doi: 10.13726/j.cnki.11-2706/tq.2014.11.029.04 Liu W Q. 2014. Interference corrosion hazards of subway stray currents on buried jet fuel metal pipeline[J]. Total Corrosion Control,28(11):29–32 (in Chinese).
马钦忠. 2014. 中外几次重要地震预测与预报结果之启示[J]. 地震学报,36(3):500–513. doi: 10.3969/j.issn.0253-3782.2014.03.015 Ma Q Z. 2014. Enlightment of the success or failure prediction for some large earthquakes at home and abroad[J]. Acta Seismologica Sinica,36(3):500–513 (in Chinese).
毛先进,杨玲英,钱家栋. 2014. 水平层状介质中深埋装置系统地电阻率影响系数特征研究[J]. 地震学报,36(4):678–685. doi: 10.3969/j.issn.0253-3782.2014.04.013 Mao X J,Yang L Y,Qian J D. 2014. Characteristics of the influence coefficient in the cases of deeply-buried configurations for geoelectrical resistivity observation[J]. Acta Seismologica Sinica,36(4):678–685 (in Chinese).
梅进武,林国松. 2017. 多列车运行情况下的地铁杂散电流分析[J]. 电气化铁道,28(4):68–70. Mei J W,Lin G S. 2017. Analysis of metro stray current under multi-train operation[J]. Electric Railway,28(4):68–70 (in Chinese).
聂永安,巴振宁,聂瑶. 2010. 深埋电极的地电阻率观测研究[J]. 地震学报,32(1):33–40. doi: 10.3969/j.issn.0253-3782.2010.01.004 Nie Y A,Ba Z N,Nie Y. 2010. Study on buried electrode resistivity monitoring system[J]. Acta Seismologica Sinica,32(1):33–40 (in Chinese).
钱复业,赵玉林,于谋明,王志贤,刘小伟,常思敏. 1982. 地震前地电阻率的异常变化[J]. 中国科学:B辑,(9):831–839. Qian F Y,Zhao Y L,Yu M M,Wang Z X,Liu X W,Chang S M. 1982. Geoelectrical resistivity anomalies before earthquake[J]. Science in China:Series B,(9):831–839 (in Chinese).
钱家栋. 1993. 与大震孕育过程有关的地电阻率变化研究[J]. 中国地震,9(4):341–350. 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):341–350 (in Chinese).
钱家栋,马钦忠,李劭秾. 2013. 汶川MS8.0地震前成都台NE测线地电阻率异常的进一步研究[J]. 地震学报,35(1):4–17. doi: 10.3969/j.issn.0253-3782.2013.01.002 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).
澹台乐琰,韩肖清,王磊,袁铁江. 2020. 多列车运行下地铁杂散电流建模仿真[J]. 电测与仪表,57(22):7–16. Tantai L Y,Han X Q,Wang L,Yuan T J. 2020. Modeling and simulation of stray current in subway with multi-train operation[J]. Electrical Measurement &Instrumentation,57(22):7–16 (in Chinese).
王崇林,马草原,王智,潘春德,王雨扬. 2007. 地铁直流牵引供电系统杂散电流分析[J]. 城市轨道交通研究,10(3):51–53. Wang C L,Ma C Y,Wang Z,Pan C D,Wang Y Y. 2007. Analysis of stray current in metro DC traction power system[J]. Urban Mass Transit,10(3):51–53 (in Chinese).
王兰炜,张宇,张世中,颜蕊,王子影,张兴国,胡哲. 2015. 我国井下地电阻率观测技术现状分析[J]. 地震地磁观测与研究,36(2):95–102. doi: 10.3969/j.issn.1003-3246.2015.02.018 Wang L W,Zhang Y,Zhang S Z,Yan R,Wang Z Y,Zhang X G,Hu Z. 2015. The status of deep-well geo-electrical resistivity observation in China[J]. Seismological and Geomagnetic Observation and Research,36(2):95–102 (in Chinese).
王兰炜,张宇,张兴国,胡哲,王子影,马小溪. 2019. 地震地电阻率交流观测方法及观测实验[J]. 大地测量与地球动力学,39(7):738–742. Wang L W,Zhang Y,Zhang X G,Hu Z,Wang Z Y,Ma X X. 2019. AC geo-electrical resistivity observation method and experimental observation[J]. Journal of Geodesy and Geodynamics,39(7):738–742 (in Chinese).
王猛. 2005. 直流牵引供电系统钢轨电位与杂散电流分析[J]. 城市轨道交通研究,8(3):24–26. doi: 10.3969/j.issn.1007-869X.2005.03.008 Wang M. 2005. Rail potential and stray current of DC traction power system[J]. Urban Mass Transit,8(3):24–26 (in Chinese).
张世中,石航,王兰炜,胡哲,刘大鹏,魏连生,鞠永. 2013. 地电台站受城市轨道交通干扰的测试分析与抗干扰措施研究[J]. 地震学报,35(1):117–124. doi: 10.3969/j.issn.0253-3782.2013.01.012 Zhang S Z,Shi H,Wang L W,Hu Z,Liu D P,Wei L S,Ju Y. 2013. Test analysis on disturbances caused by urban rail transit at geoelectric stations and measures to reduce its influence[J]. Acta Seismologica Sinica,35(1):117–124 (in Chinese).
张晓东,蒋海昆,李正媛,卢显,安艳茹. 2011. 汶川地震对地震预报工作的一些启示[J]. 工程研究:跨学科视野中的工程,3(4):309–320. Zhang X D,Jiang H K,Li Z Y,Lu X,An Y R. 2011. The revelation of Wenchuan earthquake for earthquake forecast[J]. Journal of Engineering Studies,3(4):309–320 (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 (in Chinese).
中国地震局. 2009. DB/T 33.1-2009地震地电观测方法地电阻率观测第1部分: 单极距观测[S]. 北京: 中国标准出版社: 3. China Earthquake Administration. 2009. DB/T 33.1−2009 The Method of Earthquake-Related Geoelectrical Monitoring: Geoelectrical Resistivity Observation: Part 1: Single Separation Configuration[S]. Beijing: Standards Press of China: 3 （in Chinese）.
朱峰,李嘉成,曾海波,邱日强. 2018. 城市轨道交通轨地过渡电阻对杂散电流分布特性的影响[J]. 高电压技术,44(8):2738–2745. Zhu F,Li J C,Zeng H B,Qiu R Q. 2018. Influence of rail-to-ground resistance of urban transit systems on distribution characteristics of stray current[J]. High Voltage Engineering,44(8):2738–2745 (in Chinese).
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 MS8.0 Wenchuan earthquake on the Longmen Shan tectonic zone,China[J]. Tectonophysics,675:244–257. doi: 10.1016/j.tecto.2016.03.006
-
期刊类型引用(7)
1. 刘中宪，孟思博，张妤，乔云帆，陈龙伟. 考虑建筑群-沉积盆地动力相互作用的建筑群震害评估方法. 地震学报. 2024(01): 129-143 . 
本站查看
2. 王睿，王兰民，周燕国，王刚. 土动力学与岩土地震工程. 土木工程学报. 2024(07): 71-89+105 . 百度学术
3. 赵仕兴，罗麒锐，熊峰，吴启红，夏静，杨姝姮，张敏. 复杂区域地形条件下台地地震动放大效应研究. 振动与冲击. 2024(23): 230-239+304 . 百度学术
4. 卢林彤，李孝雄，庄海洋，吴琪. 宽河谷非均匀工程场地非线性地震反应空间变异特征研究. 地震工程与工程振动. 2023(04): 216-225 . 百度学术
5. 柳国环，陈欣宇，李鑫洋. 半覆水相变V形峡谷场地P波入射下的空间地震动散射频域理论解、验证及特征. 计算力学学报. 2023(05): 764-772 . 百度学术
6. 刘中宪，周健，李程程，孟思博. 基于修正随机有限断层法的沉积盆地地震动模拟. 防灾减灾工程学报. 2023(05): 999-1008 . 百度学术
7. 蒋乐英，廖意辉，王志明，吴勇信，李梅. 河谷差异地震作用下拱式倒虹吸地震响应分析. 长江科学院院报. 2022(12): 117-121+140 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 1447
- HTML全文浏览量: 249
- PDF下载量: 161
- 被引次数: 7
