Seismic Anisotropy at the Rotokawa and Ngatamariki Geothermal Fields in the Taupo Volcanic Zone

<p>In order to investigate the cracks/fractures in the geothermal fields of Rotokawa and Ngatamariki, we measure seismic anisotropy across both fields and interpret the results in the context of stress aligned microcracks. Cracks aligned perpendicular to the direction of maximum horizontal stress close and their fluid is forced into cracks aligned with maximum horizontal stress (SHmax). Seismic anisotropy is the directional dependence of a seismic wave's velocity and provides a measure of crack orientation and density. To measure seismic anisotropy we conduct shear wave splitting measurements on 52,000 station-earthquake pairs across both Rotokawa and Ngatamariki from earthquakes recorded during 2015. Both fields are the subject of other geophysical and geological studies. Thus they are excellent subjects for studying seismic anisotropy. We cluster our measurements by their station-event path and fit the parameters from these clusters to those from theoretical crack planes. We also apply 2-D tomography to shear wave splitting time delays (𝛿t) and spatial averaging to shear wave splitting fast polarisations (∅). In addition, we compare time delays with P-wave to S-wave velocity ratios (νP / vS). Local measurements of stress within Rotokawa and regional measures of stress within the Taupo Volcanic Zone provide a comparison for the shear wave splitting measurements. We measure ∅ which agrees with the NE-SW regional direction of SHmax across Ngatamariki and parts of Rotokawa. Within Rotokawa, we observe a rotation of ∅ away from NE-SW toward N-S that agrees with borehole measurements of direction of SHmax of 023° and 030°. Spatial averaging of ∅ reveals mean orientations close to the strike of nearby active faults. The theoretical crack planes, that fit best to the shear wave splitting measurements, correspond to aligned cracks striking 045° outside of both fields, 035° within Ngatamariki, and 035° through to 0° within Rotokawa. The average percent anisotropy for the full dataset, approximately 4%, is close to the upper bound for an intact rock. Delay time tomography shows regions of higher delay time per kilometre of path length (s=km) within both fields and possibly associated with the production field fault in Rotokawa. vP =vS shows a wide range of normally distributed values, from 1.1 through to 2.4 with a mean of 1.6, indicating a mixture of gas filled and saturated cracks. A positive correlation between delay time per kilometre (𝛿tpkm) and νP /νS indicates that the majority of the cracks are saturated.</p>

Seismic Anisotropy at the Rotokawa and Ngatamariki Geothermal Fields in the Taupo Volcanic Zone

<p>In order to investigate the cracks/fractures in the geothermal fields of Rotokawa and Ngatamariki, we measure seismic anisotropy across both fields and interpret the results in the context of stress aligned microcracks. Cracks aligned perpendicular to the direction of maximum horizontal stress close and their fluid is forced into cracks aligned with maximum horizontal stress (SHmax). Seismic anisotropy is the directional dependence of a seismic wave's velocity and provides a measure of crack orientation and density. To measure seismic anisotropy we conduct shear wave splitting measurements on 52,000 station-earthquake pairs across both Rotokawa and Ngatamariki from earthquakes recorded during 2015. Both fields are the subject of other geophysical and geological studies. Thus they are excellent subjects for studying seismic anisotropy. We cluster our measurements by their station-event path and fit the parameters from these clusters to those from theoretical crack planes. We also apply 2-D tomography to shear wave splitting time delays (𝛿t) and spatial averaging to shear wave splitting fast polarisations (∅). In addition, we compare time delays with P-wave to S-wave velocity ratios (νP / vS). Local measurements of stress within Rotokawa and regional measures of stress within the Taupo Volcanic Zone provide a comparison for the shear wave splitting measurements. We measure ∅ which agrees with the NE-SW regional direction of SHmax across Ngatamariki and parts of Rotokawa. Within Rotokawa, we observe a rotation of ∅ away from NE-SW toward N-S that agrees with borehole measurements of direction of SHmax of 023° and 030°. Spatial averaging of ∅ reveals mean orientations close to the strike of nearby active faults. The theoretical crack planes, that fit best to the shear wave splitting measurements, correspond to aligned cracks striking 045° outside of both fields, 035° within Ngatamariki, and 035° through to 0° within Rotokawa. The average percent anisotropy for the full dataset, approximately 4%, is close to the upper bound for an intact rock. Delay time tomography shows regions of higher delay time per kilometre of path length (s=km) within both fields and possibly associated with the production field fault in Rotokawa. vP =vS shows a wide range of normally distributed values, from 1.1 through to 2.4 with a mean of 1.6, indicating a mixture of gas filled and saturated cracks. A positive correlation between delay time per kilometre (𝛿tpkm) and νP /νS indicates that the majority of the cracks are saturated.</p>

Characterizing microseismicity at the Rotokawa and Ngatamariki geothermal fields

<p>In this thesis, we construct a four-year (2012–2015) catalog of microearthquakes for the Ngatamariki and Rotokawa geothermal fields in the Taupō Volcanic Zone of New Zealand, and use these data to improve the knowledge of reservoir behavior. These microearthquakes occur frequently, often every few seconds, and therefore provide a tool that we use to assess reservoir properties with dense spatial and temporal resolution as well as to illuminate the underlying processes of seismogenesis. Using a matched-filter detection technique we detect and precisely relocate nearly 9000 events, from which we calculate 982 focal mechanisms. At Ngatamariki, these results constitute the first detailed analysis of seismicity at a newly-developed resource. It has been commonly assumed that induced shear on fractures increases reservoir permeability by offsetting asperities on either fracture wall, thereby propping the fracture open. During stimulation treatments of two boreholes (NM08 and NM09), borehole permeability experiences logarithmic growth. At NM08, this growth occurs for eight days in the absence of seismicity, while at NM09 only nine microearthquakes are observed during the one-month treatment. This suggests that hydro-shear, the process of inducing seismicity through increased pore pressure at critically-stressed fractures, is not the dominant mechanism of permeability increase at many geothermal wells. Instead, aseismic processes, likely thermal and overpressure induced fracture opening, dominate well stimulation in high-temperature geothermal settings. At Rotokawa, the earthquake frequency-magnitude distribution (b-value) is positively correlated with both proximity to major injection wells and depth. In an inferred pressure compartment near injection well RK23, b is ~1.18, but is <1.0 elsewhere, suggesting a connection between increased pore-fluid pressure and small-magnitude events. In addition, throughout the reservoir b increases from a value of ~1.0 at injection depth to almost 1.5 two kilometers below the reservoir, consistent with observations at volcanic areas elsewhere, but opposing the conventional wisdom that b-value is inversely proportional to differential stress. Finally, the 982 focal mechanism observations that we invert for stress show a normal faulting regime throughout both reservoirs. At Rotokawa, a lowering stress ratio, v, after reintroduction of injection well RK23 (v drops from 0.9 to 0.2 over six months) indicates that anisotropic reservoir cooling affects the reservoir stress state through a process of preferential stress reduction.</p>

Characterizing microseismicity at the Rotokawa and Ngatamariki geothermal fields

<p>In this thesis, we construct a four-year (2012–2015) catalog of microearthquakes for the Ngatamariki and Rotokawa geothermal fields in the Taupō Volcanic Zone of New Zealand, and use these data to improve the knowledge of reservoir behavior. These microearthquakes occur frequently, often every few seconds, and therefore provide a tool that we use to assess reservoir properties with dense spatial and temporal resolution as well as to illuminate the underlying processes of seismogenesis. Using a matched-filter detection technique we detect and precisely relocate nearly 9000 events, from which we calculate 982 focal mechanisms. At Ngatamariki, these results constitute the first detailed analysis of seismicity at a newly-developed resource. It has been commonly assumed that induced shear on fractures increases reservoir permeability by offsetting asperities on either fracture wall, thereby propping the fracture open. During stimulation treatments of two boreholes (NM08 and NM09), borehole permeability experiences logarithmic growth. At NM08, this growth occurs for eight days in the absence of seismicity, while at NM09 only nine microearthquakes are observed during the one-month treatment. This suggests that hydro-shear, the process of inducing seismicity through increased pore pressure at critically-stressed fractures, is not the dominant mechanism of permeability increase at many geothermal wells. Instead, aseismic processes, likely thermal and overpressure induced fracture opening, dominate well stimulation in high-temperature geothermal settings. At Rotokawa, the earthquake frequency-magnitude distribution (b-value) is positively correlated with both proximity to major injection wells and depth. In an inferred pressure compartment near injection well RK23, b is ~1.18, but is <1.0 elsewhere, suggesting a connection between increased pore-fluid pressure and small-magnitude events. In addition, throughout the reservoir b increases from a value of ~1.0 at injection depth to almost 1.5 two kilometers below the reservoir, consistent with observations at volcanic areas elsewhere, but opposing the conventional wisdom that b-value is inversely proportional to differential stress. Finally, the 982 focal mechanism observations that we invert for stress show a normal faulting regime throughout both reservoirs. At Rotokawa, a lowering stress ratio, v, after reintroduction of injection well RK23 (v drops from 0.9 to 0.2 over six months) indicates that anisotropic reservoir cooling affects the reservoir stress state through a process of preferential stress reduction.</p>

Determining seismic shear-velocity from ambient noise sources at regional and local scales

<p>We present three projects that use different bandwidths of the ambient noise spectrum to solve geophysical problems. Specifically, we use signals within the noise field to determine surface and shear wave velocities, image the shallow and deep crust, and monitor time-dependent deformation resulting from geothermal fluid injection and extraction. Harrat Al-Madinah, a Cenozoic bimodal alkaline volcanic field in west-central Saudi Arabia, is imaged using shear-velocities obtained from natural ambient seismic noise. To our knowledge, this project is the first analysis of Saudi Arabia structure using ambient noise methods. Surface wave arrivals are extracted from a year's worth of station-pair cross-correlations, which are approximations of the empirical Green's function of the interstation path. We determine group and phase velocity surface wave dispersion maps with a 0.1 decimal degree resolution and resolve a zone of slow surface wave velocity south-east of the city of Medina, which is spatially correlated with the most recent historical eruption (the 1256 CE Medina eruption). Dispersion curves are calculated at each grid-point of the surface-wave velocity maps and inverted to obtain measurements of shear-velocity with depth. The 1D velocity models are then used to produce average shear-velocity models for the volcanic field. A shear-velocity increase ranging from 0.5 to 1.0 km/s, suggesting a layer interface, is detected at approximately 20 km depth and compared to P-wave measurement from a previous refraction study. We compute cross-section profiles by interpolating the inversions into a pseudo-3D model and resolve a zone of slow shear-velocity below the 1256 CE eruption location. These areas are also spatially correlated with low values of Bouguer gravity. We hypothesize that the low shear-velocity and gravity measurements are caused by fluids and fractures created from prior volcanic eruptions. We use the coda of cross-correlations extracted from ambient noise to determine shear-velocity changes at Rotokawa and Ngatamariki, two electricity producing geothermal fields located in the North Island of New Zealand. Stacks of cross correlations between stations prior to the onset of production are compared to cross correlations of moving stacks in time periods of well stimulation and the onset of electricity production using the Moving Window Cross Spectral technique. An increase between 0.05% to 0.1% of shear-velocity is detected at Rotokawa coinciding with an increase of injection. The shear-velocity subsequently decreases by approximately 0.1% when the rate of production surpasses the rate of injection. A similar amplitude shear-velocity increase is detected at Ngatamariki during the beginning of injection. After the initial increase, the shear-velocity at Ngatamariki fluctuates in response to differences in injection and production rates. A straight-ray pseudo-tomography analysis is conducted at the geothermal fields, which reveals that localized positive velocity changes are co-located with injection wells. Lastly, we use ambient noise and active sources at the Ngatamariki geothermal field to determine the structure of the top 200 meters using the Refraction Microtremor technique. We deployed a linear 72-channel array of vertical geophones with ten meter spacing at two locations of the geothermal field and determine average 1D and 2D shear-velocity profiles. We were able to image depths between 57 to 93 meters for 2D profiles and up to 165 meters for 1D profiles. A shear-velocity anomaly was detected across one of the lines that coincided with the inferred location of a fault determined from nearby well logs. This suggests that the method can be used to cheaply and quickly constrain near-surface geology at geothermal fields, where ambient noise is abundant and typical reflection and refraction surveys require large inputs of energy and are hindered by attenuation and scattering in near-surface layers.</p>

Determining seismic shear-velocity from ambient noise sources at regional and local scales

<p>We present three projects that use different bandwidths of the ambient noise spectrum to solve geophysical problems. Specifically, we use signals within the noise field to determine surface and shear wave velocities, image the shallow and deep crust, and monitor time-dependent deformation resulting from geothermal fluid injection and extraction. Harrat Al-Madinah, a Cenozoic bimodal alkaline volcanic field in west-central Saudi Arabia, is imaged using shear-velocities obtained from natural ambient seismic noise. To our knowledge, this project is the first analysis of Saudi Arabia structure using ambient noise methods. Surface wave arrivals are extracted from a year's worth of station-pair cross-correlations, which are approximations of the empirical Green's function of the interstation path. We determine group and phase velocity surface wave dispersion maps with a 0.1 decimal degree resolution and resolve a zone of slow surface wave velocity south-east of the city of Medina, which is spatially correlated with the most recent historical eruption (the 1256 CE Medina eruption). Dispersion curves are calculated at each grid-point of the surface-wave velocity maps and inverted to obtain measurements of shear-velocity with depth. The 1D velocity models are then used to produce average shear-velocity models for the volcanic field. A shear-velocity increase ranging from 0.5 to 1.0 km/s, suggesting a layer interface, is detected at approximately 20 km depth and compared to P-wave measurement from a previous refraction study. We compute cross-section profiles by interpolating the inversions into a pseudo-3D model and resolve a zone of slow shear-velocity below the 1256 CE eruption location. These areas are also spatially correlated with low values of Bouguer gravity. We hypothesize that the low shear-velocity and gravity measurements are caused by fluids and fractures created from prior volcanic eruptions. We use the coda of cross-correlations extracted from ambient noise to determine shear-velocity changes at Rotokawa and Ngatamariki, two electricity producing geothermal fields located in the North Island of New Zealand. Stacks of cross correlations between stations prior to the onset of production are compared to cross correlations of moving stacks in time periods of well stimulation and the onset of electricity production using the Moving Window Cross Spectral technique. An increase between 0.05% to 0.1% of shear-velocity is detected at Rotokawa coinciding with an increase of injection. The shear-velocity subsequently decreases by approximately 0.1% when the rate of production surpasses the rate of injection. A similar amplitude shear-velocity increase is detected at Ngatamariki during the beginning of injection. After the initial increase, the shear-velocity at Ngatamariki fluctuates in response to differences in injection and production rates. A straight-ray pseudo-tomography analysis is conducted at the geothermal fields, which reveals that localized positive velocity changes are co-located with injection wells. Lastly, we use ambient noise and active sources at the Ngatamariki geothermal field to determine the structure of the top 200 meters using the Refraction Microtremor technique. We deployed a linear 72-channel array of vertical geophones with ten meter spacing at two locations of the geothermal field and determine average 1D and 2D shear-velocity profiles. We were able to image depths between 57 to 93 meters for 2D profiles and up to 165 meters for 1D profiles. A shear-velocity anomaly was detected across one of the lines that coincided with the inferred location of a fault determined from nearby well logs. This suggests that the method can be used to cheaply and quickly constrain near-surface geology at geothermal fields, where ambient noise is abundant and typical reflection and refraction surveys require large inputs of energy and are hindered by attenuation and scattering in near-surface layers.</p>

Development of mathematical models for calculating parameters during transportation of a steam-water mixture

В практике освоения геотермальных месторождений возникает необходимость нахождение взаимосвязи между параметрами потока. При транспортировке пароводяной смеси требуется определение условий, обеспечивающих стабильную работу трубопровода. На отечественных объектах геотермальной энергетики гидравлический расчет трубопроводов пароводяной смеси проводился с использованием компьютерной программы MODEL, разработанной авторами данной работы. Новые вызовы и появление теории устойчивости привели к созданию математической модели SWIP (Steam-Water Inclining Pipeline). В данной статье представлены этапы разработки авторами новой модели, отвечающей современным требованиям In the practice of developing geothermal fields, there is a need to find the dependence between the flow parameters. When transporting a steam-water mixture, it is necessary to determine the conditions that ensure the stability in the pipeline. At domestic geothermal energy, the hydraulic calculation of steam-water mixture pipelines was carried out using the computer program MODEL created by the authors of this work. New challenges and the emergence of the theory of stability led to the creation of a mathematical model SWIP (Stream-Water Inclining Pipeline). This article presents the stages of creating a new model by the authors that meets modern requirements.

Comparative analysis of computer programs for hydraulic calculation of steam-water mixture in pipelines

In the practice of developing geothermal fields, it is necessary to forecast the flow parameters of the well and the pressure drop in the pipeline and need to decide on the optimal design for stable operation of the pipeline. At domestic facilities of geothermal energy, the hydraulic calculation of the steam-water mixture pipelines was carried out using the computer program MODEL developed by the authors of this work. New challenges and the emergence of stability theory led to the creation of the mathematical model SWIP (Steam-Water Inclining Pipeline). This paper presents a comparative analysis of two options for implementing the new model.