通州台井下地电阻率测深初步分析

2020年7月12日唐山古冶 M_S5．1地震发生前,井下小极距地电阻率出现快速下降-折返的变化,变化形态、异常时长符合地电阻率孕震机理变化,但下降幅度远远小于地表大极距地电阻率孕震过程的变化幅度.因此,文章基于通州台地下电性结构和装置系统,采用数值分析方法,分析地表和井下小极距地电阻率的探测深度.结果表明,当底层电性变化区域介质电阻率发生某种减小时,通州台地表和井下小极距地电阻率装置系统地电阻率观测值会下降,下降的幅度随着深部介质电阻率变化区域上界面向上的扩展而增大.相比地表观测,井下小极距电阻率观测能更显著地接收到深部电阻率变化信号,对孕震有更强的反映能力.     

四平台地电阻率相反年变有限元数值分析

本文依据吉林省四平台地质剖面、地下介质电性资料,将四平台地下介质电性分布简化为沿NW走向的二维电性结构,并结合台站电测深资料建立三维有限元模型,以模型浅层介质电阻率变化模拟测区浅层介质电阻率随季节的变化,计算了四平台N50°W和N40°E两测道地电阻率的年变化形态.计算结果表明,N50°W测道地电阻率在浅层介质电阻率升高时上升,在浅层介质电阻率降低时下降;N40°E测道地电阻率在浅层介质电阻率升高时下降,在浅层介质电阻率降低时升高.两测道地电阻率年变化表现出相反的形态,且N50°W测道年变幅度要大于N40°E测道,计算结果在变化形态上符合四平台实际的观测资料,因此认为四平台测区地下介质的非均匀分布是产生两测道地电阻率年变形态相反现象的原因.     

汶川MS8.0地震前成都台NE测线 地电阻率异常的进一步研究

汶川M_S8.0地震前, 位于震源区内的成都台NE测线地电阻率显示了数年大幅度趋势下降异常, 其年变化甚至严重畸变到消失的程度, 它和多年来我国7级以上大地震记录到的异常形态有极大的相似性, 被认为是最清晰的异常之一. 然而, 长期的观测表明, 成都台NE测线地电阻率与测区潜水水位年变化呈正相关. 因此, 关注异常发展过程中测区水位变化, 成为汶川M_S8.0地震震例地电阻率异常科学总结的重要课题之一. 本文通过对成都台建台以来NE测线地电阻率的全程分析, 确认与汶川M_S8.0地震有关的地电阻率异常变化属于全程变化的第三阶段(1999—2008年). 分析了第三阶段同步观测的地下水位数据, 并进行了地电阻率和地下水位趋势变化的综合分析. 最终确认与汶川M_S8.0地震有关的地电阻率异常经历了前、 后两个时段: 前期2004年2月—2006年2月, 小幅趋势上升, 幅度约为1.9％; 后期2006年2月—2008年2月, 趋势下降, 幅度达-6.8％, 如此大幅度快速下降现象, 是成都台建台30多年来从未见到的. 基于前人开展的试验工作的启示, 推测与此下降相应的测区底层(h＞96.6 m)介质电阻率减小量可达22.3％左右. 其量级之大, 须用承载标本试验中岩石破裂前微裂隙广泛发育和排列趋同化导致的地电阻率变化来解释, 由于成都台位于汶川地震震源区内的特殊位置, 可以作为汶川地震孕震物理过程研究中的重要约束条件之一.      

分层均匀结构地电阻率影响系数一个重要特性普适性的证明

1问题的提出在地电阻率观测中,为分析观测到视电阻率的变化和地下不同区域介质真电阻率变化之间的关系,钱家栋等(1985)提出了影响系数理论。多数情况下,地电阻率观测区可近似为分层均匀的电性结构(图1),对于n层结构,设第i层电阻率为ρ,厚度为h(i=1,2.,…,n),其中底层厚度h_n→∞。     

地电阻率的数值模拟和多极距观测系统

本文利用高精度的多层水平层状介质视电阻率计算公式,模拟了不同结构剖面、不同装置极距、不同规律的真电阻率变化,以及这种变化在介质中所处不同部位条件下地表视电阻率变化的计算结果;讨论了发生于介质中的真电阻率随时间变化在上述各种不同条件下,与在地表观测到的地电阻率(即视电阻率)变化两者之间的关系;指出为了解次异常变化与干扰变化的识别问题,应当将现行单一极距观测改造为多极距观测,并研究和发展相应的地电阻率随时间变化的反演技术和方法,以促进地电阻率法预报地震从看图识字的被动状态向物理预报方向转化。      

地电阻率与地下水位,大气降水关系研究

用固定的地面对称四极装置进行的地电阻率观测,是测量地表以下一定深度范围内介质的电性变化。对于定点测量而言,横向不均匀性不是主要矛盾。因此,电测深资料的数值解释是以层状水平延伸介质为理论模型。在垂向不均匀的电性剖面中,任意一层真电阻率的改变,都会使视电...     

江宁地电台的电性结构及关于地电阻率观测值的解释

为了探究地表和井下观测的地电阻率值存在的明显差异,本文首先对江宁地电台电性结构进行了探测,并建立了水平层状电性结构模型;然后基于该模型计算了地表和井下观测装置的理论视电阻率值;最后对比分析计算值与观测值。结果表明,地表大供电极距观测装置和井下观测装置的观测值均可以由地下电性结构进行合理解释,即地表与井下观测值的差异是合理的,是由地下电性结构所决定的。      

剪切和摩擦滑动大模型的视电阻率变化幅度和各向异性

用石英砂、河砂和水泥模压制成1m×1m×0.3m的均匀和非均匀介质模型，对模型进行剪切和摩擦滑动实验，观测模型不同位置和不同方位的视电阻率变化幅度和各向异性特征. 均匀介质模型实验结果为：电阻率变化与测线的位置有关，距裂隙近的测线，电阻率变化幅度大，从百分之几到百分之几十；距裂隙远的测线，电阻率变化幅度小，从百分之几到不变. 对于非均匀介质模型，电阻率变化幅度要小一些，最大也只有百分之几. 电阻率变化除与测线的位置有关外，还与测线的方位有关，同一测点不同方位的测线，电阻率变化幅度不同，有的差别很大. 剪切与摩擦滑动两种加载方式的电阻率变化幅度数量级相同. 裂隙穿过部位及其附近测点的电阻率变化各向异性主轴方位解与剪切和摩擦滑动的实际裂隙方位吻合较好.     

两次近距离大震前成都台视电阻率重现性、相似性和各向异性变化

本文评价了四川汶川MS8.0、芦山MS7.0地震前后成都台地电观测环境,研究了该台视电阻率变化.结果为:(1)两次大震发生在该台以西的龙门山断裂带、震源机制和震源深度接近,是近距离大震,相应地,在两次地震前该台两个正交测道中的每一测道视电阻率变化均显示了中期异常及其变化过程的重现性和异常幅度的相似性;(2)在每次地震前,两个测道表现了异常变化形态、幅度和起始时间的差异性.其重现性、相似性证明这些异常与两次大震晚期孕育有关;差异性主要展示了与震源机制有直接联系的视电阻率各向异性变化,揭示了震前该台地下介质经历了强烈的电性各向异性变化的过程.     

铁轨干扰浚县地震台地电阻率观测的有限元分析

2016年5月浚县地震台地电阻率EW向测值出现明显变化,为正确评价该变化,依据浚县地震台地质资料、岩层电性资料和EW向电测深曲线,将该台地下介质电性分布简化为二维电性结构,建立三维有限元模型。通过数值模拟,计算测区内地表铁轨干扰对地电阻率观测的影响形态和幅度。结果表明,由数值模拟所得干扰幅度、形态与实际测值的异常变化较吻合,因此从数值分析的角度确定了铁轨对浚县地震台地电阻率观测的干扰。     

松辽盆地岩石圈减薄的深部动力学过程

松辽盆地作为东亚裂谷系的一部分,与华北克拉通一起经历了中生代岩石圈减薄的重大地质事件.对大陆岩石圈-软流圈状态和构造的整体认识,是研究大陆岩石圈减薄深部动力学过程的关键.在获得过松辽盆地的106个宽频和30个长周期大地电磁测深数据的基础上,完成测点数据二维偏离度、构造走向等计算与分析,进一步采用非线性共轭梯度算法,对TE和TM模式数据进行二维联合反演,获得了沿剖面的壳-幔电性结构,并依此构建了松辽盆地壳-幔结构模型.研究结果表明：（1）大兴安岭地区岩石圈厚度约为160 km,松辽盆地岩石圈厚度约为45 km,张广才岭岩石圈厚度在70~100 km之间,莫霍面与岩石圈底界面不呈镜像关系.软流圈整体表现为中、低阻异常,电阻率值在30 Ωm左右,其形态呈西倾约30°的蘑菇状异常,指示了软流圈物质上涌的形式,有别于软流圈垂直上涌的传统认识.（2）松辽盆地深部存在双层高导异常（电阻率小于5 Ωm）,上层为壳内高导层,呈"蛇"状分布,推断为岩浆底侵区,下层为幔内高导层,呈"哑铃"状,为软流圈上涌区.软流圈内存在两个"哑铃"状中、高阻异常,推断为拆沉的岩石圈地幔.具有冷的、高密度的下降物质流的堆积以及拆沉块体下插到两侧山岭是促使大兴安岭与张广才岭在中生代伸展环境中快速隆升重要原因;（3）松辽盆地经历了岩石圈减薄事件,与大兴安岭岩石圈厚度相比,松辽盆地岩石圈厚度减薄了近100 km,与东侧张广才岭相比减薄了70 km,而与中生代华北地台100 km的岩石圈厚度相比,减薄了近50 km,其经历了岩石圈伸展期、裂解期、拆沉期和增长期的动力学过程.     

Stream bottom resistivity tomography to map ground water discharge

This study investigates the effectiveness of direct current electrical resistivity as a tool for assessing ground water/surface water interactions within streams. This research has shown that patterns of ground water discharge can be mapped at the meter scale, which is important for understanding stream water quality and ecosystem function. Underwater electrical resistivity surveys along a 107-m stream section within the Burd Run Watershed in South Central Pennsylvania identified three resistivity layers: a resistive (100 to 400 Ωm) surface layer corresponding to the streambed sediments, a conductive (20 to 100 Ωm) middle layer corresponding to residual clay sediments, and a resistive (100 to 450 Ωm) bottom layer corresponding to the carbonate bedrock. Tile probing to determine the depth to the bedrock and resistivity test box analysis of augered sediment samples confirmed these interpretations of the resistivity data. Ground water seeps occurred where the resistivity data showed that the residual clays were thinnest and bedrock was closest to the streambed. Plotting the difference in resistivity between two surveys, one conducted during low-stage and the other during high-stage stream conditions, showed changes in the conductivity of the pore fluids saturating the sediments. Under high-stream stage conditions, the top layer showed increased resistivity values for sections with surface water infiltration but showed nearly constant resistivity in sections with ground water seeps. This was expressed as difference values less than 50 Ωm in the area of the seeps and greater than 50 Ωm change for the streambed sediments saturated by surface water. Thus, electrical resistivity aided in characterizing ground water discharge zones by detecting variations in subsurface resistivity under high- and low-stream stage conditions as well as mapping subsurface heterogeneities that promote these exchanges.     

Dissolved barium as a tracer of Kuroshio incursion in the Kuroshio region east of Taiwan Island and the adjacent East China Sea

From May to June 2014, the geochemical characteristics of dissolved barium(Ba) in sea water and its influx from the Kuroshio into the East China Sea(ECS) were studied by investigation of the Kuroshio mainstream east of Taiwan Island and the adjacent ECS. This allowed for the scope and extent of the Kuroshio incursion to be quantitatively described for the first time by using Ba as a tracer. The concentration of Ba in the Kuroshio mainstream increased gradually downward from the surface in the range 4.91–19.2 μg L.1. In the surface layer of the ECS, the Ba concentration was highest in coastal water and gradually decreased seaward, while it was higher in coastal and offshore water but lowest in middle shelf for bottom layer. The influx of Ba from Kuroshio into the ECS during May to October was calculated to be 2.19×108 kg by a water exchange model, in which the subsurface layer had the largest portion. The distribution of Ba indicated that Kuroshio upwelled in the sea area northeast of Taiwan Island. The north-flowing water in the Taiwan Strait restrained the incursion of Kuroshio surface water onto the ECS shelf, while Kuroshio subsurface water gradually affected the bottom of the ECS from outside. The results of end member calculation, using Ba as a parameter, showed that the Kuroshio surface water had little impact on the ECS, while the Kuroshio subsurface water formed an intrusion current by climbing northwest along the bottom of the middle shelf from the sea area northeast of Taiwan Island into the Qiantang Estuary, of which the volume of Kuroshio water was nearly 65%. Kuroshio water was the predominant part of the water on the outer shelf bottom and its proportion in areas deeper than the 100 m isobath could reach more than 95%. In the DH9 section(north of Taiwan Island), Kuroshio subsurface water intruded westward along the bottom from the shelf edge and then rose upward(in lower proportion). Kuroshio water accounted for 95% of the ocean volume could reach as far as 122°E. Ba was able to provide detailed tracing of the Kuroshio incursion into the ECS owing to its geochemical characteristics and became an effective tracer for revealing quantitative interactions between the Kuroshio and the ECS.     

华南东部吉安—福州剖面岩石圈电性结构研究

为了研究华南东部地区岩浆活动的深部构造背景,对吉安一福州宽频大地电磁测深剖面数据进行了系统的分析和处理,并利用非线性共轭梯度法进行二维反演,得到了武夷隆起带及周缘地区的岩石圈电性结构;结合区域重磁资料,详细分析了研究区内地壳、上地幔电性结构特征及地质含义.结果表明:华南东部地区岩石圈电性结构存在明显的分区性,并且壳内普遍发育不同成因的高导层,揭示出华南东部地区不同构造单元内的岩浆活动具有不同的成岩构造背景.其中,东南沿海褶皱带深部热侵蚀活跃,岩石圈物质和结构被强烈改造,电阻率普遍较低,软流圈上涌并伴随玄武岩浆底侵,导致岩石圈、地壳剧烈减薄;而武夷隆起带岩石圈电阻率相对较高,印支-燕山早期陆内挤压变形的构造形迹明显,晚中生代岩石圈拉张伸展作用对该地区岩石圈的物质结构有一定的改造.     

Subduction of continental margins and the uplift of high-pressure metamorphic rocks

The mechanism by which high-pressure metamorphosed continental material is emplaced at high structural levels is a major unsolved problem of collisional orogenesis. We suggest that the emplacement results from partial subduction of the continental margin which, because of its high flexural rigidity, produces a rapid change in the trajectory of the descending slab. We assume a two-fold increase in effective elastic thickness of the lithosphere as the continental margin approaches the subduction zone, and calculate the flexural profile of a thin plate for progressive downward migration of the zone of increased rigidity. We assess the effect of changes in the flexural profile on the overlying accretionary prism and mantle wedge as the continent approaches by estimating the extra stresses that are imposed on the wedge due to the bending moment exerted by the continental part of the plate. The wedges overlying the subduction zones, and the subducting slab itself, experience substantial extra compressional stress at depths of around 100 km, and extensional stress at shallower depths, as the continental margin passes through the zone of maximum curvature. The magnitudes of such extra stresses are probably adequate to effect significant deformation of the wedge and/or the descending plate, and are experienced in a time interval of less than 5 m.y. for typical subduction rates. The spatial variation of yield stresses in the region of the wedge and descending slab indicates that much of this deformation may be taken up in the crustal part of the descending slab, which is the weakest region in the deeper parts of the subduction zone. This may result in rapid upward migration of the crust of the partially subducted continental margin, against the flow of subduction. High-pressure metamorphosed terranes emplaced by the mechanism envisaged in this paper would be bounded by thrust faults below and normal faults above. Movement on the faults would have been coeval, and would have resulted in rapid unroofing of the high-pressure terranes, synchronous with arrival of the continental margin at the subduction zone and, therefore, relatively early in the history of a collisional orogen.     

过套管地层电阻率曲线环境影响校正方法研究

过套管电阻率测井是在金属套管井中测量地层电阻率的一种电法测井技术,其测量结果受套管、层厚/围岩、水泥环等井眼环境因素的影响,所以在进行过套管电阻率测井资料解释时必须进行测井环境影响因素校正.本文在论述过套管电阻率测井原理的基础上,针对大庆油田引进的俄罗斯过套管电阻率测井仪器ECOS,分别考察了层厚/围岩、水泥环对测量结果的响应规律,利用正演模拟方法研制完成了层厚/围岩影响校正图版以及水泥环影响校正图版.在计算机上实现了层厚/围岩、水泥环影响校正的算法,并利用模型算例和测井实例验证了本文方法的合理性.该方法能够减少水泥环的影响,提高层厚小于1.0 m薄层电阻率的测量精度,为过套管电阻率测井资料正确解释提供了保障.     

Forward modeling of marine DC resistivity method for a layered anisotropic earth

Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell’s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.     

井下地电阻率观测影响系数分析——以江宁地震台为例

采用水平层状均匀介质中点电源位于任意深度时的电位解析表达式,以江宁台3层电性结构为例,分析了井下对称四极地电阻率观测时各层影响系数随深度、极距的变化,并结合探测深度探讨了实施井下观测时影响系数在选择供电极距和电极埋深时的作用。结果表明,对于"K"型电性结构,江宁台井下观测对地表、浅层干扰有较强的抑制作用,其短极距观测对地表、浅层干扰的抑制能力显著优于长极距观测;长极距观测在电极埋深H小于100m时对地表介质季节性的干扰具有放大作用;浅层影响系数一定时,电极埋深和供电极距需同时增加;江宁台井下观测供电极距AB/2取100～150m、电极埋深H为250m较为合理。     

Geoelectrical structure and hydrogeological investigations of the southern Rif Cordillera (Morocco)

An electrical prospecting survey is conducted in the Rharb basin, a semi‐arid region in the southern part of the Rifean Cordillera (Morocco) to delineate characteristics of the aquifer and the groundwater affected by the marine intrusion related to Atlantic Ocean. Analysis and interpretations of electrical soundings, bi‐logarithmic diagrams and the geoelectrical sections highlight a monolayer aquifer in the southern part, a multilayer system in the northern part of the Rharb basin and lenticular semi‐permeable formations. Several electrical layers have been deduced from the analysis of bi‐logarithmic diagrams: resistant superficial level (R₀), conducting superficial level (C₀), resistant level (R), intermediary resistant level (R′), conducting level (Cp) and intermediary layer of resistivity (AT). Spatial distribution of the resistivity deduced from the interpretation of apparent resistivity maps (AB = 400 and 1000 m) and the decreasing of resistivity values (35–10 Ωm), in particular in the coastal zone show that this heterogeneity is related to several anomalies identified in the coastal area, which result from hydraulic and geological processes: (i) heterogeneous hydraulic conductivity in particular in the southern part of the Rharb; (ii) lateral facies and synsedimentary faulting and (iii) the relationship between the electrical conductivity and chloride concentration of groundwater shows that salinity is the most important factor controlling resistivity. The distribution of fresh/salt‐water zones and their variations in space along geoelectrical sections are established through converting subsurface depth‐resistivity models. Copyright © 2010 John Wiley & Sons, Ltd.