A simple inversion algorithm to estimate a linearly increasing velocity model for microseismic monitoring

2018 ◽  
Vol 49 (5) ◽  
pp. 647-654
Author(s):  
Woochang Choi ◽  
Wonsik Kim ◽  
Sukjoon Pyun
2018 ◽  
Vol 6 (3) ◽  
pp. SH39-SH48 ◽  
Author(s):  
Wojciech Gajek ◽  
Jacek Trojanowski ◽  
Michał Malinowski ◽  
Marek Jarosiński ◽  
Marko Riedel

A precise velocity model is necessary to obtain reliable locations of microseismic events, which provide information about the effectiveness of the hydraulic stimulation. Seismic anisotropy plays an important role in microseismic event location by imposing the dependency between wave velocities and its propagation direction. Building an anisotropic velocity model that accounts for that effect allows for more accurate location of microseismic events. We have used downhole microseismic records from a pilot hydraulic fracturing experiment in Lower-Paleozoic shale gas play in the Baltic Basin, Northern Poland, to obtain accurate microseismic events locations. We have developed a workflow for a vertical transverse isotropy velocity model construction when facing a challenging absence of horizontally polarized S-waves in perforation shot data, which carry information about Thomsen’s [Formula: see text] parameter and provide valuable constraints for locating microseismic events. We extract effective [Formula: see text], [Formula: see text] and [Formula: see text], [Formula: see text] for each layer from the P- and SV-wave arrivals of perforation shots, whereas the unresolved [Formula: see text] is retrieved afterward from the SH-SV-wave delay time of selected microseismic events. An inverted velocity model provides more reliable location of microseismic events, which then becomes an essential input for evaluating the hydraulic stimulation job effectiveness in the geomechanical context. We evaluate the influence of the preexisting fracture sets and obliquity between the borehole trajectory and principal horizontal stress direction on the hydraulic treatment performance. The fracturing fluid migrates to previously fractured zones, while the growth of the microseismic volume in consecutive stages is caused by increased penetration of the above-lying lithologic formations.


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCB47-WCB55 ◽  
Author(s):  
Donghong Pei ◽  
John A. Quirein ◽  
Bruce E. Cornish ◽  
Dan Quinn ◽  
Norman R. Warpinski

To accurately locate microearthquakes that are genetically related to hydraulic fracture stimulation, a thorough knowledge of the velocity structure between monitoring and fracturing treatment wells is essential. Very fast simulated annealing (VFSA) is implemented to invert for a flat-layered velocity model between wells using perforation or string-shot data. A two-point ray-tracing method is used to find the ray parameter [Formula: see text] for a ray traveling from a source to a receiver. The original traveltime-calculation formula is modified to account for the borehole source-receiver geometry. VFSA is used as a tool to optimize P- and S-wave velocities simultaneously. Unlike previous applications of VFSA, two improvements result from a new study: (1) both P- and S-wave arrival-time misfits are considered in a joint-objective function, and (2) P- and S-wave velocities are perturbed simultaneously during annealing. The inverted velocities follow the true values closely with a very small root-mean-square error, indicating the inverted model is close to the global minimum solution whose rms error should be zero for synthetic examples. Data noise contaminates inverted models, but not substantially in synthetic test results. A comparison of models inverted using VFSA and Occam’s inversion technique indicates that inverted models using VFSA are superior to those using Occam’s method in terms of velocity accuracy.


2016 ◽  
Vol 135 ◽  
pp. 317-327 ◽  
Author(s):  
Haiyu Jiang ◽  
Zhongren Wang ◽  
Xiaoxian Zeng ◽  
Hao Lü ◽  
Xiaohua Zhou ◽  
...  

Geosciences ◽  
2020 ◽  
Vol 10 (9) ◽  
pp. 327
Author(s):  
Zhiyong Zhang ◽  
Diego Arosio ◽  
Azadeh Hojat ◽  
Luigi Zanzi

To monitor the stability of a mountain slope in northern Italy, microseismic monitoring technique has been used since 2013. Locating microseismic events is a basic step of this technique. We performed a seismic tomographic survey on the mountain surface above the rock face to obtain a reliable velocity distribution in the rock mass for the localization procedure. Seismic travel-time inversion showed high heterogeneity of the rock mass with strong contrast in velocity distribution. Low velocities were found at shallow depth on the top of the rock cliff and intermediate velocities were observed in the most critical area of the rock face corresponding to a partially detached pillar. Using the 3D velocity model obtained from inversion, localization tests were performed based on the Equal Differential Time (EDT) localization method. The results showed hypocenter misfits to be around 15 m for the five geophones of the microseismic network and the error was significantly decreased compared to the results produced by a constant velocity model. Although the localization errors are relatively large, the accuracy is sufficient to distinguish microseismic events occurring in the most critical zone of the monitored rock mass from microseismic events generated far away. Thus, the 3D velocity model will be used in future studies to improve the classification of the recorded events.


Geophysics ◽  
2009 ◽  
Vol 74 (3) ◽  
pp. R15-R24 ◽  
Author(s):  
Taeyoung Ha ◽  
Wookeen Chung ◽  
Changsoo Shin

Waveform inversion faces difficulties when applied to real seismic data, including the existence of many kinds of noise. The [Formula: see text]-norm is more robust to noise with outliers than the least-squares method. Nevertheless, the least-squares method is preferred as an objective function in many algorithms because the gradient of the [Formula: see text]-norm has a singularity when the residual becomes zero. We propose a complex-valued Huber function for frequency-domain waveform inversion that combines the [Formula: see text]-norm (for small residuals) with the [Formula: see text]-norm (for large residuals). We also derive a discretized formula for the gradient of the Huber function. Through numerical tests on simple synthetic models and Marmousi data, we find the Huber function is more robust to outliers and coherent noise. We apply our waveform-inversion algorithm to field data taken from the continental shelf under the East Sea in Korea. In this setting, we obtain a velocity model whose synthetic shot profiles are similar to the real seismic data.


Geophysics ◽  
2020 ◽  
pp. 1-79
Author(s):  
Can Oren ◽  
Jeffrey Shragge

Accurately estimating event locations is of significant importance in microseismic investigations because this information greatly contributes to the overall success of hydraulic fracturing monitoring programs. Full-wavefield time-reverse imaging (TRI) using one or more wave-equation imaging conditions offers an effective methodology for locating surface-recorded microseismic events. To be most beneficial in microseismic monitoring programs, though, the TRI procedure requires using accurate subsurface models that account for elastic media effects. We develop a novel microseismic (extended) PS energy imaging condition that explicitly incorporates the stiffness tensor and exhibits heightened sensitivity to isotropic elastic model perturbations compared to existing imaging conditions. Numerical experiments demonstrate the sensitivity of microseismic TRI results to perturbations in P- and S-wave velocity models. Zero-lag and extended microseismic source images computed at selected subsurface locations yields useful information about 3D P- and S-wave velocity model accuracy. Thus, we assert that these image volumes potentially can serve as the input into microseismic elastic velocity model building algorithms.


Sign in / Sign up

Export Citation Format

Share Document