Structural delineation at the Los Humeros geothermal field, Mexico, by P-wave reflection retrieval from noise

Author(s):  
Arie Verdel ◽  
Boris Boullenger ◽  
Joana E. Martins ◽  
Anne Obermann ◽  
Tania Toledo ◽  
...  

<p>The overall purpose of the recently finalized GEMex project*, a European-Mexican collaboration, has been to gain an improved understanding of the subsurface at two unconventional geothermal sites: for EGS development at Acoculco and for a superhot resource near Los Humeros. Providing a more precise description of both the geological structure and the geothermal reservoir behavior for these two sites form important requirements for achieving that goal.</p><p>For delineating the main structural features at geothermal reservoir level, reflection retrieval from ambient seismic noise can be considered interesting because of its relatively low-cost and low environmental impact as compared to more conventional, controlled-source, seismic surveying practice, where (expensive) active sources are required.</p><p>In this study, we present results from the application of ambient noise seismic interferometry (ANSI) to retrieve zero-offset reflected P-waves from continuous seismic data recorded during the second half of 2017 at the Los Humeros geothermal field, Mexico. It is known from noise interferometry theory that reflected P-waves can provide local structural detail at locations directly underneath the employed seismic stations.</p><p>We address various data selection and processing aspects related to the retrieval of these reflected P-waves. The reflections are thereafter compared with modelled reflectivities at station locations with sufficient data availability, data quality and proximity to a location at which seismic interval velocity information is available from the literature.</p><p>From our study it can be concluded that the ANSI auto-correlation technique that was applied for zero-offset reflectivity retrieval at the Los Humeros site indeed can provide relatively high structural detail: for near-horizontal reflectors in the close vicinity of the selected stations, local depth-estimates of seismic velocity-contrasts were determined. This information can be used to constrain both the geological structure and geothermal reservoir property description.</p><p>As such, results from this passive-seismic method may partially complement and partially confirm subsurface information derived from active-seismic, that can only be acquired at a higher cost, which is more labor-intensive and which has more impact on the environment.</p><p>We thank the Mexican GEMex team around Angel Figueroa Soto from UMSNH and Marco Calo from UNAM for setting up the seismic network and station maintenance as well as data retrieval. The Comisión Federal de Electricidad (CFE) kindly provided us with access to their geothermal field and permission to install the seismic stations. OGS is thanked for providing us the location details of the four active seismic lines. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 727550 and the Mexican Energy Sustainability Fund CONACYT-SENER, project 2015-04-68074.</p><p>* http://www.gemex-h2020.eu/index.php?option=com_content&view=featured&Itemid=101&lang=en</p>

2019 ◽  
Vol 486 (2) ◽  
pp. 237-242
Author(s):  
I. P. Kuzin ◽  
L. I. Lobkovskiy ◽  
K. A. Dozorova

The results of coseismic GPS observations in the epicentral area of 2013 Sea-of- Okhotsk earthquake are presented and specific features of seismic waves amplitudes variations with distance are detected basing on the records of Russian and international seismic stations. Global propagation of P-waves for the Sea-of-Okhotsk and Bolivian (09.06.1994) earthquakes was studied and their amplitudes on teleseismic distances were estimated.


2012 ◽  
Vol 550-553 ◽  
pp. 2472-2477
Author(s):  
Yu Chun Bai ◽  
Yong Li Li ◽  
Fu Li Qi ◽  
Feng Long Zhang

Heiyu Lake zone of Daqing is located in the southwest hollow borderland of Heiyu Lake and on the arching transitional zone of Daqing placanticline. Based on the geological background of Heiyu Lake, this paper analyzes the landform, the regional geological structure, the formation lithology and the irruptive rock and other metallogenic conditions in detail. The indispensable geological conditions for forming geothermal field in layers were summed up. Combining with the development characteristics and geophysical data of formation, the bore hole site of geothermal well and target stratum were ascertained. The four major elements of forming geothermal resources in this region were confirmed by carrying out geothermal drilling.


Geophysics ◽  
1988 ◽  
Vol 53 (1) ◽  
pp. 76-84 ◽  
Author(s):  
E. L. Majer ◽  
T. V. McEvilly ◽  
F. S. Eastwood ◽  
L. R. Myer

In a pilot vertical seismic profiling study, P-wave and cross‐polarized S-wave vibrators were used to investigate the potential utility of shear‐wave anisotropy measurements in characterizing a fractured rock mass. The caprock at The Geysers geothermal field was found to exhibit about an 11 percent velocity variation between SH-waves and SV-waves generated by rotating the S-wave vibrator orientation to two orthogonal polarizations for each survey level in the well. The effect is generally consistent with the equivalent anisotropy expected from the known fracture geometry.


Geophysics ◽  
2021 ◽  
Vol 86 (3) ◽  
pp. T155-T164
Author(s):  
Wanting Hou ◽  
Li-Yun Fu ◽  
José M. Carcione ◽  
Zhiwei Wang ◽  
Jia Wei

Thermoelasticity is important in seismic propagation due to the effects related to wave attenuation and velocity dispersion. We have applied a novel finite-difference (FD) solver of the Lord-Shulman thermoelasticity equations to compute synthetic seismograms that include the effects of the thermal properties (expansion coefficient, thermal conductivity, and specific heat) compared with the classic forward-modeling codes. We use a time splitting method because the presence of a slow quasistatic mode (the thermal mode) makes the differential equations stiff and unstable for explicit time-stepping methods. The spatial derivatives are computed with a rotated staggered-grid FD method, and an unsplit convolutional perfectly matched layer is used to absorb the waves at the boundaries, with an optimal performance at the grazing incidence. The stability condition of the modeling algorithm is examined. The numerical experiments illustrate the effects of the thermoelasticity properties on the attenuation of the fast P-wave (or E-wave) and the slow thermal P-wave (or T-wave). These propagation modes have characteristics similar to the fast and slow P-waves of poroelasticity, respectively. The thermal expansion coefficient has a significant effect on the velocity dispersion and attenuation of the elastic waves, and the thermal conductivity affects the relaxation time of the thermal diffusion process, with the T mode becoming wave-like at high thermal conductivities and high frequencies.


1991 ◽  
Vol 81 (2) ◽  
pp. 508-523
Author(s):  
Jim Mori

Abstract Event record sections, which are constructed by plotting seismograms from many closely spaced earthquakes recorded on a few stations, show multiple free-surface reflections (PP, PPP, PPPP) of the P wave in the Imperial Valley, California. The relative timing of these arrivals is used to estimate the strength of the P-wave velocity gradient within the upper 5 km of the sediment layer. Consistent with previous studies, a velocity model with a value of 1.8 km/sec at the surface increasing linearly to 5.8 km/sec at a depth of 5.5 km fits the data well. The relative amplitudes of the P and PP arrivals are used to estimate the source depth for the aftershock distributions of the Elmore Ranch and Superstition Hills main shocks. Although the depth determination has large uncertainties, both the Elmore Ranch and Superstition Hills aftershock sequences appear to have similar depth distribution in the range of 4 to 10 km.


Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. C175-C185 ◽  
Author(s):  
Ivan Pšenčík ◽  
Véronique Farra

We have developed approximate nonhyperbolic P-wave moveout formulas applicable to weakly or moderately anisotropic media of arbitrary anisotropy symmetry and orientation. Instead of the commonly used Taylor expansion of the square of the reflection traveltime in terms of the square of the offset, we expand the square of the reflection traveltime in terms of weak-anisotropy (WA) parameters. No acoustic approximation is used. We specify the formulas designed for anisotropy of arbitrary symmetry for the transversely isotropic (TI) media with the axis of symmetry oriented arbitrarily in the 3D space. Resulting formulas depend on three P-wave WA parameters specifying the TI symmetry and two angles specifying the orientation of the axis of symmetry. Tests of the accuracy of the more accurate of the approximate formulas indicate that maximum relative errors do not exceed 0.3% or 2.5% for weak or moderate P-wave anisotropy, respectively.


2012 ◽  
Vol 56 (4) ◽  
pp. 631-635 ◽  

Abstract The electrocardiographic examination was performed in 33 training horses (2-16 years of age, 11 males and 22 females). Einthoven and precordial leads (I, II, III, aVR, aVL, aVF, CV1, CV2, CV4) were used. The ECG was performed in resting horses and immediately after exercise (10 min walk, 15 min trot, 10 min canter) using a portable Schiller AT-1 3-channel electrocardiograph, with a paper speed of 25 mm sec-1 and a sensitivity of 10 mm.mV-1. The heart rate, wave amplitudes, and duration time were estimated manually. All horses presented a significant increase in heart rate after exercise (rest 43.83 ±10.33 vs. exercise 73.2 ±14.8). QT intervals were significantly shortened in most of the leads. In resting horses, all P waves in the lead I were positive and almost all II, III and CV4 leads were positive. Simple negative P wave dominated in aVR and only simple negative T wave was found in the leads I. The biphasic shape was observed. After exercise, the amplitude of P and T waves rose, however, clear changes were not observed in wave polarisation and form. In the absence of specific racial characteristics of the electrocardiogram in the Polish Anglo- Arabians, electrocardiographic findings can be interpreted according to ECG standards adopted for horses.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. C171-C180 ◽  
Author(s):  
Qifan Liu ◽  
Ilya Tsvankin

Tilted orthorhombic (TOR) models are typical for dipping anisotropic layers, such as fractured shales, and can also be due to nonhydrostatic stress fields. Velocity analysis for TOR media, however, is complicated by the large number of independent parameters. Using multicomponent wide-azimuth reflection data, we develop stacking-velocity tomography to estimate the interval parameters of TOR media composed of homogeneous layers separated by plane dipping interfaces. The normal-moveout (NMO) ellipses, zero-offset traveltimes, and reflection time slopes of P-waves and split S-waves ([Formula: see text] and [Formula: see text]) are used to invert for the interval TOR parameters including the orientation of the symmetry planes. We show that the inversion can be facilitated by assuming that the reflector coincides with one of the symmetry planes, which is a common geologic constraint often employed for tilted transversely isotropic media. This constraint makes the inversion for a single TOR layer feasible even when the initial model is purely isotropic. If the dip plane is also aligned with one of the symmetry planes, we show that the inverse problem for [Formula: see text]-, [Formula: see text]-, and [Formula: see text]-waves can be solved analytically. When only [Formula: see text]-wave data are available, parameter estimation requires combining NMO ellipses from a horizontal and dipping interface. Because of the increase in the number of independent measurements for layered TOR media, constraining the reflector orientation is required only for the subsurface layer. However, the inversion results generally deteriorate with depth because of error accumulation. Using tests on synthetic data, we demonstrate that additional information such as knowledge of the vertical velocities (which may be available from check shots or well logs) and the constraint on the reflector orientation can significantly improve the accuracy and stability of interval parameter estimation.


Diagnostics ◽  
2021 ◽  
Vol 11 (9) ◽  
pp. 1694
Author(s):  
Dimitrios Tachmatzidis ◽  
Dimitrios Filos ◽  
Ioanna Chouvarda ◽  
Anastasios Tsarouchas ◽  
Dimitrios Mouselimis ◽  
...  

Early identification of patients at risk for paroxysmal atrial fibrillation (PAF) is essential to attain optimal treatment and a favorable prognosis. We compared the performance of a beat-to-beat (B2B) P-wave analysis with that of standard P-wave indices (SPWIs) in identifying patients prone to PAF. To this end, 12-lead ECG and 10 min vectorcardiogram (VCG) recordings were obtained from 33 consecutive, antiarrhythmic therapy naïve patients, with a short history of low burden PAF, and from 56 age- and sex-matched individuals with no AF history. For both groups, SPWIs were calculated, while the VCG recordings were analyzed on a B2B basis, and the P-waves were classified to a primary or secondary morphology. Wavelet transform was used to further analyze P-wave signals of main morphology. Univariate analysis revealed that none of the SPWIs performed acceptably in PAF detection, while five B2B features reached an AUC above 0.7. Moreover, multivariate logistic regression analysis was used to develop two classifiers—one based on B2B analysis derived features and one using only SPWIs. The B2B classifier was found to be superior to SPWIs classifier; B2B AUC: 0.849 (0.754–0.917) vs. SPWIs AUC: 0.721 (0.613–0.813), p value: 0.041. Therefore, in the studied population, the proposed B2B P-wave analysis outperforms SPWIs in detecting patients with PAF while in sinus rhythm. This can be used in further clinical trials regarding the prognosis of such patients.


2022 ◽  
Vol 9 ◽  
Author(s):  
Yifan Fan ◽  
Shikuan Zhang ◽  
Yonghui Huang ◽  
Zhonghe Pang ◽  
Hongyan Li

Recoverable geothermal resources are very important for geothermal development and utilization. Generally, the recovery factor is a measure of available geothermal resources in a geothermal field. However, it has been a pre-determined ratio in practice and sustainable utilization of geothermal resources was not considered in the previous calculation of recoverable resources. In this work, we have attempted to develop a method to calculate recoverable geothermal resources based on a numerical thermo-hydraulic coupled modeling of a geothermal reservoir under exploitation, with an assumption of sustainability. Taking a geothermal reservoir as an example, we demonstrate the effectiveness of the method. The recoverable geothermal resources are 6.85 × 1018 J assuming a lifetime of 100 years in a well doublet pattern for geothermal heating. We further discuss the influence of well spacing on the recoverable resources. It is found that 600 m is the optimal well spacing with maximum extracted energy that conforms to the limit of the pressure drop and no temperature drop in the production well. Under the uniform well distribution pattern for sustainable exploitation, the recovery factor is 26.2%, which is higher than the previous value of 15% when depending only on lithology. The proposed method for calculating the recoverable geothermal resources is instructive for making decisions for sustainable exploitation.


Sign in / Sign up

Export Citation Format

Share Document