Cause for abnormal decline of DC apparent resistivity at Tanggu station in Tianjin
-
摘要: 依据塘沽台所在区域的水文环境、地质构造和动力学等特征，讨论了地表储水能力、地下静水位及构造变形等因素对视电阻率的影响和塘沽台直流视电阻率下降异常的成因机制。结果表明，长期构造位移变化与地下静水位下降是视电阻率呈下降趋势的主要原因；短期降水量大且雨季长，加之有利的地表储水能力是造成塘沽台视电阻率破年变的主要原因。虽然对于塘沽视电阻率异常中是否掺有前兆信息尚缺少确凿证据，但本文采取的分析方法仍可为典型的异常形态提供有益的参考。Abstract: The decline of DC apparent resistivity is a common abnormal variation in electromagnetic forecast, the reason is the change of underground electrical structure caused by rock stress change. Therefore, it is very important for timely identification and scientific judgment of earthquake precursory anomalies to effectively analyze their genetic properties. The EW direction of the resistivity of Tanggu station showed a downward trend since 2013, and the low value in 2021 was significantly lower than previous years, which was a breaking anomaly in the downward trend. The low range was 0.104 Ω·m, accounting for about 1% of the background value, and the low value time was delayed compared with previous years. Based on the hydrologic environment, geological structure and dynamic characteristics of the area where Tanggu station is located, we discussed the influence and genetic mechanism of surface water storage capacity, underground static water level and tectonic deformation. The results show that the long-term tectonic displacement change and the decline of groundwater static water level are the main reasons for the decline of the resistivity trend, and the large short-term precipitation and long rainy season, together with the favorable surface water storage capacity, are the main reasons for the formation of the break year change pattern. Although there is no conclusive evidence that is precursory information in the apparent resistivity anomaly of Tanggu station, the analytical method adopted in this paper can still provide useful reference for the typical anomaly pattern.
-
Keywords:
- apparent resistivity /
- Tanggu station /
- decline change /
- electrical structure /
- lag effect
-
引言
气体构造地球化学方法已广泛应用于构造活动监测、隐伏型活动断裂探测等地震地质相关研究之中（Fu et al，2005，2017；Weinlich et al，2006；Lombardi，Voltattorni，2010；Li et al，2013；Neri et al，2016）。地下气体组分变化能够客观、灵敏地反映地壳的应力、应变状态（李营等，2009），Rn作为土壤气体的重要组分，被认为是构造地球化学观测的有效指标之一。222Rn是一种无色、无味的放射性稀有气体，主要来自于土壤或岩层中238U衰变系列中226Ra的自然衰变（Nazaroff，1992；Jaishi et al，2014）。自然界中222 Rn自然丰度较高且半衰期大（3.823天），而在封闭环境中含量最高（Sciarra et al，2018）。研究表明，Rn是断裂和地震活动最有效的示踪剂之一，对地震构造过程有重要的指示作用（Ciotoli et al，2007），Rn浓度的重大变化与地壳构造活动之间存在某种内在联系。通常在地震活动前后，在相关活动断裂带上能捕捉到Rn浓度异常信息（Barkat et al，2018；Tareen et al，2019；Papachristodoulou et al，2020）。同时，Rn的时空分布特征能够很好地指示隐伏型断裂的构造位置及反映地震活动前的构造应力积累状态（Burton et al，2004；Fu et al，2008；Walia et al，2008；Yuce et al，2017；Li et al，2019）。但是，土壤Rn含量不仅受到断裂带岩层的渗透性、构造活动性等地质条件的制约，还受到沉积物厚度、土壤颗粒组成、孔隙度、放射性物质含量以及温度、气压、降水量等气象要素的影响（Perrier et al，2009；Kumar et al，2017；Sciarra et al，2017；Kandari et al，2018；Sun et al，2018；Yang et al，2018；Chowdhury et al，2019）。通常不同断裂、不同区域，以及同一区域不同断裂构造之间土壤Rn含量及空间分布特征迥异。汤阴地堑处于地层隆起区与沉降区的过渡地带，地质环境特殊，构造复杂。目前，在该地堑边界断裂—汤东、汤西以及安阳南断裂开展的土壤Rn研究结果表明，土壤Rn含量能够敏感指示断裂位置，且三条断裂带之间活动背景差异明显（胡宁等，2019；王明亮等，2019）。为进一步了解不同断裂带间活动背景的差异，需对汤阴地堑土壤Rn的空间分布特征进行深入研究。
鉴于此，本文拟通过在汤阴地堑南部网格化布点并野外监测土壤Rn浓度，分析土壤Rn浓度的空间分布特征，阐明不同断裂带和不同区域气体活动背景及活动程度，揭示其与断裂构造、岩性、覆盖层厚度等之间的联系，和以期为本区域的气体地球化学监测及防震减灾相关工作提供科学依据。
1. 研究区概况和研究方法
汤阴地堑位于太行山隆起区与华北平原沉降区的交接部位（图1），是太行山隆起和内黄隆起间的一个北东向凹陷，其东西边界分别受汤东和汤西断裂控制，南北分别以新商断裂和安阳断裂为界。受NNE−NE向断裂控制，新构造时期形成北东向地堑，基底为三叠系地层，主要发育于古近纪，最大沉积厚度达2500 m；新近纪以后，地层在整体下沉的基础上，具继承性下降的趋势，最大沉积厚度约800 m，其中第四系沉积厚度较薄（中国地震局地球物理勘探中心，2016）。汤阴地堑及其两侧的地层岩性和沉积层厚度差别巨大，西部太行山区，地表主要出露地层为太古界、元古界震旦系、古生界寒武系、奥陶系、石炭系和二叠系，局部为中生界地层；东部丘陵及平原区，主要出露地层为新近系和第四系（刘保金等，2012）。
![]()
图 1 汤阴地堑地质构造及监测点分布图F1：汤西断裂；F2：汤中断裂；F3：汤东断裂；F4：新商断裂；F5：盘古寺断裂；F6：凤凰岭断裂；F7：朱营断裂；F8：薄壁断裂；F9：九里山断裂；F10：百泉断裂，下同。底图引自中国地震局地球物理勘探中心（2016）Figure 1. Geological map and sampled points plot of Tangyin grabenF1：Tangxi fault；F2：Tangzhong fault；F3：Tangdong fault；F4：Xinshang fault；F5：Pangusi fault；F6：Fenghuangling fault； F7：Zhuying fault；F8：Bobi fault；F9：Jiulishan fault；F10：Baiquan fault，the same below. Modified after Geophysical Exploration Center，China Earthquake Administration （2016）新乡市活断层探测资料表明（中国地震局地球物理勘探中心，2016），研究区内主要断裂带为：汤西断裂带、汤中断裂带、汤东断裂带以及新商断裂带，除汤西断裂带有局部露头之外，其余均为隐伏型断裂。其中：汤西断裂带是一条走向NNE、倾向SE的高角度正断层在研究区内卫辉市西代村十里河桥下有出露，该处槽探结果表明其最新活动时代为中更新世，在隐伏区的浅层人工地震显示其可分辨的上断点埋深为65—200 m；汤中断裂带为汤阴地堑内部一条走向NNE、倾向SE的正断层。浅层人工地震探测结果表明该断裂可分辨的上断点埋深为135—180 m，其最新活动时代为早、中更新世；汤东断裂带为汤阴地堑东边界，是一条走向NNE、倾向NW的铲型正断裂。跨断层联合钻孔结果表明上断点埋深为77—103 m，最新活动时代为晚更新世早期。此外，深地震反射剖面表明，汤东断裂带下方存在一条高角度断裂带，其切割了岩石圈底界，属于岩石圈尺度的深大断裂，存在大地震发生的构造背景。344年卫辉东M6、1737年封丘M5
${\raise0.5ex\hbox{$\scriptstyle 1$}\kern-0.1em/\kern-0.15em\lower0.25ex\hbox{$\scriptstyle 2$}}$ 、1814年汤阴M5${\raise0.5ex\hbox{$\scriptstyle 1$}\kern-0.1em/\kern-0.15em\lower0.25ex\hbox{$\scriptstyle 2$}} $ 及2008年封丘ML4.8等地震活动均与汤阴地堑边界断裂活动密切相关。在汤阴地堑南部布设土壤气观测网格，测点间距约1.5 km，共计布设观测点380个（图1）。采用野外流动观测的方法于2019年3月底至4月底进行土壤Rn浓度测量。野外观测时段研究区气象条件相对稳定，新乡市辖区平均气温为13.1—18.8℃ ，降水量为11.7—60 mm，观测期间避开降雨时段，以减小降水对观测结果的影响。测量前按照仪器维护要求进行校准，全部测点使用同一台测氡仪（Alpha GUARD P2000）观测，仪器灵敏度为50 cpm/（kBq·m−3），仪器校准误差小于3%。测量时先用钢钎在测点并列打两个深度约为80 cm的钻孔，将气体取样器置于孔内，封住孔口，用硅胶管连接测氡仪，现场测定Rn浓度。由于土层气体释放不稳定，监测过程中观测数据不断变化，因此在每个测点（15分钟内）进行10次Rn数据观测。在数据分析处理时取后8个观测数据的平均值作为观测值并用ArcGIS10.2，EXCEL和SPSS13.0软件对数据进行统计分析处理，不同断裂带间Rn浓度的比较采用方差分析。
2. 结果与分析
2.1 土壤Rn浓度及其异常阈值
汤阴地堑南部380个土壤Rn观测数据统计分析结果如表1所示。可以看出Rn浓度介于3.09—78.54 kBq/m3，且Rn浓度主要分布于19.10—35.56 kBq/m3。Q-Q （分位数）图分析结果表明（图2），全部观测数据可分为3部分，其中，B部分为气体活动背景，C部分为高于背景的异常值，Rn浓度异常阈值下限为48.40 kBq/m3，背景均值（B部分平均值）为27.22 kBq/m3。A，B和C数据占比分别为5.5%，87.7% 和6.8%。其中，Rn浓度小于10 kBq/m3的低值异常测点，一半以上分布在断裂带及其附近，受断裂破碎带的影响，Rn浓度被大气稀释，浓度降低；而另外一部分低值异常数据主要分布在研究区东部的非断裂带区域，可能受到表层土壤湿度的影响（此部分测点地表大气湿度显著高于其它测点），表层土壤孔隙部分被水分封闭，气体逸出性降低，使原本属于低值背景的测点偏离低值背景。
表 1 汤阴地堑南部土壤Rn浓度分布特征Table 1. The soil radon concentration statistical characteristics of southern Tangyin graben测点数 最大值
/（kBq·m−3）最小值
/（kBq·m−3）平均值
/（kBq·m−3）中值
/（kBq·m−3）下四分位
/（kBq·m−3）上四分位
/（kBq·m−3）标准差
/（kBq·m−3）背景值
/（kBq·m−3）异常阈值
/（kBq·m−3）异常
衬度全部测点 380 78.54 3.09 28.27 27.15 19.10 35.56 12.96 27.22 48.40 2.13 西部测点 111 67.61 3.09 34.11 33.76 24.16 43.55 14.41 33.36 58.51 1.92 东部测点 269 78.54 4.88 25.86 25.18 18.20 32.69 11.50 25.40 44.59 2.11 汤西断裂 36 58.58 7.19 32.85 31.88 25.75 41.85 12.02 32.96 54.80 1.72 汤中断裂 43 51.91 7.88 26.43 24.93 19.78 30.41 11.13 26.59 46.87 1.89 汤东断裂 68 71.40 7.16 23.28 22.58 14.34 29.19 11.80 22.38 38.64 2.38 ![]()
图 2 汤阴地堑南部的土壤Rn分布Q-Q图Figure 2. The Q-Q plots of soil radon concentration in southern Tangyin graben2.2 土壤Rn空间分布及主要断裂的异常衬度
土壤Rn浓度空间分布如图3所示，可以看出，Rn浓度大于异常阈值下限的测点主要分布于西部覆盖层较薄的区域，其次主要分布于汤中断裂与汤东断裂及其邻近区域。以第四系等厚线为依据将研究区划分为西部（等厚线≤50 m）和东部（等厚线＞50 m）两个区域，以分析地质条件对观测结果的影响，其分区统计结果列于表1。由于分区测点Rn浓度与沿断裂带测点Rn浓度均符合正态分布，因此采用平均值加2倍均方差的方法确定Rn浓度异常值，将去除异常值与低背景值（图2中A部分）后的平均值作为背景值。西部测点Rn浓度介于3.09—67.61 kBq/m3，背景均值为33.36 kBq/m3。从图3还可以看出，有近1/3的大于异常阈值的测点分布于汤西断裂带邻近区域，其余测点则呈斑块状分布于断裂带外围，汤西断裂对Rn浓度分布的约束能力在一定程度上被较高的背景所掩盖。东部测点Rn浓度介于4.88—78.54 kBq/m3，背景均值为25.40 kBq/m3。Rn浓度大于异常阈值的测点主要沿汤中断裂、汤东断裂以及新商断裂分布，构造对气体高值异常点的分布具有重要的控制作用。此外，在东部地区还有部分大于异常阈值的测点分布于断裂带外围，这些测点和汤中、汤东断裂附近的部分测点组成与第四系等厚线近似平行的异常条带，推测在此区域存在NE向隐伏断裂。综合分析两个区域的土壤Rn分布特征，其背景均值、异常阈值下限等均呈自西向东递减的变化趋势，且东部地区的高值异常主要沿断裂带分布。结果表明，研究区内自西向东随着覆盖层厚度的增加，由覆盖层厚度较薄导致的岩层出露对土壤Rn浓度的影响逐步减弱，而断裂构造对气体含量高值的控制作用逐步增强。
![]()
图 3 汤阴地堑土壤Rn浓度空间分布Figure 3. The spatial distribution plot of soil radon concentration in Tangyin graben沿断裂带走向，以断裂带为中心，选取东西两侧各1.5 km范围内的测点的Rn浓度，计算气体异常衬度（大于异常阈值测点Rn浓度的平均值与背景均值的比值），来进一步考察各主要断裂带Rn的异常程度。统计分析结果列于表1。研究区内三条主要断裂的异常衬度表现为汤东（2.38）大于汤西（1.72）、汤中（1.89），表明汤东断裂带Rn浓度异常程度最高，这与汤东断裂是晚更新世以来的活动断裂，且在其下方存在一条高角度、切割了岩石圈底界的深大断裂的构造背景相吻合。
2.3 土壤Rn浓度地统计分析
运用ArcGIS10.2地统计分析模块，采用径向基函数规则样条空间插值方法进行地统计分析，生成Rn浓度空间等值线图（图4）。由图4可见，研究区西部一部分等值线高值中心走向与汤西断裂带走向一致，另一部分高值中心受到出露地层的影响在断裂带外围呈斑块状分布。在研究区中东部地区，等值线高值中心连线分布规律表现为两类：① 高值中心连线与断裂带（汤中、汤东等）走向基本一致，反映出断裂构造对土壤气释放的控制作用；② 在新乡市与卫辉市之间，高值中心连线与第四系等厚线近似平行，推测此区域应该存在一条规模较大的隐伏断裂。至于该隐伏断裂存在的合理性在其它方面也有据可循：① 从地质构造图（图1）可以看出，异常条带所在区域第四系等厚线较密集、梯度较大，表明第四系地层厚度变化较剧烈，可能存在断裂构造活动；② 在研究区及其邻近区域，历史上共发生M6地震两次，分别为公元344年12月的卫辉东和1587年4月的修武东地震，两次地震分别位于Rn异常条带两端，距离研究区边缘分别为5和15 km，两次地震震中连线与异常条带近乎完全重合，地震活动或许与推测隐伏断裂有关；③ 从构造应力场方面进行分析，研究区存在出现NE-SW向次级断裂的力学基础。研究区及其附近中小地震震源机制力轴张量反演结果（中国地震局地球物理勘探中心，2016）表明，该区域的最大与最小主应力方位角分别为77°和163°，且最大和最小应力轴接近水平构造应力场。在此构造应力作用下，NNE 向的汤西断裂、汤东断裂发生右旋走滑变形，在地堑内部必然产生近NS 向的拉分作用，从而在新乡与卫辉之间易发生ENE 向张性走滑构造变形。
![]()
图 4 汤阴地堑南部土壤Rn空间等值线及结果解释Figure 4. The contours map of soil Rn concentration and geophysical interpretation for southern Tangyin graben3. 讨论与结论
3.1 Rn浓度空间差异
土壤Rn浓度不仅与构造密切相关而且受到沉积层厚度及岩石、土壤放射性物质含量的影响。本文的监测结果表明，在汤阴地堑西部、东部土壤Rn浓度存在显著差异，西部Rn浓度显著高于东部（P＜0.01）。这种差异与研究区地质背景密切相关，研究区西部处于太行山隆起与汤阴地堑交接的过渡部位，第四系覆盖层厚度较薄，寒武系和奥陶系灰岩大面积出露，与西部毗邻的太行山山前的前寒武系地层中片岩、片麻岩、麻粒岩、混合花岗岩等不同程度的变质岩广泛分布，局部地区零星出露新生代岩浆岩。Choubey等（1999）在加瓦尔喜马拉雅的研究表明，土壤Rn浓度受岩石、构造、铀含量以及土壤厚度等因素的控制，其中Rn浓度和铀含量之间呈显著正相关关系。Kumar等（2017）的研究也表明，Rn浓度与活动断裂的位置密切相关，且Rn的析出率与土壤镭含量以及土壤孔隙度显著正相关。土壤Rn主要来源于铀、钍等放射性物质衰变。岩石类型及分布是评估Rn浓度的重要影响因素（Choubey et al，1999）。通常，铀、镭、钾等放射性物质在三大岩类中的含量依次为岩浆岩＞变质岩＞沉积岩（Barnet et al，2018；Giustini et al，2019），土壤Rn含量在花岗岩、变质岩地区最高，而在沉积岩地区较低，Rn含量与铀含量显著正相关（Ramola et al，2006）。研究区西部及其毗邻地区大面积出露的沉积岩及变质岩中放射性物质衰变产生的Rn通过岩石裂隙或覆盖的土壤层孔隙以扩散和平流迁移的方式直接向地表释放，从而造成在西部浅覆盖地区土壤Rn高值不仅沿断裂带分布，还受覆盖层不连续的影响，在断裂带外围呈斑块状分布，从而造成了土壤Rn浓度背景值的整体抬升。同时，西部土壤Rn浓度的分布特征还受到人类活动的影响。奥陶系灰岩广泛出露的自然背景，为石料开采提供了便利条件。在研究区西部存在不同规模的采石场（图4），在灰岩的开采过程中，爆破作业使周围基岩产生大量裂隙，增加了薄覆盖层下伏岩层的渗透性，为Rn的运移提供了通道，从而使采石场周围的土壤Rn浓度呈现斑块状高值。
东部地区是汤阴地堑主体发育区，在新构造运动的作用下，在三叠系基底的基础上沉积了最大厚度达2500 m的古近系地层；而新近纪以后，在整体下沉的基础上，保持继承性下降的趋势，最大沉积厚度达800 m，其中第四系厚度普遍在100 m以上。巨厚的沉积层抑制了Rn的逸散，而断裂构造的存在大大增加了岩层的渗透性，为Rn释放提供了有利条件，从而造成在中东部地区土壤Rn浓度一部分高值点基本沿着汤中断裂、汤东断裂以及新商断裂等构造带分布，另一部分高值点与第四系等厚线近似平行的格局。Choubey等（2007）在喜马拉雅构造活动区分析岩石和断裂构造对土壤Rn浓度的影响时也得出，在变质程度较低的岩石类型中，在无断裂构造影响下，土壤Rn浓度相对较低；而在石英岩（控制铀的矿化）伴随花岗岩侵入地区，在深部断裂的作用下土壤Rn的浓度达到最高。
3.2 断裂构造与土壤Rn浓度
断裂带的存在增强了岩层及土壤的渗透性，为地下不同来源气体的迁移提供了良好的通道。岩石破裂实验显示，在应力作用下，随着宏观裂隙的发展，为岩石Rn析出不断提供新的表面，使观测Rn浓度显著增加（Tuccimei et al，2010；Mollo et al，2011）。断层气观测实践表明，在断裂带附近可观测到土壤Rn等气体组分浓度异常（Fu et al，2005；Kumar et al，2017；Yuce et al，2017），断裂构造对土壤Rn释放具有重要的控制作用。本文观测结果也表明，在断裂带附近均出现了土壤Rn高值异常，等值线空间分布也显示，断裂带走向与高值中心点连线走向基本一致，显示了断裂带对土壤气释放的支配作用。“新乡市活断层探测与地震危险性评价”项目在本次土壤气监测区内布设的浅层人工地震测线（图4，5）（中国地震局地球物理勘探中心，2016）。对比分析表明，在研究区中东部地区气体高值中心与人工地震显示的断点之间有较好的对应关系，但个别断点或断裂带附近Rn浓度值较低，这可能是断裂带局部地段受到大气的稀释作用，或者是由右旋走滑的构造变形使破裂面局部被断层泥等细粒物质胶结闭合，气体不易逸散所致。从人工地震勘探与气体监测结果的角度分析，中东部地区Rn浓度高值与主断裂、次级断裂或者推测隐伏断裂之间联系密切。
![]()
图 5 浅层人工地震剖面Ⅱ（a）和Ⅲ（b）解释断点Figure 5. Interpretation of faults located on shallow artificial seismic profiles Ⅱ （a） and Ⅲ （b）从研究区内各主要断裂带土壤Rn释放强度来看，汤西断裂Rn浓度显著高于汤东断裂（P＜0.05），Rn浓度平均值、背景均值均表现为汤西断裂＞汤中断裂＞汤东断裂。但仅仅从断裂带气体浓度值来评估气体构造活动性往往存在一定的局限性，Rn浓度值不仅受到断裂构造的控制，还受到Rn释放源、岩层的物理属性、气象条件以及人类活动等因素的综合影响，采用Rn浓度相对值（衬度）来评价断裂带的气体构造活动程度更为合理，评价结果可用于不同断裂之间活动性的对比。Rn浓度衬度结果表明汤东断裂带的气体异常衬度最高，汤西断裂带较低，且整个汤阴地堑南部气体异常衬度东部高于西部。此现象是汤阴地堑南部构造活动背景的反映。“新乡市活断层探测与地震危险性评价”（中国地震局地球物理勘探中心，2016）的探测结果表明，汤东断裂带上浅层人工地震显示的主断裂断点FP27，FX9，FP16 （图4，5）均错断中更新世底界，相应的上断点埋深分别约为90，110和90 m。而断点FP13错断晚更新世底界地层，相应的上断点埋深约为50 m。该断裂带上的两条跨断层联合钻孔（图4）结果表明，南部张河村剖面主断裂错断层位为中更新统第4段底界，上断点埋深为103 m。北部邢李庄剖面（图4）主断裂错断层位为上更新统第1段底界，上断点埋深为77 m。表明汤东断裂带最新活动时代为晚更新世早期。同时，深地震反射剖面（图4中Ⅵ剖面）反映的汤阴地堑南部深部孕震环境显示，汤东断裂下方存在一条高角度、切割了岩石圈底界的岩石圈尺度的深大断裂，汤西断裂与汤中断裂分别在8—10 km和5—6 km深度处收敛归并到汤东断裂上。而汤西断裂带上，西代村十里河桥下的断层出露处的探槽（图4）探测结果揭示，该处中更新统地层与汤西断裂呈断层接触，全新统覆盖于断裂之上，该断裂在中更新世有过活动。断裂南部（研究区之外）的浅层人工地震结果显示，汤西主断裂错断第四系底界，可分辨的上断点埋深约为75 m。汤中断裂带浅层人工地震表明，主断裂断点FP23，FE9 （图4，5）分别错断早更新世底界与第四系底界，上断点埋深分别约为145和180 m，表明汤中断裂带活动时代为早、中更新世。因此，从浅部构造展布和深部孕震环境来看，汤东断裂带规模大、埋藏深度小、切割深度深，其活动性最强，344年卫辉东M6地震震中位于该断裂带东部10 km处。汤东断裂下伏深大断裂的存在为深部来源气体的运移提供了良好通道，深部来源的Rn在载气（CO2等）的作用下，长距离移动至地表释放，显示出较高的气体异常衬度。因此，土壤Rn异常衬度的空间分布特征是研究区地质构造背景的地球化学表现。
综上，汤阴地堑南部土壤Rn浓度，在研究区西部（以第四系等厚线50 m为界），受到岩石单元和人类石料开采活动的影响，Rn浓度背景值高于东部。在西部高浓度背景影响下，Rn浓度高值异常点除沿汤西断裂带分布外，还在断裂带外围呈斑块状分布，断裂带对气体释放的控制作用在一定程度上被掩盖。而东部地区，覆盖层较厚，Rn浓度背景值较低，部分高值异常点主要沿汤中、汤东断裂带分布，显示出构造对气体迁移的控制作用，另一部分高值异常点与第四系等厚线近似平行分布，推测此处存在一条规模较大的隐伏断裂。此外，从研究区主要断裂带Rn异常衬度来看，汤东断裂带高于汤西、汤中断裂带。综合分析研究区地质背景和深部孕震环境认为，Rn异常衬度的这一表现是汤阴地堑南部构造活动背景的反映。土壤Rn能有效用于汤东活动断裂带的构造活动监测，而对位于隆起区与沉降区过渡地带、局部出露于地表，且受人类活动影响较大的汤西断裂带要充分考虑环境背景的影响。
-
![]()
图 9 GNSS基线值在日本MW9.0地震前后的变化
Figure 9. Changes of GNSS baseline values before and after the MW9.0 earthquake in Japan
![]()
图 1 区域构造背景与电性结构
（a） 区域主要断裂、台站位置及2000—2022年M＞1.0地震分布图；（b） 观测区视电阻率布极图；（c） NS向与EW向电测深曲线图
Figure 1. Regional tectonic background and electrical structure
（a） Distribution of major faults，station locations and M＞1.0 earthquakes during 2000−2022； （b） The electrode layout of apparent resistivity at the observation area ；（c） Electrical sounding curves in NS and EW directions
![]()
图 2 塘沽台视电阻率多年日均值曲线图
A，B为往年数据低值对比辅助线，C，D，E，F，G为历年低值时间对比辅助线
Figure 2. Curves of multi-year daily mean value for apparent resistivity at Tanggu station
A and B are contrast auxiliary lines of low-values of previous years，and C，D，E，F， and G are contrast auxiliary lines of low-value time of previous years，respectively
![]()
图 3 塘沽视电阻率夜静时段（1:00—4:00）平均值数据
虚线为理想的虚拟年变趋势线
Figure 3. Apparent resistivity data of Tanggu station at night （1:00—4:00 am）. Dashed line is an ideal virtual annual trend line
![]()
图 4 塘沽台NS向和EW向测线附近的水体位置分布图
长方形与不规则形状表示池塘或水渠，点划线为河流，直线为布极方向线，虚线为金属管线
Figure 4. Water locations near the NS and EW survey lines at Tanggu station
Rectangles and irregular shapes represent ponds or canals，dot dash lines are rivers，straight lines are survey lines，dotted lines are metal pipelines
![]()
图 5 2017—2021年塘沽台降雨量统计图
黑色矩形为最大降雨月份后降雨量持续影响视电阻率变化的部分；红色竖虚线代表每年最大降雨量月份对应的视电阻率位置；红色横虚线为往年视电阻率低值辅助线
Figure 5. Statistical map of rainfall at Tanggu station from 2017 to 2021
The black box is the part where the rainfall continues to affect the change of apparent resistivity after the maximum rainfall month； the vertical red dotted line represents the position of apparent resistivity corresponding to the month of maximum rainfall of each year；the horizontal red dotted line is the low value auxiliary line of apparent resistivity in previous years
![]()
图 6 2017—2021年7—10月塘沽台降雨量统计图
线A表示2021年降雨量，线B表示2017—2020年降雨量均值
Figure 6. Statistical diagram of rainfall at Tanggu station from July to October during 2017 to 2021
Line A refers to the rainfall in 2021，line B refers to the mean rainfall during 2017—2020
![]()
图 7 塘沽台EW向视电阻率与静水位对比图
Figure 7. Comparison of EW direction apparent resistivity and static water level at Tanggu station
![]()
图 8 2002—2022年塘沽台视电阻率长趋势曲线形态
Figure 8. Morphology of long trend curves of apparent resistivity at Tanggu station during 2002−2022
![]()
图 10 塘沽台2018年以来视电阻率EW向数据图像
黑色曲线为视电阻率EW向数据的多项式拟合曲线；虚线A为多年低点趋势性下降的低值辅助线；虚线B为2021年非拟合低值的实际位置，该点略低于A；虚线C为2019年非拟合低值的实际位置，该点略高于A
Figure 10. EW trend data images of apparent resistivity at Tanggu station since 2018
The black line is the polynomial fitting curve of the EW direction data of the apparent resistivity；the dotted line A is the moving average of multi-year low values trending downward；the dotted line B is the actual location of the unfitted low value in 2021，which is slightly below A；the dotted line C is the actual location of the unfitted low value in 2019，which is slightly above A
![]()
图 11 无人机航拍照片与土方堆积的实地拍摄照片
航拍图像叠加在百度地图上
Figure 11. Aerial photos of UAV and field photos of earthwork accumulation
Aerial images are overlaid on Baidu Maps
![]()
图 12 1976年唐山MS7.8地震前后塘沽台视电阻率变化
Figure 12. Apparent resistivity variation of Tanggu station before and after the Tangshan MS7.8 earthquake in 1976
![]()
图 13 塘沽台视电阻率EW向与地震时间对应关系
地震b，c相距34天可视为一组，地震d，e相距21天可视为一组，图中字母对应表2中的地震事件
Figure 13. Correspondence between EW direction of apparent resistivity at Tanggu station and historical earthquake time
Earthquakes b and c can be considered as a group of 34 days apart，so as the earthquakes d and e of 21 days apart，letters mean the events listed in Table 2
表 1 塘沽台降雨量数据统计表
Table 1 Statistical table of rainfall data at Tanggu station
月份 降雨量/mm 2017年 2018年 2019年 2020年 2017—2020年均值 2021年 1月 0.5 0 0 0 0.1 0 2月 0 0 0 0 0 0 3月 4.4 0 0.7 8.1 3.3 10 4月 2.0 31.1 12.7 15.9 15.4 18.6 5月 15.7 21.9 14.8 57.3 27.4 8.5 6月 56.8 56.2 18.1 5.9 34.3 30.7 7月 60.7 130.1 155.7 98.9 111.4 298.9 8月 236.4 139.3 58.3 176.0 152.5 88.7 9月 0 13.9 18.9 147.6 45.1 128.3 10月 62.5 14.1 4.1 7.0 21.9 86.5 11月 0 3.3 1.3 30 8.7 7.2 12月 0 0 1.9 0 0.5 0 总和 439.0 409.9 286.5 546.7 420.5 677.4 
下载: 导出CSV
表 2 2002年以来塘沽台120 km范围地震目录
Table 2 List of earthquakes around Tanggu station with 120 km range since 2002
地震编号 发震日期
年-月-日发震地点 M 震中距/km 走向/° 倾角/° 滑动角/° 震前阶段 震源机制解 A 2006-07-04 河北文安 5.1 112.8 215 87 −115 年变下降段 
B 2020-07-12 河北古冶 5.1 105.6 61 73 −180 年变下降段 
a 2002-05-19 河北丰南 4.2 45.7 318 79 −32 年变下降段 
b 2010-03-06 河北唐山 4.2 99.9 354 45 −33 高值转折段 
c 2010-04-09 河北丰南 4.1 59.7 324 90 −27 高值转折段 
d 2012-05-28 河北唐山 4.8 100.1 314 76 −20 年变下降段 
e 2012-06-18 天津宝坻 4.1 55.6 131 77 −25 年变下降段 
f 2016-09-10 河北唐山 4.2 92.7 152 80 8 低值转折段 
g 2018-02-12 河北永清 4.3 83.9 52 62 −140 年变上升段 
h 2019-12-05 河北丰南 4.5 46.0 325 54 −170 低值转折段 
注：2008年以后地震目录来源于地震编目系统，2002—2007年地震目录来自天津地震局历史地震目录总结
下载: 导出CSV
-
陈峰,马麦宁,安金珍. 2013. 承压介质电阻率变化的方向性与主应力的关系[J]. 地震学报,35(1):84–93. Chen F,Ma M N,An J Z. 2013. Relation between directional characteristics of resistivity changes and principal stress[J]. Acta Seismologica Sinica,35(1):84–93 (in Chinese).
杜学彬. 2010. 在地震预报中的两类视电阻率变化[J]. 中国科学:地球科学,40(10):1321–1330. 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
杜学彬,阮爱国,范世宏,郝臻. 2001. 强震近震中区地电阻率变化速率的各向异性[J]. 地震学报,23(3):289–297. Du X B,Ruan A G,Fan S H,Hao Z. 2001. Anisotropy of the variation rate of apparent resistivity near the epicentral region of strong earthquakes[J]. Acta Seismologica Sinica,23(3):289–297 (in Chinese).
杜学彬,李宁,叶青,马占虎,闫睿. 2007. 强地震附近视电阻率各向异性变化的原因[J]. 地球物理学报,50(6):1802–1810. Du X B,Li N,Ye Q,Ma Z H,Yan R. 2007. A possible reason for the anisotropic changes in apparent resistivity near the focal region of strong earthquake[J]. Chinese Journal of Geophysics,50(6):1802–1810 (in Chinese).
杜学彬,孙君嵩,陈军营. 2017. 地震预测中的地电阻率数据处理方法[J]. 地震学报,39(4):531–548. Du X B,Sun J S,Chen J Y. 2017. Processing methods for the observation data of georesistivity in earthquake prediction[J]. Acta Seismologica Sinica,39(4):531–548 (in Chinese).
高战武,徐杰,宋长青,孙建宝. 2000. 华北沧东断裂的构造特征[J]. 地震地质,22(4):395–404. Gao Z W,Xu J,Song C Q,Sun J B. 2000. Structural characters of the Cangdong fault in North China[J]. Seismology and Geology,22(4):395–404 (in Chinese).
何康,程鑫,李军辉,洪德全. 2010. 安徽省数字化地电阻率干扰与短临异常研究[J]. 地震地磁观测与研究,31(4):86–91. He K,Cheng X,Li J H,Hong D Q. 2010. Research on the short-term and impending anomalies and interference of digital earth resistivity in Anhui Province[J]. Seismological and Geomagnetic Observation and Research,31(4):86–91 (in Chinese).
李金铭. 2005. 地电场与电法勘探[M]. 北京: 地质出版社: 137–140. Li J M. 2005. Geoelectric Field and Electrical Exploration[M]. Beijing: Geology Press: 137–140 （in Chinese）.
邵永新,任峰,陈宇坤,王志胜. 2020. 沧东断裂第四纪活动性探测结果[J]. 地震工程学报,42(4):927–933. Shao Y X,Ren F,Chen Y K,Wang Z S. 2020. Exploration results of Quaternary activities along the Cangdong fault[J]. China Earthquake Engineering Journal,42(4):927–933 (in Chinese).
沈红会,叶碧文,孙春仙,祝涛. 2020. 地电阻率与水位关系的机理分析[J]. 地震,40(4):183–190. Shen H H,Ye B W,Sun C X,Zhu T. 2020. Analysis on the relationship between changes of ground resistivity and water level[J]. Earthquake,40(4):183–190 (in Chinese).
王同利,崔博闻,叶青,李菊珍,王丽红,童琼. 2020. 九寨沟MS7.0地震地电阻率变化时空演化分析[J]. 地球物理学报,63(6):2345–2356. Wang T L,Cui B W,Ye Q,Li J Z,Wang L H,Tong Q. 2020. Temporal-spatial distribution of apparent resistivity before and after the Jiuzhaigou MS7.0 earthquake[J]. Chinese Journal of Geophysics,63(6):2345–2356 (in Chinese).
解滔,卢军. 2016. 地表固定干扰源影响下地电阻率观测随时间变化特征分析[J]. 地震地质,38(4):922–936. Xie T,Lu J. 2016. Apparent resistivity temporal variation characteristics affected by the fixed disturbance source on surface of measuring area[J]. Seismology and Geology,38(4):922–936 (in Chinese).
解滔,卢军. 2020. 含裂隙介质中的视电阻率各向异性变化[J]. 地球物理学报,63(4):1675–1694. doi: 10.6038/cjg2020M0620 Xie T,Lu J. 2020. Apparent resistivity anisotropic variations in cracked medium[J]. Chinese Journal of Geophysics,63(4):1675–1694 (in Chinese).
解滔,于晨,王亚丽,李美,卢军. 2020a. 基于断层虚位错模式讨论2008年汶川MS8.0地震前视电阻率变化[J]. 中国地震,36(3):492–501. Xie T,Yu C,Wang Y L,Li M,Lu J. 2020a. Apparent resistivity variations before 2008 Wenchuan MS8.0 earthquake based on fault virtual dislocation model[J]. Earthquake Research in China,36(3):492–501 (in Chinese).
解滔,王同利,肖武军,胡毅涛,李然,卢军. 2020b. 2020年7月12日唐山MS5.1地震前通州台井下地电阻率变化[J]. 中国地震,36(3):375–382. Xie T,Wang T L,Xiao W J,Hu Y T,Li R,Lu J. 2020b. Apparent resistivity variations before Tangshan MS5.1 earthquake on July 12,2020 recorded at Tongzhou station[J]. Earthquake Research in China,36(3):375–382 (in Chinese).
张国苓,任印国,乔子云,贾立峰,张素欣. 2021. 降雨对阳原地震台地电阻率的影响[J]. 山西地震,(3):22–26. Zhang G L,Ren Y G,Qiao Z Y,Jia L F,Zhang S X. 2021. Influence of rainfall on earth resistivity of Yangyuan[J]. Earthquake Research in Shanxi,(3):22–26 (in Chinese).
张国民, 傅征祥, 桂燮泰. 2001. 地震预报引论[M]. 北京: 科学出版社: 234–238. Zhang G M, Fu Z X, Gui X T. 2001. Introduction to Earthquake Prediction[M]. Beijing: Science Press: 234–238 （in Chinese）.
张明东,张文蕾,曹井泉,王建国,谭青,张玮. 2014. 基于Copula理论挖掘地磁Z分量震磁信息[J]. 地震学报,36(5):930–943. Zhang M D,Zhang W L,Cao J Q,Wang J G,Tan Q,Zhang W. 2014. Mining for seismomagnetic information of geomagnetic Z component based on Copula theory[J]. Acta Seismologica Sinica,36(5):930–943 (in Chinese).
张群伟,朱守彪. 2019. 华北地区主要断裂带上的库仑应力变化及地震活动性分析[J]. 地震地质,41(3):649–669. doi: 10.3969/j.issn.0253-4967.2019.03.008 Zhang Q W,Zhu S B. 2019. The Coulomb stress changes and seismicity on some major faults in North China[J]. Seismology and Geology,41(3):649–669 (in Chinese).
张文朋,张春丽,高武平,闫成国,王志胜. 2022. 用浅层地震勘探资料研究蓟运河断裂的第四纪活动特征[J]. 地震工程学报,44(1):183–191. Zhang W P,Zhang C L,Gao W P,Yan C G,Wang Z S. 2022. Quaternary activity characteristics of the Jiyunhe fault revealed by shallow seismic prospecting data[J]. China Earthquake Engineering Journal,44(1):183–191 (in Chinese).
赵玉林,钱复业. 1978. 唐山7.8级强震前震中周围形变电阻率的下降异常[J]. 地球物理学报,21(3):181–190. Zhao Y L,Qian F Y. 1978. Electrical resistivity anomaly observed in and around the epicentral area prior to the Tangshan earthquake of 1976[J]. Acta Geophysica Sinica,21(3):181–190 (in Chinese).
赵玉林,卢军,张洪魁,钱卫,钱复业. 2001. 电测量在中国地震预报中的应用[J]. 地震地质,23(2):277–285. Zhao Y L,Lu J,Zhang H K,Qian W,Qian F Y. 2001. The application of electrical measurements to earthquake prediction in China[J]. Seismology and Geology,23(2):277–285 (in Chinese).
Brace W F,Orange A S. 1968. Electrical resistivity changes in saturated rocks during fracture and frictional sliding[J]. J Geophy Res,73(4):1433–1445. doi: 10.1029/JB073i004p01433
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
Yamazaki Y. 1965. Electrical conductivity of strained rocks,the first paper,laboratory experiments on sedimentary rock[J]. Bull Earthq Res Inst,43(4):783–802.
计量
- 文章访问数: 239
- HTML全文浏览量: 79
- PDF下载量: 72
