scholarly journals Seismic Anisotropy and  Stress of the Canterbury  Plains

2021 ◽  
Author(s):  
◽  
Rob Holt

<p>The Mѡ=7.1 Darfield (Canterbury) earthquake struck on 4 September 2010, approximately 45 km west of Christchurch, New Zealand. It revealed a previously unknown fault (the Greendale fault) and caused billions of dollars of damage due to high peak ground velocities and extensive liquefaction. It also triggered the Mw=6.3 Christchurch earthquake on 22 February 2011, which caused further damage and the loss of 185 lives. The objective of this research was to determine the relationship between stress and seismic properties in a seismically active region using manually-picked P and S wave arrival times from the aftershock sequence between 8 September 2010-13 January 2011 to estimate shear-wave splitting (SWS) parameters, VP =VS-ratios, anisotropy (delay-time tomography), focal mechanisms, and tectonic stress on the Canterbury plains. The maximum horizontal stress direction was highly consistent in the plains, with an average value of SHmax=116 18 . However, the estimates showed variation in SHmax near the fault, with one estimate rotating by as much as 30° counter-clockwise. This suggests heterogeneity of stress at the fault, though the cause remains unclear. Orientations of the principal stresses predominantly indicate a strike-slip regime, but there are possible thrust regimes to the west and north/east of the fault. The SWS fast directions (ø) on the plains show alignment with SHmax at the majority of stations, indicating stress controlled anisotropy. However, structural effects appear more dominant in the neighbouring regions of the Southern Alps and Banks Peninsula.</p>

2021 ◽  
Author(s):  
◽  
Rob Holt

<p>The Mѡ=7.1 Darfield (Canterbury) earthquake struck on 4 September 2010, approximately 45 km west of Christchurch, New Zealand. It revealed a previously unknown fault (the Greendale fault) and caused billions of dollars of damage due to high peak ground velocities and extensive liquefaction. It also triggered the Mw=6.3 Christchurch earthquake on 22 February 2011, which caused further damage and the loss of 185 lives. The objective of this research was to determine the relationship between stress and seismic properties in a seismically active region using manually-picked P and S wave arrival times from the aftershock sequence between 8 September 2010-13 January 2011 to estimate shear-wave splitting (SWS) parameters, VP =VS-ratios, anisotropy (delay-time tomography), focal mechanisms, and tectonic stress on the Canterbury plains. The maximum horizontal stress direction was highly consistent in the plains, with an average value of SHmax=116 18 . However, the estimates showed variation in SHmax near the fault, with one estimate rotating by as much as 30° counter-clockwise. This suggests heterogeneity of stress at the fault, though the cause remains unclear. Orientations of the principal stresses predominantly indicate a strike-slip regime, but there are possible thrust regimes to the west and north/east of the fault. The SWS fast directions (ø) on the plains show alignment with SHmax at the majority of stations, indicating stress controlled anisotropy. However, structural effects appear more dominant in the neighbouring regions of the Southern Alps and Banks Peninsula.</p>


2013 ◽  
Vol 5 (2) ◽  
pp. 1191-1257 ◽  
Author(s):  
T. Satsukawa ◽  
B. Ildefonse ◽  
D. Mainprice ◽  
L. F. G. Morales ◽  
K. Michibayashi ◽  
...  

Abstract. This study presents a unique database of 170 plagioclase Crystallographic Preferred Orientations (CPO) of variously deformed gabbroic rocks. The CPO characteristics as a function of the deformation regime (magmatic or crystal-plastic) are outlined and discussed. The studied samples are dominantly from slow- and fast-spread present-day ocean crust, as well as from the Oman ophiolite. Plagioclase is the dominant mineral phase in the studied samples. Plagioclase CPOs are grouped in three main categories: Axial-B, a strong point alignment of (010) with a girdle distribution of [100]; Axial-A, a strong point maximum concentration of [100] with parallel girdle distributions of (010) and (001); and P-type, point maxima of [100], (010), and (001). A majority of CPO patterns are Axial-B and P-type, in samples showing either magmatic or crystal-plastic deformation textures. Axial-A CPOs are less common; they represent 21% of the samples deformed by crystal-plastic flow. Although fabric strength (ODF J-index) does not show any consistent variation as a function of the CPO patterns, there is a significant difference in the relationship between the ODF and pole figures J-indices; the magmatic type microstructures have high (010) pole figures J-indices, which increase linearly with ODF J-index, whereas the high [100] pole figures J-indices of plastically deformed samples vary in a more scattered manner with ODF J-index. The multistage nature of plastic deformation superposed on a magmatic structure compared with magmatic flow, and the large number of possible slip-systems in plagioclase probably account for these differences. Calculated seismic properties (P wave and S wave velocities and anisotropies) of plagioclase aggregates show that anisotropy (up to 12% for P wave and 14% for S wave) tends to increase as a function of ODF J-index. In comparison with the olivine 1998 CPO database, the magnitude of P wave anisotropy for a given J-index is much less than olivine, whereas it is similar for S wave anisotropy. Despite a large variation of fabric patterns and geodynamic setting, seismic properties of plagioclase-rich rocks have similar magnitudes of anisotropy. There is a small difference in the aggregate elastic symmetry, with magmatic microstructures having higher orthorhombic and hexagonal components, whereas plastic deformation microstructures have a slightly higher monoclinic component, possibly correlated with predominant monoclinic simple shear flow in plastically-deformed samples. Overall, plots for CPO strength (ODF J-index), pole figure strength, CPO symmetry and seismic anisotropy show significant scattering. This could be related to sampling statistics, although our database is a factor of ten higher than the olivine database of 1998, or it could be related to the low symmetry (triclinic) structure of plagioclase resulting in the addition of degrees of freedom in the processes creating the CPOs.


Solid Earth ◽  
2014 ◽  
Vol 5 (1) ◽  
pp. 45-63 ◽  
Author(s):  
V. Baptiste ◽  
A. Tommasi

Abstract. We calculated the seismic properties of 47 mantle xenoliths from 9 kimberlitic pipes in the Kaapvaal craton based on their modal composition, the crystal-preferred orientations (CPO) of olivine, ortho- and clinopyroxene, and garnet, the Fe content of olivine, and the pressures and temperatures at which the rocks were equilibrated. These data allow constraining the variation of seismic anisotropy and velocities within the cratonic mantle. The fastest P and S2 wave propagation directions and the polarization of fast split shear waves (S1) are always subparallel to olivine [100] axes of maximum concentration, which marks the lineation (fossil flow direction). Seismic anisotropy is higher for high olivine contents and stronger CPO. Maximum P wave azimuthal anisotropy (AVp) ranges between 2.5 and 10.2% and the maximum S wave polarization anisotropy (AVs), between 2.7 and 8%. Changes in olivine CPO symmetry result in minor variations in the seismic anisotropy patterns, mainly in the apparent isotropy directions for shear wave splitting. Seismic properties averaged over 20 km-thick depth sections are, therefore, very homogeneous. Based on these data, we predict the anisotropy that would be measured by SKS, Rayleigh (SV) and Love (SH) waves for five endmember orientations of the foliation and lineation. Comparison to seismic anisotropy data from the Kaapvaal shows that the coherent fast directions, but low delay times imaged by SKS studies, and the low azimuthal anisotropy with with the horizontally polarized S waves (SH) faster than the vertically polarized S wave (SV) measured using surface waves are best explained by homogeneously dipping (45°) foliations and lineations in the cratonic mantle lithosphere. Laterally or vertically varying foliation and lineation orientations with a dominantly NW–SE trend might also explain the low measured anisotropies, but this model should also result in backazimuthal variability of the SKS splitting data, not reported in the seismological data. The strong compositional heterogeneity of the Kaapvaal peridotite xenoliths results in up to 3% variation in density and in up to 2.3% variation of Vp, Vs, and Vp / Vs ratio. Fe depletion by melt extraction increases Vp and Vs, but decreases the Vp / Vs ratio and density. Orthopyroxene enrichment due to metasomatism decreases the density and Vp, strongly reducing the Vp / Vs ratio. Garnet enrichment, which was also attributed to metasomatism, increases the density, and in a lesser extent Vp and the Vp / Vs ratio. Comparison of density and seismic velocity profiles calculated using the xenoliths' compositions and equilibration conditions to seismological data in the Kaapvaal highlights that (i) the thickness of the craton is underestimated in some seismic studies and reaches at least 180 km, (ii) the deep sheared peridotites represent very local modifications caused and oversampled by kimberlites, and (iii) seismological models probably underestimate the compositional heterogeneity in the Kaapvaal mantle root, which occurs at a scale much smaller than the one that may be sampled seismologically.


2021 ◽  
Vol 2021 ◽  
pp. 1-13
Author(s):  
Ahmed Gowida ◽  
Ahmed Farid Ibrahim ◽  
Salaheldin Elkatatny ◽  
Abdulwahab Ali

The least principal stresses of downhole formations include minimum horizontal stress (σmin) and maximum horizontal stress (σmax). σmin and σmax are substantial parameters that significantly affect the design and optimization of the drilling process. These stresses can be estimated using theoretical equations in addition to some field tests, i.e., leak-off test to include the effect of tectonic stress. This approach is associated with many technical and financial issues. Therefore, the objective of this study is to provide a novel machine learning-based solution to estimate these stresses while drilling. First, new models were developed using artificial neural network (ANN) to directly predict σmin and σmax from the drilling data; which are injection rate (Q), standpipe pressure (SPP), weight on bit (WOB), torque (T), and rate of penetration (ROP). Such data are always available while drilling, and hence, no additional cost is required. Actual data from a Middle Eastern field were collected, statistically analyzed, and fed to the models. First, the models’ predictions showed a significant match with the actual stress values with a correlation coefficient (R-value) exceeding 0.90 and a mean absolute average error (MAPE) of 0.75% as a maximum. Second, new empirical equations were generated based on the developed ANN-based models. The new equations were then validated using another unseen dataset from the same field. The predictions had an R-value of 0.98 and 0.93 in addition to MAPE of 0.36% and 0.96% for σmin and σmax models, respectively. The results demonstrated the outperformance of the developed ANN-based equations to estimate the least principal stresses from the drilling data with high accuracy in a timely and economically effective way.


Solid Earth ◽  
2013 ◽  
Vol 4 (2) ◽  
pp. 511-542 ◽  
Author(s):  
T. Satsukawa ◽  
B. Ildefonse ◽  
D. Mainprice ◽  
L. F. G. Morales ◽  
K. Michibayashi ◽  
...  

Abstract. This study presents a unique database of 172 plagioclase Crystallographic Preferred Orientations (CPO) of variously deformed gabbroic rocks. The CPO characteristics as a function of the deformation regime (magmatic or crystal-plastic) are outlined and discussed. The studied samples are dominantly from slow- and fast-spread present-day ocean crust, as well as from the Oman ophiolite. Plagioclase is the dominant mineral phase in the studied samples. Plagioclase CPOs are grouped into three main categories: Axial-B, a strong point alignment of (010) with a girdle distribution of [100]; Axial-A, a strong point maximum concentration of [100] with parallel girdle distributions of (010) and (001); and P-type, point maxima of [100], (010), and (001). A majority of CPO patterns are Axial-B and P-type, in samples showing either magmatic or crystal-plastic deformation textures. Axial-A CPOs are less common; they represent 21% of the samples deformed by crystal-plastic flow. Although fabric strength (ODF J index) does not show any consistent variation as a function of the CPO patterns, there is a significant difference in the relationship between the ODF and pole figures J indices; the magmatic type microstructures have high (010) pole figures J indices, which increase linearly with ODF J index, whereas the high [100] pole figures J indices of plastically deformed samples vary in a more scattered manner with ODF J index. The multistage nature of plastic deformation superposed on a magmatic structure compared with magmatic flow, and the large number of possible slip-systems in plagioclase probably account for these differences. Calculated seismic properties (P wave and S wave velocities and anisotropies) of plagioclase aggregates show that anisotropy (up to 12% for P wave and 14% for S wave) tends to increase as a function of ODF J index. In comparison with the olivine 1998 CPO database, the magnitude of P wave anisotropy for a given J index is much less than olivine, whereas it is similar for S wave anisotropy. Despite a large variation of fabric patterns and geodynamic setting, seismic properties of plagioclase-rich rocks have similar magnitudes of anisotropy. There is a small difference in the aggregate elastic symmetry, with magmatic microstructures having higher orthorhombic and hexagonal components, whereas plastic deformation microstructures have a slightly higher monoclinic component, possibly correlated with predominant monoclinic simple shear flow in plastically deformed samples. Overall, plots for CPO strength (ODF J index), pole figure strength, CPO symmetry and seismic anisotropy show significant scattering. This could be related to sampling statistics, although our database is a factor of ten higher than the olivine database of 1998, or it could be related to the low symmetry (triclinic) structure of plagioclase resulting in the addition of degrees of freedom in the processes creating the CPOs.


Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D245-D261 ◽  
Author(s):  
Jaime Meléndez-Martínez ◽  
Douglas R. Schmitt

We obtained the complete set of dynamic elastic stiffnesses for a suite of “shales” representative of unconventional reservoirs from simultaneously measured P- and S-wave speeds on single prisms specially machined from cores. Static linear compressibilities were concurrently obtained using strain gauges attached to the prism. Regardless of being from static or dynamic measurements, the pressure sensitivity varies strongly with the direction of measurement. Furthermore, the static and dynamic linear compressibilities measured parallel to the bedding are nearly the same whereas those perpendicular to the bedding can differ by as much as 100%. Compliant cracklike porosity, seen in scanning electron microscope images, controls the elastic properties measured perpendicular to the rock’s bedding plane and results in highly nonlinear pressure sensitivity. In contrast, those properties measured parallel to the bedding are nearly insensitive to stress. This anisotropy to the pressure dependency of the strains and moduli further complicates the study of the overall anisotropy of such rocks. This horizontal stress insensitivity has implications for the use of advanced sonic logging techniques for stress direction indication. Finally, we tested the validity of the practice of estimating the fracture pressure gradient (i.e., horizontal stress) using our observed elastic engineering moduli and found that ignoring anisotropy would lead to underestimates of the minimum stress by as much as 90%. Although one could ostensibly obtain better values or the minimum stress if the rock anisotropy is included, we would hope that these results will instead discourage this method of estimating horizontal stress in favor of more reliable techniques.


2020 ◽  
Author(s):  
Felix Kästner ◽  
Simona Pierdominici ◽  
Judith Elger ◽  
Christian Berndt ◽  
Alba Zappone ◽  
...  

&lt;p&gt;Deeply rooted thrust zones are key features of tectonic processes and the evolution of mountain belts. Exhumed and deeply-eroded orogens like the Scandinavian Caledonides allow to study such systems from the surface. Previous seismic investigations of the Seve Nappe Complex have shown indications for a strong but discontinuous reflectivity of this thrust zone, which is only poorly understood. The correlation of seismic properties measured on borehole cores with surface seismic data can help to constrain the origin of this reflectivity. In this study, we compare seismic velocities measured on cores to in situ velocities measured in the borehole. The core and downhole velocities deviate by up to 2 km/s. However, velocities of mafic rocks are generally in close agreement. Seismic anisotropy increases from about 5 to 26 % at depth, indicating a transition from gneissic to schistose foliation. Differences in the core and downhole velocities are most likely the result of microcracks due to depressurization of the cores. Thus, seismic velocity can help to identify mafic rocks on different scales whereas the velocity signature of other lithologies is obscured in core-derived velocities. Metamorphic foliation on the other hand has a clear expression in seismic anisotropy. To further constrain the effects of mineral composition, microstructure and deformation on the measured seismic anisotropy, we conducted additional microscopic investigations on selected core samples. These analyses using electron-based microscopy and X-ray powder diffractometry indicate that the anisotropy is strongest for mica schists followed by amphibole-rich units. This also emphasizes that seismic velocity and anisotropy are of complementary importance to better distinguish the present lithological units. Our results will aid in the evaluation of core-derived seismic properties of high-grade metamorphic rocks at the COSC-1 borehole and elsewhere.&lt;/p&gt;


2015 ◽  
Author(s):  
Fen Yang ◽  
Larry K. Britt ◽  
Shari Dunn-Norman

Abstract Since the late 1980's when Maersk published their work on multiple fracturing of horizontal wells in the Dan Field, the use of transverse multiple fractured horizontal wells has become the completion of choice and become the “industry standard” for unconventional and tight oil and tight gas reservoirs. Today approximately sixty percent of all wells drilled in the United States are drilled horizontally and nearly all of them are multiple fractured. Because a horizontal well adds additional cost and complexity to the drilling, completion, and stimulation of the well we need to fully understand anything that affects the cost and complexity. In other words, we need to understand the affects of the principal stresses, both direction and magnitude, on the drilling completion, and stimulation of these wells. However, little work has been done to address and understand the relationship between the principal stresses and the lateral direction. This paper has as its goal to fundamentally address the question, in what direction should I drill my lateral? Do I drill it in the direction of the maximum horizontal stress (longitudinal) or do I drill it in the direction of the minimum horizontal stress (transverse)? The answer to this question relates directly back to the title of this paper and please "Don't let your land man drive that decision." This paper focuses on the horizontal well's lateral direction (longitudinal or transverse fracture orientation) and how that direction influences productivity, reserves, and economics of horizontal wells. Optimization studies using a single phase fully three dimensional numeric simulator including convergent non-Darcy flow were used to highlight the importance of lateral direction as a function of reservoir permeability. These studies, conducted for both oil and gas, are used to identify the point on the permeability continuum where longitudinal wells outperform transverse wells. The simulations compare and contrast the transverse multiple fractured horizontal well to longitudinal wells based on the number of fractures and stages. Further, the effects of lateral length, fracture half-length, and fracture conductivity were investigated to see how these parameters affected the decision over lateral direction in both oil and gas reservoirs. Additionally, how does completion style affect the lateral direction? That is, how does an open hole completion compare to a cased hole completion and should the type of completion affect the decision on in what direction the lateral should be drilled? These simulation results will be used to discuss the various horizontal well completion and stimulation metrics (rate, recovery, and economics) and how the choice of metrics affects the choice of lateral direction. This paper will also show a series of field case studies to illustrate actual field comparisons in both oil and gas reservoirs of longitudinal versus transverse horizontal wells and tie these field examples and results to the numeric simulation study. This work benefits the petroleum industry by: Establishing well performance and economic based criteria as a function of permeability for drilling longitudinal or transverse horizontal wells,Integrating the reservoir objectives and geomechanic limitations into a horizontal well completion and stimulation strategy,Developing well performance and economic objectives for horizontal well direction (transverse versus longitudinal) and highlighting the incremental benefits of various completion and stimulation strategies.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. B183-B195 ◽  
Author(s):  
K. De Meersman ◽  
J.-M. Kendall ◽  
M. van der Baan

We relocate 303 microseismic events recorded in 1998 by sensors in a single borehole in the North Sea Valhall oil field. A semiautomated array analysis method repicks the P- and S-wave arrival times and P-wave polarizations, which are needed to locate these events. The relocated sources are confined predominantly to a [Formula: see text]-thick zone just above the reservoir, and location uncertainties are half those of previous efforts. Multiplet analysis identifies 40 multiplet groups, which include 208 of the 303 events. The largest group contains 24 events, and five groups contain 10 or more events. Within each multiplet group, we further improve arrival-time picking through crosscorrelation, which enhances the relative accuracy of the relocated events and reveals that more than 99% of the seismic activity lies spatially in three distinct clusters. The spatial distribution of events and wave-form similarities reveal two faultlike structures that match well with north-northwest–south-southeast-trending fault planes interpreted from 3D surface seismic data. Most waveform differences between multiplet groups located on these faults can be attributed to S-wave phase content and polarity or P-to-S amplitude ratio. The range in P-to-S amplitude ratios observed on the faults is explained best in terms of varying source mechanisms. We also find a correlation between multiplet groups and temporal variations in seismic anisotropy, as revealed by S-wave splitting analysis. We explain these findings in the context of a cyclic recharge and dissipation of cap-rock stresses in response to production-driven compaction of the underlying oil reservoir. The cyclic nature of this mechanism drives the short-term variations in seismic anisotropy and the reactivation of microseismic source mechanisms over time.


Sign in / Sign up

Export Citation Format

Share Document