Monitoring of relative changes in the electrical conductivity of rocks from observations of the magnetotelluric apparent resistivity (numerical modeling)

A method for direct conversion of observed variations in the magnetotelluric (MT) apparent resistivity ρ_a into relative variations in the resistivity of elements of a well-studied geoelectric structure is proposed. The method is tested on a 1-D model structure consisting of seven horizontal layers in three of which the conductivity can vary within certain limits. It is inferred that the frequency range and the accuracy of methods of magnetotelluric sounding presently applied to the construction of transfer operators are sufficient for determining relative changes in the resistivity of rocks; the latter, as distinct from ρ_a, can serve as an effective prognostic parameter. The method can be extended to more complex geoelectric structures.     

Modelling the DC electrical response of fully and partially saturated Permo-Triassic sandstone

Results from a laboratory investigation into the electrical properties of fully and partially saturated Wildmoor Triassic Sandstone have been modelled using the Archie, Waxman–Smits and Hanai–Bruggeman equations. The results demonstrate the limitation of using simple relationships to describe samples when the matrix resistivity ρ_r is not significantly greater than the saturating electrolyte resistivity ρ_w. In these situations Archie's parameters m and n are not accurately determined. Conversely, the more sophisticated Waxman–Smits and Hanai–Bruggeman models provide parameters that better describe the electrical properties of the rock and are able to identify heterogeneity between samples that would otherwise be missed. The ranges of values for matrix resistivity (49 < ρ_r < 161 Ωm) and cementation factor (1.6 < m < 2.1) obtained from the Hanai–Bruggeman model indicate significant variation between samples. Comparison of laboratory‐determined values for cation exchange capacity (0.06 < Q_v < 0.51 meq/mL) and those obtained from the Waxman–Smits model (0.09 < Q_v < 0.55 meq/mL) indicates a very strong correlation, suggesting this model is appropriate for describing the rock. There is good agreement between parameters modelled using fully and partially saturated versions of both the Hanai–Bruggeman and Waxman–Smits equations, indicating that the data are consistent with these models and that the assumptions made are appropriate.     

机械排烟作用下无机防烟卷帘防烟特性的实验研究

采用全尺寸模拟实验,通过测量机械排烟作用下无机防烟卷帘背火面壁面的平均温升、卷帘缝隙处的烟气流动等,对无机防烟卷帘在机械排烟作用下的防烟特性展开研究。结果表明：无机防烟卷帘在发生火灾时自动下降至地面,在机械排烟作用下,在20 min的实验时间内,基本上能将烟气阻挡在模拟着火房间内,只有极少量烟气进入模拟中庭,卷帘背火面壁面温度接近环境温度;该无机防烟卷帘具有良好的挡烟耐热性能,火灾时可以保证非着火区域人员的安全疏散。本文研究可为该种新型无机防烟卷帘在复杂中庭防火设计中的应用提供实验支持。     

A NEW PARAMETER FOR THE INTERPRETATION OF INDUCED POLARIZATION FIELD PROSPECTING (time-domain) *

In a previous paper it has been shown that we can relate the transient IP electric field E_p , existing in a rock after a step wave of polarizing current, with the steady-state current density J_ss during the current step wave as follows: E_p =ρ' J_ss . This relation may be interpreted as a generalized Ohm's law, valid in linear cases, in which ρ’(fictitious resistivity) is defined as the product of the true resistivity ρ with the chargeability m. Supposing E _p=— grad U_p and applying the divergence condition div J_ss = o, one can, for a layered earth, obtain a general expression for the depolarization potential U_p as a solution of Laplace's equation ?²U_p= o. Since the mathematical procedure for the solution of this last equation is identical to that used in resistivity problems, we propose now the introduction of an apparent fictitious resistivity ρ'_a (defined as the product of the apparent resistivity ρ_a with the apparent chargeability m_a) as a new parameter for the interpretations of IP soundings carried out over layered structures with a common electrode array. The most general expression of ρ'_a as a function of the electrode distance turns out to be mathematically identical to the general expression of ρ'_a. Therefore it is possible to interpret a ρ'_a field curve using the same standard graphs for resistivity prospecting with the usual method of complete curve matching. In this manner a great deal of work is saved since there is no need to construct proper m_a graphs for the interpretation of IP soundings, as it has been done up to now. Finally some field examples are reported.     

A compact representation of magnetotelluric responses for two-layer models

In this paper the locations where ρ_app = ρ₁ and ? = π/4 and where these parameters reach an extreme value in two-layer magnetotelluric (MT) sounding curves are summarized in an extremely compact form. The key parameters over two-layer models with conductivities σ₁, σ₂ and upper layer thickness h are the real S and α, where S is the conductivity contrast and α is the distance between the observation site and the conductivity interface, normalized to the half skindepth in the first layer. If the impedance components, various resistivity definitions ( ρ_{Re Z}, ρ_{Im Z} and ρ_|Z|, based on different parts of the complex impedance Z ) and the magnetotelluric phase ? are derived as a function of S and α, then the conditions for the apparent resistivity ρ_app and the phase ? are that they either satisfy ρ_app = ρ₁ and ? = π/4 or attain extreme values which can be given in terms of simple algebraic equations between S and α. All equations are valid for observation sites at any depth 0 ≤ z ≤ h in the first layer. The set of equations, presented in a tabular form, may make it possible to determine a layer boundary from the short period part of the sounding curves, in particular the ρ_{Re Z} and the ?_MT curves.     

THE QUANTITATIVE ANALYSIS OF FREQUENCY-DOMAIN INDUCED POLARIZATION SOUNDINGS OVER HORIZONTAL BEDS*

Following up our recent study of an indirect procedure for the practical determination of the maximum frequency-effect, defined as fe = 1 ? pρ_∞/ρdc with ρ_∞ the resistivity at infinite frequency, we show at first how, through the Laplace transform theory, ρ_∞ can be related to stationary field vectors in the simple form of Ohm's law. Then applying the equation of continuity for stationary currents with a suitable set of boundary conditions, we derive the integral expression of the apparent resistivity at infinite frequency ρ_∞,a in the case of a horizontally layered earth. Finally, from the definition of the maximum apparent frequency-effect, analytical expressions of fe_α are obtained for both Schlumberger and dipole arrays placed on the surface of the multi-layered earth section in the most general situation of vertical changes in induced polarization together with dc resistivity variations not at the same interfaces. Direct interpretation procedures are suggested for obtaining the layering parameters directly from the analysis of the sounding curves.     

碳酸盐岩储层孔隙结构对电阻率的影响研究

碳酸盐岩储层孔隙类型多样,各种孔隙的尺寸变化范围可以跨越几个数量级,孔隙结构非常复杂,这种复杂孔隙结构和不均匀分布的多元孔隙空间使得储层电性呈现明显非阿尔奇特性.为了了解影响电阻率变化的控制因素,本次研究选取中三叠世雷口坡组的8块全直径碳酸盐岩岩样,开展了核磁共振、岩电实验、孔渗实验、压汞实验及薄片等实验,并利用数字图像分析法定量分析了孔隙结构特征.研究结果表明：①孔隙度是影响电阻率高低的重要因素,但并非唯一因素,除孔隙度以外,孔隙尺寸和数量、孔隙网络复杂程度远比吼道大小对电阻率的影响大;②在孔隙度一定的条件下,胶结指数m随储层中孤立大孔隙占比的增多而增大,当孔隙度增大到一定程度后,胶结指数m又随大孔隙占比的增多而减小,微裂缝起重要沟通作用;③在给定孔隙度时,以简单大孔隙为主的岩样表现为胶结指数m值较大,而以复杂孔隙网络、细小孔隙为主的岩样表现为胶结指数m值较小,具分散、孤立大孔隙的岩样,胶结指数m值最高;④依据孔隙几何参数与电阻率和胶结指数之间的关系,可以利用测井资料间接判别储层类型,从而提高储层有效性和含水饱和度评价精度.

     

Geophysical methods in drinkwater protection of near-surface reservoirs

Two case histories are presented in which near-surface water-bearing formations and their overburden are examined. This survey utilized a combination of traditional geoelectric methods (direct current resistivity and time domain induced polarization) and a special shallow-depth engineering geophysical sounding method. This latter method was developed in Hungary to investigate near-surface unconsolidated formations such as clay, silt, sand, gravel and other similarly “penetratable” formations.For a gravel terrace, in addition to the usual parameter maps (resistivity, thickness, polarizability), combined multiparametric characterization maps have been plotted. These maps illustrate the hydrogeologic value of the water-bearing formations and the protecting capacity of the overburden better than the conventional parameter maps.Surface measurement results together with filtration coefficient values show how an inexpensive geophysical survey can successfully be applied in estimating the dangerous or protected situation of near-surface reservoirs and how such a survey can be used for siting reservoirs. This work enabled the near-surface geological structure to be determined, hereby providing orientation for more expensive, highly detailed surveys.     

THE POSSIBILITIES OF GEOPHYSICAL METHODS APPLIED FOR INVESTIGATING THE IMPACT OF MAN ON THE GEOLOGICAL MEDIUM*

Geophysical methods can be applied to investigate the harmful effect of man's activities on the environment: the study of specific electrical resistivities and natural electric fields of filtrational origin makes it possible to control the penetration of exogenic pollutants into soil as well as desalinization and secondary salinization of soils; electrometric and seismometric methods allow to observe the groundwater level near water reservoirs thus evaluating the harmful effect of backing up natural groundwater flow; contrasting properties of bedrocks and rocks in a landslide body permit application of seismic and electrical prospecting methods. Observation of the changes in specific electrical resistivities with time on the slopes of quarries is an effective method of assessing slope stability and predicting landslide hazard. Mining activities, groundwater pumping, and oil extraction are the main causes of endogenic pollution of geological medium; surface and borehole geophysical methods make it possible to assess vertical and horizontal displacements of the interfaces between salty and fresh subterranean waters resulting from an extensive exploitation of water supply sources; activation of geodynamic processes associated with mining activities is determined from the data of ‘regime’ electrometric, seismometric, gravimetric, and inclinometric observations. Geophysical methods are applied widely for studying the intensification of physico-geological processes under the impact of man. A specific object of electrometric investigations is degradation of permafrost and intensification of karst processes in soluble rocks. The principal advantage of geophysical investigations lies in the possibility of creating high spatial and temporal density of observations permitting an extensive employment of statistical methods in the assessment of the impact of man on the geological medium.     

低电阻率油气层物理参数变化机理研究

研究了淡水环境下富含粘土矿物、细颗粒砂岩低电阻率油气层. 分析了粘土的附加导电性,认为粘土质量百分数高且富含伊利石和蒙脱石的地层中,若地层水质量浓度小于15g/L等效NaCl盐溶液时,粘土的附加导电性是造成油气层电阻率低的首要原因. 同时分析了细颗粒、粘土质量百分数不高的储层,高的毛细管压力使束缚水饱和度很高,同样能形成低电阻率油气层. 储层中粘土质量百分数的多少对储层岩石电性特征(I-Sw曲线形态)有重要影响,可用于指导测井资料解释模型的选择.     

低电阻率油气层物理参数变化机理研究

研究了淡水环境下富含粘土矿物、细颗粒砂岩低电阻率油气层. 分析了粘土的附加导电性,认为粘土质量百分数高且富含伊利石和蒙脱石的地层中,若地层水质量浓度小于15g/L等效NaCl盐溶液时,粘土的附加导电性是造成油气层电阻率低的首要原因. 同时分析了细颗粒、粘土质量百分数不高的储层,高的毛细管压力使束缚水饱和度很高,同样能形成低电阻率油气层. 储层中粘土质量百分数的多少对储层岩石电性特征(I-Sw曲线形态)有重要影响,可用于指导测井资料解释模型的选择.     

碳酸盐岩储层孔隙结构对电阻率的影响研究

碳酸盐岩储层孔隙类型多样,各种孔隙的尺寸变化范围可以跨越几个数量级,孔隙结构非常复杂,这种复杂孔隙结构和不均匀分布的多元孔隙空间使得储层电性呈现明显非阿尔奇特性.为了了解影响电阻率变化的控制因素,本次研究选取中三叠世雷口坡组的8块全直径碳酸盐岩岩样,开展了核磁共振、岩电实验、孔渗实验、压汞实验及薄片等实验,并利用数字图像分析法定量分析了孔隙结构特征.研究结果表明：①孔隙度是影响电阻率高低的重要因素,但并非唯一因素,除孔隙度以外,孔隙尺寸和数量、孔隙网络复杂程度远比吼道大小对电阻率的影响大;②在孔隙度一定的条件下,胶结指数m随储层中孤立大孔隙占比的增多而增大,当孔隙度增大到一定程度后,胶结指数m又随大孔隙占比的增多而减小,微裂缝起重要沟通作用;③在给定孔隙度时,以简单大孔隙为主的岩样表现为胶结指数m值较大,而以复杂孔隙网络、细小孔隙为主的岩样表现为胶结指数m值较小,具分散、孤立大孔隙的岩样,胶结指数m值最高;④依据孔隙几何参数与电阻率和胶结指数之间的关系,可以利用测井资料间接判别储层类型,从而提高储层有效性和含水饱和度评价精度.     

Lithological and physical properties of core samples from the Sagara oil field: Oil occurrence in Sagara core samples

Abstract   A drilling and coring investigation of the Sagara oil field, central Honshu, Japan, was conducted to contribute to the understanding of hydrocarbon migration processes in a forearc basin. Core samples were analyzed to determine lithology, physical properties (specifically gas permeability) and the characteristics of oil occurrence. Gas permeability values greater than approximately 10⁻¹¹ m² constitute the basic precondition for any lithology to serve as a potential fluid conduit or reservoir in the Sagara oil field. Cores recovered from the 200.6-m-deep borehole were primarily composed of alternating siltstone, sandstone and conglomerate, all of which are correlated to the late Miocene Sagara Group. Both sandstone and conglomerate can be classified into two types, carbonate-cemented and poorly to non-cemented, based on matrix material characteristics. Oil stains are generally absent in the former lithology and more common in the latter. Variations in physical properties with respect to gas permeability values are directly related to the presence and character of carbonate cement, with higher permeabilities common in poorly to non-cemented rocks. The relationships between lithology, oil-staining, cementation and permeability indicate that cementation preceded oil infiltration and that cementation processes exerted significant control on the evolution of the reservoir.     

THEORETICAL AND EXPERIMENTAL INVESTIGATIONS FOR CAVITY RESEARCH WITH GEOELECTRICAL RESISTIVITY METHODS*

The exact localization of subterranean cavities and the determination of their dimensions is very important for the planning of geotechnical and mining activities. It is a complicated geophysical task often at the limit of detection. Nevertheless geophysical investigation is the only alternative to a dense and expensive grid of boreholes. This report tests the usefulness of geoelectrical resistivity methods for cavity detection under some new aspects. The basis for evaluation was a theoretical analysis of different conventional and focussing measuring arrays and of special arrays for a geoelectrical research between two boreholes. The limit of detectability of a cylindrical cavity of defined cross-section and depth was calculated for the different measuring arrays on the basis of computation of the apparent resistivity ρa. Furthermore, the influence of possible errors (current supply of the electrodes and the distance between the electrodes) is discussed for focussed systems. The second part of the article is directed at the behaviour of the apparent resistivity ρa, the disturbing potential δV_d caused by the cavity and the normal potential δV₀ of the measuring array all in relation to a homogeneous earth. Some new results are presented. In the last part of the article theoretical results are compared with some field measurements.     

Impact of grout curtains on karst groundwater behaviour: an example from the Dinaric karst

Grout curtains are vertical grout walls installed in the ground. In karst terrains, their construction is primarily connected with dams and reservoirs. Their main role is to increase water tightness and to prevent progressive erosion, blocking possible seepage paths along karst fissures and conduits. In this article, changes in the behaviour of the groundwater level (GWL) and the water temperature in nine deep piezometers, which were caused by the construction of a grout curtain at the ?ale Reservoir on the Cetina River (Croatia), were analysed. The total length of the grout curtain is 3966 m. It spreads 120 m below the dam. The most analysed data are from the period after the dam had been built. Only few data and figures concern the comparison between pre‐ and post‐dam periods. The hourly data of the GWL and the water temperature were analysed for the period between 1 September 2008 at 02:00 h to 31 December 2009 at 23:00 h (11 687 h total) in six deep piezometers (marked in the text and figures as 1, 2, 3, 4, 5 and 6). For three piezometers (marked in the text and figures as A, B and C), some discontinuous measurements of the GWL and the water temperature were available for analysis. The construction of the grout curtain made strong, sudden and possibly dangerous changes to the characteristics of the aquifer and the circulation of groundwater in the local area. Special attention is paid to analyses of the behaviour of the hourly GWL data measured in the piezometers pairs (two neighbouring piezometers, one inside and the other outside of the grout curtain). During more than 80% of the analysed period, the GWL was higher in the piezometer inside the grout curtain than the one outside of it. The intensity and range of the dynamics of GWL was higher in piezometer outside the grout curtain than the inside ones. After the construction of the grout curtain, the maximum measured hydrostatic pressure on some parts of the grout curtain was approximately 40 m. It changes quickly in both time and direction. The water temperature was found to be similar in all of the measured piezometers, and it varies between 10.2 and 15.7 °C with an average value of 12.7 °C. Copyright © 2011 John Wiley & Sons, Ltd.     

差异性成岩过程对百口泉组砂砾岩岩石物理特征的影响

现有针对储层岩石地震岩石物理性质的研究通常不考虑地质过程的影响,致使研究结果缺乏明确的地质意义,制约了储层地震预测和评价的效果.本文通过对玛湖凹陷百口泉组砂砾岩岩石物理性质的测量,并结合岩石学特征分析及沉积环境的划分,在沉积环境与成岩石作用过程框架下分析岩石物理特征的变化规律.研究结果表明,不同微相砂砾岩初始沉积物矿物组成影响了后期成岩过程,差异性的成岩过程形成不同的岩石微观结构特征.三角洲平原辫状河道与泥石流微相砂砾岩塑性组分含量较高,在机械压实阶段已致密化,成岩作用主要表现为成岩早期的淡水渗滤黏土与火山碎屑蚀变黏土的岩屑颗粒包裹式胶结,具有以塑性黏土作为岩石受力骨架的特征.三角洲前缘砂坝与水下分流河道砂砾岩中刚性花岗岩岩屑、石英以及长石的含量较高,机械压实作用相对较弱,原生粒间孔隙保存较好,同时溶蚀孔隙也较发育;砂砾岩成岩作用表现为从前缘砂坝微相的硅质胶结为主转变为水下分流河道微相的硅质+黏土胶结,具有以刚性石英或石英+黏土作为岩石受力骨架的特征.前缘碎屑流砂砾岩样品则表现为快速堆积的钙质基底式胶结,以钙质胶结物作为岩石受力骨架.由砂砾岩颗粒微观接触特征(胶结特征)所决定的岩石受力骨架控制了宏观岩石物理性质变化规律,表现为具有相同颗粒边界力学特征的微相样品具有一致的速度-密度、纵波速度-横波速度以及孔隙度-速度变化关系,也造成不同微相样品在纵、横波速度比及电阻率值上具有明显的差异,具有分区分布的特征,速度比与电阻率可作为百口泉组砂砾岩储层评价的敏感岩石物理参数.研究结果可为该组砂砾岩储层的测井解释和“甜点”地震预测提供依据,并可作为示例形成在地质过程中研究岩石物理特征演化的新方法.

     

视电阻率法测矿井岩石应变实验的初步总结

本实验目的是揭示受力岩石中应力与电阻率的关系,以验证形变电阻率法能否用于地震预报工作。 实验在一离地面深100多米、正在开采的矿坑中进行。视电阻率的测量用四极对称装置,干电池供电,最高电压13.5伏,最大观测均方误差为1-2％。五个测点上贴了电阻丝应变片,以测定岩石主应变的大小、方向以与视电阻率极大值进行对比。应变计的读数值最大误差不超过3-4με。实验中考虑了岩石有严重极化现象的影响。 实验的初步结果表明:岩石的视电阻率和应变值的变化,总趋势相当一致,所以用形变电阻率法来测量应变是极有希望的。对1973年在四川省南坪发生的几次地震进行震例分析的结果,得到同样的结论。 文中对实验中存在的问题,也作了一些讨论。     

水泥砂浆面层加固砌体墙抗震可靠指标研究

采用水泥砂浆面层加固方法加固砌体房屋是一种简单有效的方法。对比分析了《砌体结构加固设计规范》(GB 50702-2011)^[3]和《建筑抗震加固技术规程》(JGJ 116-2009)^[4]中钢筋网水泥砂浆面层加固砌体墙的可靠性;对10片未加固低强度砖墙和20片单面钢筋网水泥砂浆面层加固低强度砖墙进行拟静力试验,分析了未加固和加固墙体的破坏模式和机理,并提出了适用于钢筋网水泥砂浆面层加固砌体墙的抗震验算公式。研究表明：两标准的可靠水平差别较大,实际应用时,易产生矛盾;未加固墙主要发生沿阶梯形斜裂缝受剪破坏,加固墙主要发生沿通缝受剪破坏和沿阶梯形斜裂缝受剪破坏;建议公式的计算值与试验值吻合良好,未加固墙抗震可靠指标为2.2,加固墙抗震可靠指标为2.5～3.1;砂浆面层加固砌体结构可以显著提高结构的抗震性能。     

水合物赋存形式对石英砂电性特征影响的数值模拟研究

天然气水合物储层精确评价关系到储层的生产与开发过程.由于储层微观结构具有非均质性, 且水合物合成分解过程中分布形式复杂, 水合物储层物性变化规律尚不明晰, 从而难以对储层水合物饱和度进行准确预测.为明确水合物赋存模式对电性特征的影响规律, 本文基于石英砂样品采用随机模拟的方法实现了颗粒包裹型、孔隙填充型和团簇状分布三类赋存模式的水合物数字岩心构建, 并进一步提取孔径分布分析不同水合物分布模式下的微观结构变化特征, 利用有限元方法对含水合物数字岩心进行电性模拟, 结合阿尔奇公式得到适用于含水合物石英砂介质的经验参数取值, 并采用实验室内水合物合成实验电阻率测量数据与实际储层电阻率测井数据对比验证了阿尔奇经验参数取值的正确性及有效性.结果表明, 随着水合物饱和度逐渐增加, 在相同水合物饱和度下, 颗粒包裹型生长引起的微观结构变化对电传导过程的影响最为显著; 孔隙填充型水合物生长引起的孔隙结构变化对电传导过程影响最小; 不同水合物赋存形式下, 沉积物孔径分布变化规律从微观角度解释了三类分布模式的电性差异.研究结果有助于研究储层内含水合物沉积物合成分解过程中的电性特征演化规律.

     

一种评价致密砂岩储层孔隙结构的新方法及其应用(英文)

致密砂岩储层的孔隙结构对其渗透性和电性影响显著,是此类复杂储层岩石物理研究的关键。针对仅从连通喉道半径评价渗透率的多解性以及储层孔隙结构与电性关系研究欠缺等不足,综合影响物性的主要因素,提出了一种同时考虑孔隙度、最大连通喉道半径及分选性三种因素的新型孔隙结构参数δ的计算公式。利用岩心及实测数据对比分析表明,δ值能够较连通喉道半径等传统方法更精确地刻画致密砂岩储层渗透性,同时它与储层电性具有密切关系,可用于估算地层因素F和胶结指数m。据此提出将孔隙结构对电阻率的影响进行归一化校正以及基于核磁共振测井预测储层完全含水电阻率R0的评价方法,从而突出储层流体性质变化引起的电性变化,并提供了一种新的致密砂岩储层流体识别思路,研究结果得到了实验资料和实际测井试油资料的验证。