The effect of a deviatoric stress on physical rock properties

Petrophysical measurements were carried out on dry specimens of mica-gneiss, amphibolite and serpentinite from KTB core samples and samples of surface outcrops in order to determine the effect that a deviatoric stress field, as observed at the KTB area, may have on the in-situ rock properties. Simulating the variation of the actual principal stresses and temperature with depth, seismic wave velocities, densities, linear and volumetric strain (porosity) have been measured, taking into account the overall spatial orientation of the foliation at the KTB area with respect to the principal stress axes. Comparison with respective data evaluated for lithostatic pressure conditions revealed that the stress-related (crack-related) effect on wave velocities respectively on velocity anisotropy is in the range 1–3%, due to microcracks which are selectively closed or kept open by the deviatoric stress. The effect of the deviatoric stress is particularly documented by shear wave splitting due to microfractures that are oriented normal to the minimum principal stress axis.     

Research on dependence of resistivity changing anisotropy on microcracks extending in rock with experiment

Intr0ducti0nTheproblemsaboutcrustalrocksfracturingandfricti0nalslippingalongthefracturesurfaCearecloselyrelatedtotheresearchesofearthquakemechanismandgeol0gichazards,andalsoareoneofthosebeinggivenmoreconcernbyscientistsinearthscience-Discussingab0utthegenera-tionanddevelopmentprocess0fintraPlateearthquakes,Wang,etal(l999)proPOsedthatearth-quake'spregnant,mechanicspropertyandgeometricstructureofmatterinearthquake'sgestatingareaaredifferinthousandswaysf0rdifferentplaceinsideslab-Thus,thepregna…     

Magnetic fabrics and superimposed deformations: example of Dalradian rocks from the southwest Highlands of Scotland

Analysis of anisotropy differences in specimens close in the same sample permits the separation of anisotropies due to ferrimagnetic minerals and their dia-paramagnetic matrix. This method, previously described, has been applied to rocks, affected by superimposed deformations, from the Dalradian of Scotland. It appears that ferrimagnetics anisotropy is often found in relation with the secondary deformation, even though matrix anisotropy is related to the primary deformation. This last anisotropy has here a dominant effect on the orientation of the resultant fabric of the whole rock.     

利用主成分分析法研究余干地震台接收函数

对模型合成的理论接收函数进行主成分分析,研究在倾斜界面和各向异性影响下主成分形态变化规律,分析倾斜界面倾角和各向异性强度对接收函数R分量主成分贡献率的影响。对江西余干地震台的接收函数R分量主成分进行分析,结果发现:Ps震相的平均走时为3.4 s;台站下方介质的各向异性和倾斜界面同时存在,各向异性快轴大致呈NE向,倾斜界面倾向约170°。     

Field-induced anisotropy of magnetic susceptibility in rocks

Summary An induced anisotropy of magnetic susceptibility results from the domain alignment which is produced by treating stationary specimens in a strong alternating field. Appreciable domain re-orientation occurs in fields as low as 50 oersteds and the effect must therefore normally be an important part of the process of alternating field demagnetization. Induced anisotropy has been measured in a number of igneous rocks with a range of palaeomagnetic stabilities and in magnetite powders of controlled grain sizes, dispersed in plaster or kaolin specimens which were mechanically deformed to produce instrinsic magnetic anisotropy by grain alignment. The saturation magnitude of the induced anisotropy is not a function of grain size but the saturating field required increases with decreasing grain size. In the larger grains, induced anisotropy is a function of grain orientation.     

Processing of magnetotelluric data

The processing of magnetotelluric data involves concepts from electromagnetic theory, time series analysis and linear systems theory for reducing natural electric and magnetic field variations recorded at the earth's surface to forms suitable for studying the electrical properties of the earth's interior.The electromagnetic field relations lead to either a scalar transfer impedance which couples an electric component to an orthogonal magnetic component at the surface of a plane-layered earth, or a tensor transfer impedance which couples each electric component to both magnetic components in the vicinity of a lateral inhomogeneity.A number of time series spectral analysis methods can be used for estimating the complex spectral coefficients of the various field quantities. These in turn are used for estimating the nature of the transfer function or tensor impedance. For two dimensional situations, the tensor impedance can be rotated to determine the principal directions of the electrical structure.In general for real data, estimates of the apparent resistivity are more stable when calculated from the tensor elements rather than from simple orthogonal field ratios (Cagniard estimates), even when the fields are measured in the principal coordinates.     

Determination of Stress State from Focal Mechanisms of Microseismic Events Induced During Hydraulic Injection at the Hijiori Hot Dry Rock Site

—?The stress state at the Hijiori hot dry rock site was estimated based on the inversion from focal mechanisms of microseismic events induced during hydraulic injection experiments. The best fit stress model obtained by inverting 58 focal mechanisms of seismic events simultaneously indicates that the maximum principal stress σ₁ is vertical, while the minimum principal stress σ₃ is horizontal and trends north-south. The average misfit between the stress model and all the data is 6.8°. The inversion results show that the average misfit is small enough to satisfy the assumption of homogeneity in the focal mechanism data and that the 95% confidence regions of σ₁ and σ₃ are well constrained, i.e., they do not overlap, suggesting that the inversion results are acceptable. The stress estimates obtained by the focal mechanism inversion essentially agree with other stress estimates previously obtained. It is therefore concluded that the focal mechanism inversion method provides a useful tool for estimating the stress state. The hypocentral distributions of microseismic events associated with the hydraulic fracturing experiments are distributed around the plane that spreads to almost east–west from the injection wells and declines to the north at a high angle. The vertical orientation and east–west strike of the seismic events are essentially coplanar with the caldera ring-fault structure in the southern portion of the Hijiori Caldera. This indicates that tensile fractures of intact rock were not being created, but pre-existing fractures were being re-opened and developed in the direction of the maximum horizontal principal stress, although microseismic events were caused by shear failures.     

Inversion of swell frequency from a 1-year HF radar dataset collected in Brittany (France)

This article presents long period ocean wave (swell) frequencies inverted from a 13-month dataset of high-frequency (HF) phased array radars and an assessment of these estimates by comparison with WAVEWATCH III model data. The method of swell frequency inversion from high-frequency radar sea echo Doppler spectra is described. Radar data were collected from a two-site HF Wellen Radar (WERA) radar system on the west coast of Brittany (France) operating at 12 MHz. A standard beam-forming processing technique has been used to obtain Doppler spectra of processed radar cells. Swell frequencies are obtained from the frequencies of particular spectral peaks of the second-order continuum in hourly averaged Doppler spectra. The data coverage of effective Doppler spectra considered for swell frequency estimates shows the influence of islands and shallow water effects. Swell estimates from both radar stations are in good agreement. The comparison of radar-derived results to WAVEWATCH III (WW3) estimates shows that radar measurements agree quite well with model results. The bias and standard deviation between two estimates are very small for swells with frequency less than 0.09 Hz (period >11 s), whereas radar estimates are generally lower than model estimates for shorter swells, along with higher standard deviation. Statistical analysis suggests that radar measurement uncertainty explains most of the difference between radar and model estimates. For each swell event, time series of frequency exhibits a quasi-linear frequency increase which is associated with the dispersive property of wave phase velocity. The use of swell frequency estimates from both radars on common radar cells only slightly increases the accuracy of swell frequency measurement.     

龙马溪页岩弹性各向异性与矿物分布之间的关系探讨

岩石弹性各向异性特征是普遍存在的,但导致岩石弹性各向异性的原因复杂且仍然存在一定争议.本研究以龙马溪页岩为例,试图建立页岩弹性各向异性和矿物分布之间关系.研究使用无损超声波探测获取岩石弹性各向异性参数,并使用背散射技术获取岩石矿物分布特征.研究通过引进变异系数来描述矿物或孔裂隙在不同方向的差异性,并通过2个正交方向的变异系数建立微观异质性指标,用于表征微观尺度上矿物或孔裂隙在不同方向的差异程度.微观背散射图像显示,龙马溪组页岩主要由石英和伊利石构成,且这两种矿物和孔裂隙在X和Y面上都有定向排列特征.相应地,它们的变异系数在X面和Y面上均表现出随角度增大而降低的特征;而在Z面,石英、伊利石和孔裂隙变异系数变化不明显,这与Z面上矿物和孔裂隙无明显方向性的特点一致.无损超声波探测结果显示,波速在X和Y面上随角度增加而减小,这与主要矿物和孔裂隙变异系数变化趋势相同;而在Z面,波速变化不大,与主要矿物和孔裂隙变异系数变化不明显的特征一致.以上观测结果说明,宏观波速与矿物的微观变异系数明显相关,暗示岩石弹性各向异性与矿物分布直接相关.     

Inversion of arrival-times in a region of dilatancy anisotropy

Shear-wave splitting has been identified in many three-component seismograms from two separate field experiments on a section of the North Anatolian Fault in North-West Turkey. These observations are consistent with shear-wave propagation through a zone of extensive-dilitancy anisotropy. A preliminary attempt has been made to confirm this interpretation by simultaneously inverting suites of arrival-times for hypocentral locations and for parameters describing an anisotropic halfspace. Although the inversion procedure is not globally convergent, it is possible to recognize the true solution by systematically varying the initial conditions. Applied to selected data sets, the inversion defines several anisotropic models that fit the data significantly better than a simple isotropic model, and display the anisotropy required by the shear-wave splitting. However, most of these anisotropic models are not superior when they are used to individually locate events in a much larger data set. However, for each experiment, there is a single model that produces clearly superior locations for the larger data sets than those of other anisotropic or simple isotropic models. Both models display similar velocity variations which are characteristic of propagation through distributions of biplanar cracks displaying orthorhombic symmetry. The principal axes of the two models are oriented in similar directions and are within 20° of the principal axis of regional stress derived from fault-plane solutions. The solutions indicate low velocities close to the tensional axis, as would be expected in extensive-dilatancy anisotropy.     

Anisotropy Versus Heterogeneity in Continental Solid Earth Electromagnetic Studies: Fundamental Response Characteristics and Implications for Physicochemical State

Electrical anisotropy, the effect of current density in a medium being a function of the orientation of the electric field, is being recognized increasingly as an important effect in explaining Earth electromagnetic observations. A consideration of anisotropy, however, in most cases is an admission of spatial aliasing in earth structure, wherein the averaging volume of diffusive EM fields may be greater than the characteristic dimensions of a family of oriented structures, thus leading to a response which is equivalent to a bulk anisotropic medium. Even for two-dimensional geometries, there can be strong non-parallelism of principal axes of vertical magnetic field relative to the impedance over broad areas, as well as impedance phase variations which leave normal quadrants, if there are multiple directions of anisotropy or anisotropy strike distinct from bulk geometric (2D) strike. This paper concentrates on experience with regional field studies in continental settings where bulk anisotropy is apparent. Upper crustal anisotropy may result from preferred orientations of fracture porosity, or lithologic layering, or oriented heterogeneity. Lower crustal anisotropy may result from preferred orientations of fluidized/melt-bearing or graphitized shear zones, but does not necessarily reflect current state of stress per se. In the upper mantle, the prior causes all may act in pertinent domains, but added to these is the possibility of strong electrical anisotropy due to hydrous defects within shear-aligned olivine crystals (solid-state conduction). Several field examples from continental MT investigations will be discussed, which roughly fall into active transpressional, active transtensional, and fossil transpressional regimes. A general challenge in interpreting data with apparent anisotropic effects is to establish the tradeoff between heterogeneity and anisotropy in the inversion of EM responses.     

Performance of the steady-state dipole flow test in layered aquifers

Steady-state numerical simulations of the dipole flow test in layered aquifers demonstrate that the test produces a good estimate of the equivalent hydraulic conductivity anisotropy ratio for the part of the aquifer spanned by the well chambers. The effects of chamber size, different conductivity of layers and layer location on the estimated anisotropy ratios are presented. The steady-state dipole flow test, when performed at different levels in the well, can yield estimates of the down-hole anisotropy ratio distribution. Numerical simulations also illustrate that the skin effect can significantly distort the anisotropy estimates produced by the dipole flow test. © 1997 John Wiley & Sons, Ltd.     

The effects of anisotropy of marine magnetic anomalies on the Curie point depth estimates from spectral analysis

Azimuthally averaged power spectra are widely used in the Curie point depth (CPD) estimation with the implicit assumption that the magnetization distribution is random and uncorrelated. However, the marine magnetic anomalies are caused by bands of normal and reverse magnetization and show obvious trends. To investigate the effects of the anisotropy of marine magnetic anomalies on the CPD estimates, we develop 3D fractal striped magnetization models to produce lineated marine magnetic anomalies for the first time. We analyze the spectra anisotropy of the lineated magnetic anomalies of the synthetic fractal striped magnetization models and investigate its effects on the CPD estimates. The synthetic models and actual data show that the spectra of the lineated marine magnetic anomalies are directionally anisotropic. The amplitude response is strong and the slope of the logarithmic spectrum is large in a direction perpendicular to the stripes of magnetic anomalies, whereas the amplitude response is weak and the slope of the logarithmic spectrum is small in a direction parallel to the stripes of magnetic anomalies. The depth estimates in the perpendicular direction are close to the actual values, whereas the depths estimates in the parallel direction are significantly lower than the actual values. The actual marine magnetic anomalies of the South China Sea exhibit an anisotropic power spectrum that is consistent with the spectral anisotropy of magnetic anomalies of the synthetic fractal striped magnetization models.     

Changes of anisotropy of the magnetic susceptibility of rocks induced by a magnetic field

Summary The changes of the anisotropy of magnetic susceptibility of igneous rocks, induced by a magnetic field, are studied. It is proved that changes in the degree of anisotropy of susceptibility and of the orientation of the susceptibility ellispoid of specimens occur due to the configuration of the domain structure under the effect of the magnetic field. The influence of this effect on the total anisotropy of rocks depends on the degree of anisotropy due to the shape factor and on the stability of the domain structure. A model concept is presented, explaining the qualitatively different pattern of the changes of the anisotropy of susceptibility under the effect of the magnetic field in various directions of the specimens.     

A case example of near-surface correction for multicomponent VSPS

A deterministic near-surface correction procedure is developed for multicomponent VSP data, shot using directional sources and recorded using three-component receivers. The method is capable of removing unwanted effects of acquisition such as unequal source strengths or misorientations, but may also remove near-surface multiples and anisotropy. This is of considerable benefit for obtaining accurate and consistent estimates of subsurface anisotropy from different source combinations. Application of the technique is illustrated using a dataset from the Romashkino reservoir in Russia, where three or four different directional sources are used at the same source locations. The technique corrects for the large discrepancies which exist between the estimates obtained using different source combinations. Application of the technique to three wells in the survey region reveals a nearly isotropic subsurface, except for a few isolated zones of moderate to high (2 to 8%) anisotropy which lie close to the expected depth for the reservoir. Although there is no significant correlation with the production figures for each well, the qS1 polarization azimuth within the reservoir does vary at each well location, suggesting that this may be a more sensitive indicator of reservoir Drocess.     

Sources of magnetic susceptibility in a slate

Magnetic susceptibility and its anisotropy in the Borrowdale Volcanic slates at Kentmere in the English Lake District are attributed largely to preferred orientation of a paramagnetic chlorite of diabantite-ripidolite composition. In units of 10⁻⁶ cgs/g, the principal susceptibilities for the slates are 9.61; 9.42; 8.69 and for the chlorite grains the minimum anisotropy is represented by principal susceptibilities of 11.57; 11.22 and 9.15. Because the magnetic susceptibility is carried by a tightly packed, matrix-forming mineral that has recrystallised during the deformation it is not possible to imagine simple grain rotation as being responsible for the anisotropy of susceptibility.     

磁铁矿的低温磁学性质研究进展

磁铁矿(Fe3O4)是自然界沉积物和岩石天然剩磁的最主要载磁矿物.本文主要评述低温条件下(2\|300K)磁铁矿的晶体结构特征、Verwey相变以及磁晶各向异性等,重点介绍了两种重要的剩磁随温度变化曲线.最后简单地介绍磁铁矿的低温磁学研究在地学中的应用以及氧化和超顺磁效应对低温磁学性质的影响.     

任意各向异性介质中三维可控源音频大地电磁正演模拟

在一些地层层理发育的地区,地下介质存在显著的电各向异性,此时基于各向同性模型解释含各向异性效应的可控源音频大地电磁（CSAMT）测深观测数据会导致错误的结果.本文通过引入3×3的对称正定张量表征电导率各向异性,采用非结构四面体网格和矢量有限元方法离散电场满足的矢量Helmholtz方程,并将电磁场源等效为系列电偶极子,实现任意各向异性介质中CSAMT高效数值模拟.本文首先通过层状各向异性模型检验三维有限元算法的精度和有效性,进一步建立三维地电模型研究异常体各向异性和围岩各向异性对CSAMT响应的影响,最后使用视电阻率极性图来识别各向异性电导率主轴方向.数值模拟结果表明,各向异性电导率对CSAMT视电阻率幅值及分布规律都有很大影响,视电阻率极性图能够很好地识别各向异性主轴方向.     

Scaling behavior and the effects of heterogeneity on shallow seismic imaging of mineral deposits: A case study from Brunswick No. 6 mining area,Canada

We have studied the scaling behavior of compressional-wave velocity and density logs from an exploration borehole that extends down to about 700 m depth in the Brunswick No. 6 mining area, Bathurst Mining Camp, Canada. Using statistical methods, vertical and horizontal scale lengths of heterogeneity were estimated. Vertical scale length estimates from the velocity, density and calculated acoustic impedance are 14 m, 33 m, and about 20 m, respectively. Although the estimated scale length for the acoustic impedance implies a weak scattering environment, elastic finite difference modeling of seismic wave propagation in 2D heterogeneous media demonstrates that even this weak scattering medium can mask seismic signals from small, but yet economically feasible, massive sulfide deposits. Further analysis of the synthetic seismic data suggests that in the presence of heterogeneity, lenticular-shaped targets may only exhibit incomplete diffraction signals whereby the down-dip tails of these diffractions are mainly visible on the stacked sections. Therefore, identification of orebody generated diffractions is much easier on the unmigrated stacked sections than on migrated stacked sections. The numerical seismic modeling in 2D heterogeneous media indicates that in the presence of large horizontal, but small vertical scale lengths (structural anisotropy), identification of massive sulfide deposits is possible, but their delineation at depth requires detailed velocity modeling and processing algorithms which can handle the anisotropy.     

Interpretation of transmissivity estimates from single-well pumping aquifer tests

Interpretation of single-well tests with the Cooper-Jacob method remains more reasonable than most alternatives. Drawdowns from 628 simulated single-well tests where transmissivity was specified were interpreted with the Cooper-Jacob straight-line method to estimate transmissivity. Error and bias as a function of vertical anisotropy, partial penetration, specific yield, and interpretive technique were investigated for transmissivities that ranged from 10 to 10,000 m(2)/d. Cooper-Jacob transmissivity estimates in confined aquifers were affected minimally by partial penetration, vertical anisotropy, or analyst. Cooper-Jacob transmissivity estimates of simulated unconfined aquifers averaged twice the known values. Transmissivity estimates of unconfined aquifers were not improved by interpreting results with an unconfined aquifer solution. Judicious interpretation of late-time data consistently improved estimates where transmissivity exceeded 250 m(2)/d in unconfined aquifers.