1,721,037 research outputs found
Fracture clustering effect on amplitude variation with offset and azimuth analyses
Traditional amplitude variation with offset and azimuth (AVOAz) analysis for fracture characterization extracts fracture properties through analysis of reflection AVOAz to determine anisotropic parameters (e.g., Thomsen’s parameters) that are then related to fracture properties. The validity of this method relies on the basic assumption that a fractured unit can be viewed as an equivalent anisotropic medium. As a rule of thumb, this assumption is taken to be valid when the fracture spacing is less than λ/10. Under the effective medium assumption, diffractions from individual fractures destructively interfere and only specular reflections from boundaries of a fractured layer can be observed in seismic data. The effective medium theory has been widely used in fracture characterization, and its applicability has been validated through many field applications. However, through numerical simulations, we find that diffractions from fracture clusters can significantly distort the AVOAz signatures when a fracture system has irregular spacing even though the average fracture spacing is much smaller than a wavelength (e.g., ≪λ/10). Contamination by diffractions from irregularly spaced fractures on reflections can substantially bias the fracture properties estimated from AVOAz analysis and may possibly lead to incorrect estimates of fracture properties. Additionally, through Monte Carlo simulations, we find that fracture spacing uncertainty inverted from amplitude variation with offset (AVO) analysis can be up to 10%–20% when fractures are not uniformly distributed, which should be the realistic state of fractures present in the earth. Also, AVOAz and AVO analysis gives more reliable estimates of fracture properties when reflections at the top of the fractured layer are used compared with those from the bottom of the layer
Investigation of borehole cross-dipole flexural dispersion crossover through numerical modeling
Crossover of the dispersion of flexural waves recorded in borehole cross-dipole measurements is interpreted as an indicator of stress-induced anisotropy around a circular borehole in formations that are isotropic in the absence of stresses. We have investigated different factors that influence flexural wave dispersion. Through numerical modeling, we determined that for a circular borehole surrounded by an isotropic formation that is subjected to an anisotropic stress field, the dipole flexural dispersion crossover is detectable only when the formation is very compliant. This might happen only in the shallow subsurface or in zones having high pore pressure. However, we found that dipole dispersion crossover can also result from the combined effect of formation intrinsic anisotropy and borehole elongation. We found that a small elongation on the wellbore and very weak intrinsic anisotropy can result in a resolvable crossover in flexural dispersion that might be erroneously interpreted as borehole stress-induced anisotropy. A thorough and correct interpretation of flexural dispersion crossover thus has to take into account the effects of stress-induced and intrinsic anisotropy and borehole cross-sectional geometry
Simulation of the effect of stress-induced anisotropy on borehole compressional wave propagation
Formation elastic properties near a borehole may be altered from their original state due to the stress concentration around the borehole. This can lead to an incorrect estimation of formation elastic properties measured from sonic logs. Previous work has focused on estimating the elastic properties of the formation surrounding a borehole under anisotropic stress loading. We studied the effect of borehole stress concentration on sonic logging in a moderately consolidated Berea sandstone using a two-step approach. First, we used an iterative approach, which combines a rock-physics model and a finite-element method, to calculate the stress-dependent elastic properties of the rock around a borehole subjected to an anisotropic stress loading. Second, we used the anisotropic elastic model obtained from the first step and a finite-difference method to simulate the acoustic response of the borehole. Although we neglected the effects of rock failure and stress-induced crack opening, our modeling results provided important insights into the characteristics of borehole P-wave propagation when anisotropic in situ stresses are present. Our simulation results were consistent with the published laboratory measurements, which indicate that azimuthal variation of the P-wave velocity around a borehole subjected to uniaxial loading is not a simple cosine function. However, on field scale, the azimuthal variation in P-wave velocity might not be apparent at conventional logging frequencies. We found that the low-velocity region along the wellbore acts as an acoustic focusing zone that substantially enhances the P-wave amplitude, whereas the high-velocity region caused by the stress concentration near the borehole results in a significantly reduced P-wave amplitude. This results in strong azimuthal variation of P-wave amplitude, which may be used to infer the in situ stress state
3D weak-dispersion reverse time migration using a stereo-modeling operator
Reliable 3D imaging is a required tool for developing models of complex geologic structures. Reverse time migration (RTM), as the most powerful depth imaging method, has become the preferred imaging tool because of its ability to handle complex velocity models including steeply dipping interfaces and large velocity contrasts. Finite-difference methods are among the most popular numerical approaches used for RTM. However, these methods often encounter a serious issue of numerical dispersion, which is typically suppressed by reducing the grid interval of the propagation model, resulting in large computation and memory requirements. In addition, even with small grid spacing, numerical anisotropy may degrade images or, worse, provide images that appear to be focused but position events incorrectly. Recently, stereo-operators have been developed to approximate the partial differential operator in space. These operators have been used to develop several weak-dispersion and efficient stereo-modeling methods that have been found to be superior to conventional algorithms in suppressing numerical dispersion and numerical anisotropy. We generalized one stereo-modeling method, fourth-order nearly analytic central difference (NACD), from 2D to 3D and applied it to 3D RTM. The RTM results for the 3D SEG/EAGE phase A classic data set 1 and the SEG Advanced Modeling project model demonstrated that, even when using a large grid size, the NACD method can handle very complex velocity models and produced better images than can be obtained using the fourth-order and eighth-order Lax-Wendroff correction (LWC) schemes. We also applied 3D NACD and fourth-order LWC to a field data set and illustrated significant improvements in terms of structure imaging, horizon/layer continuity and positioning. We also investigated numerical dispersion and found that not only does the NACD method have superior dispersion characteristics but also that the angular variation of dispersion is significantly less than for LWC.
Read More: http://library.seg.org/doi/abs/10.1190/geo2013-0472.1National Natural Science Foundation (China) (Grant 41230210)Massachusetts Institute of Technology. Earth Resources Laboratory (Founding Members Consortium
Investigation of collar properties on data-acquisition scheme for acoustic logging-while-drilling
We have used the wavenumber integration, velocity-time semblance, and dispersion methods to investigate the influence of collar properties including velocities, density, and attenuation on acoustic logging-while-drilling wavefields. We have found that when the velocities of the collar wave and the P-wave of the formation are similar, they interfere. However, the interference disappears when the velocity difference increases. Having a collar with large velocities (especial large shear velocity) and density makes the direct P-velocity determination possible in a fast formation even without isolators. For a slow formation, the interference of the collar flexural wave with the formation flexural and leaky P-waves is slight for a dipole tool when collar velocities are large. For this case, the S velocity can be determined by the flexural formation wave at low frequency (approximately 2 kHz). Based on these observations, we propose that the measurement of the P- and S-velocities can be easier if the collar is made of an advanced composite material that has high compressional and shear velocities as well as density. This is a direct and easy change to implement and a new idea for an acoustic logging-while-drilling tool design
Investigation of the high-frequency wavefield of an off-center monopole acoustic logging-while-drilling tool
During logging-while-drilling (LWD) operations, complex drill string movements and the weight of the drill pipe often lead to a measurement tool that is not centralized. Therefore, studies of the response of an off-center acoustic LWD tool are essential to facilitate better interpretation of measurements made in an actual drilling environment. Such studies will be helpful for tool design and data processing. We used a finite-difference method to simulate the response of a noncentralized monopole acoustic LWD tool at high frequency (10 kHz). We analyzed the effects on the waveforms for receivers at different azimuths caused by an off-center tool with differing amounts of offset. We used velocity-time semblance and dispersion analysis methods to help us to understand the modes in the waveforms at different azimuth receivers for different tool offsets. We have found that the waveforms in the direction of the tool offset, that is, where the fluid column is smallest, were affected the most. Waveforms in the orthogonal direction were less affected by tool offset. Collar flexural and collar quadrupole modes appear when the tool is off center. In addition, the formation flexural and quadrupole modes contaminate the Stoneley wave. Waveforms in a fast formation are more strongly affected by the offset of the tool than those in a slow formation. In a fast formation, the new collar modes make it difficult to determine the P-wave velocity in the direction of tool offset whereas it is easier in the orthogonal direction. However, P-waves are less contaminated by new modes in a slow formation. Due to the significant changes in waveforms with azimuth when the tool is off center, the simple addition of all waveforms from an azimuthal distribution of receivers will not result in a clean waveform that is sensitive to only the surrounding formation.National Natural Science Foundation (China) (41174118)National Natural Science Foundation (China) (41404100)S&T Special Project (2011ZX05020-009)China Postdoctoral Science Foundation (2013M530106)International Postdoctoral Exchange Fellowship Progra
Converted-wave seismic imaging: Amplitude-balancing source-independent imaging conditions
We have developed crosscorrelational and deconvolutional forms of a source-independent converted-wave imaging condition (SICW-IC) and show the relationship between them using a concept of conversion ratio coefficient, a concept that we developed through reflection, transmission, and conversion coefficients. We applied the SICW-ICs to a two half-space model and the synthetic Marmousi I and II models and show the sensitivity of the SICW-ICs to incorrect wave speed models. We also compare the SICW-ICs and source-dependent elastic reverse time migration. The results of SICW-ICs highlight the improvements in spatial resolution and amplitude balancing with the deconvolutional forms. This is an attractive alternative to active and passive source elastic imaging
Estimation of elastic constants for HTI media using Gauss-Newton and full-Newton multiparameter full-waveform inversion
In seismic full-waveform inversion (FWI), subsurface parameters are estimated by iteratively minimizing the difference between the modeled and the observed data. We have considered the problem of estimating the elastic constants of a fractured medium using multiparameter FWI and modeling naturally fractured reservoirs as equivalent anisotropic media. Multiparameter FWI, although promising, remains exposed to a range of challenges, one being the parameter crosstalk problem resulting from the overlap of Fréchet derivative wavefields. Parameter crosstalk is strongly influenced by the form of the scattering pattern for each parameter. We have derived 3D radiation patterns associated with scattering from a range of elastic constants in general anisotropic media. Then, we developed scattering patterns specific to a horizontal transverse isotropic (HTI) medium to draw conclusions about parameter crosstalk in FWI. Bare gradients exhibit crosstalk, as well as artifacts caused by doubly scattered energy in the data residuals. The role of the multiparameter Gauss-Newton (GN) Hessian in suppressing parameter crosstalk is revealed. We have found that the second-order term in the multiparameter Hessian, which is associated with multiparameter second-order scattering effects, can be constructed with the adjoint-state technique. We have examined the analytic scattering patterns for HTI media with a 2D numerical example. We have examined the roles played by the first- and second-order terms in multiparameter Hessian to suppress parameter crosstalk and second-order scattering artifacts numerically. We have also compared the multiparameter GN and full-Newton methods as methods for determining the elastic constants in HTI media with a two-block-layer model.Consortium for Research in Elastic Wave Exploration SeismologyNatural Sciences and Engineering Research Council of Canada (Grant CRDPJ 461179-13
Going Beyond Counting First Authors in Author Co-citation Analysis
The present study examines one of the fundamental aspects of author co-citation analysis (ACA) - the way co-citation
counts are defined. Co-citation counting provides the data on which all subsequent statistical analyses and mappings
are based, and we compare ACA results based on two different types of co-citation counting - the traditional type that
only counts the first one among a cited work's authors on the one hand and a non-traditional type that takes into
account the first 5 authors of a cited work on the other hand. Results indicate that the picture produced through this non-traditional author co-citation counting contains more coherent author groups and is therefore considerably clearer. However, this picture represents fewer specialties in the research field being studied than that produced through the traditional first-author co-citation counting when the same number of top-ranked authors is selected and analyzed. Reasons for these effects are discussed
- …
