Rock strength and elastic properties are primary rock parameters. Many rocks are inherently anisotropic. Knowledge of strength anisotropy and elastic anisotropy has great practical application in stimulation design. Twelve oriented pyrophyllite core plugs were tested to analyze the effects of anisotropy. The orientation angles between core axes and the anisotropy axis range from 0° to 90° with 15° increments. Ultrasonic pulse transmission technique is employed to measure P-wave velocity and S-wave velocities. Uniaxial compression tests give rock strength and static modulus from stress-strain data. Static moduli are compared with dynamic moduli calculated from wave velocities.
The results indicate that compressive strength, wave velocities and elastic moduli are affected significantly by anisotropy. Pyrophylllite has an asymmetrical U-shaped strength anisotropy behavior. The greatest strength is measured parallel to fabric, 108.7 MPa. The lowest strength is 88.3 MPa found for samples cored at 45° orientation to the fabric. The strength anisotropy behavior of pyrophyllite can be predicted by knowing the strength values at 0°, 90° and the orientation angle of the minimum strength (approximately 45° in this study). P-wave and S-wave velocities, dynamic and static Young’s moduli are strongly anisotropic. The variation of measured P-wave and S-wave phase velocities agree with the predicted model (Thomsen 1986). Dynamic Young’s modulus varies from 38.6 GPa to 53.6 GPa while the static modulus varies from 13.9 GPa to 32.5 GPa. Dynamic modulus is generally higher than the static for all the orientations. The ratio of dynamic to static moduli ranges from 1.65 to 2.83. A cross-plot of strength anisotropy ratio versus static modulus anisotropy ratio for pyrophyllite and other rock types reported in the literature is presented. No obvious correlation is found between these two parameters, which may be due to the limited number of published data..
Today, the problem of surface-related multiples, especially in shallow water, is not fully solved. Although surface-related multiple elimination (SRME) method has proved to be successful on a large number of data cases, the involved adaptive subtraction acts as a weak link in this methodology, where primaries can be distorted due to their interference with multiples. Therefore, recently, SRME has been redefined as a large-scale inversion process, called estimation of primaries by sparse inversion (EPSI). In this process the multi-dimensional primary impulse responses are considered as the unknowns in a large-scale inversion process. By parameterizing these impulse responses as spikes in the space-time domain, and using a sparsity constraint in the update step, the algorithm looks for those primaries that, together with their associated multiples, explain the total input data. As the objective function in this minimization process truly goes to zero, the tendency for distorting primaries is greatly reduced. An additional advantage is that imperfections in the data can be included in the forward model and resolved simultaneously, such as the missing near offsets. In this paper it is demonstrated that the ghost effect can also be included in the EPSI formulation after which a ghost-free primary estimate can be obtained, even in the case the ghost notch is within the desired spectrum.
Full wavefield migration (FWM) is an inversion-based imaging algorithm that utilizes the complete reflection measurements: primaries as well as all multiples, both surface and internal. Using multiples in the imaging can extend the illumination of the subsurface. In this paper we concentrate the study on the internal scattering that can be helpful in imaging structures from below that are otherwise difficult to image by primaries, i.e. an undershooting setting. This can be fruitful in the case of the obstacles, like oil-production facilities, or in the case of poor illumination, like sub-salt imaging or near-surface complexities. We demonstrate such approach on two synthetic examples.
Tomographic full waveform inversion (TFWI) provides a framework to invert the seismic data that is immune to cycle-skipping problems. This is achieved by extending the wave equation and adding an offset axis to the velocity model. However, this extension makes the propagation considerably more expensive because each multiplication by velocity becomes a convolution. We provide an alternative formulation which computes the backscattering and the forward scattering components of the gradient separately. To maintain high resolution results of TFWI, the two components of the gradient are first mixed and then separated based on a Fourier domain scale separation. This formulation is based on the Born approximation where the medium parameters are broken into a long-wavelength and short-wavelength components. The inversion setup includes two steps that maintain the high resolution results of TFWI. First, the linearized residual are updated in a nested inversion scheme. This step corrects for the underlying assumption that the data contain primaries only without multiples. Second, the two components of the gradient are first mixed and then separated based on a Fourier domain scale separation to allow for a fully simultaneous inversion of model scales. After deriving the equations, we test the theory with two synthetic examples. The results of both the Marmousi and BP models show that convergence is possible even with large errors in the initial model that would have prevented convergence to conventional FWI.
Han, Jianfa (Tarim Oilfield Company, CNPC) | Sun, Sam Zandong (China University of Petroleum) | Zhang, Xinchao (China University of Petroleum) | Zhang, Yuanyin (China University of Petroleum) | Chen, Jun (Tarim Oilfield Company, CNPC) | Pan, Yangyong (Tarim Oilfield Company, CNPC) | Liu, Xin (Tarim Oilfield Company, CNPC) | Zhao, Haitao (Tarim Oilfield Company, CNPC)
In order to improve the SNR of seismic data, the multi-coverage technique is developed and it is a common practice to stack all trace information of a CRP as zero-incident data. In this process, although the SNR is greatly enhanced, yet the resolution is inevitably lost since the AVO effects are smeared. Moreover, the traditional post-stack seismic data is not the true pure P-wave data. By contrast, the pre-stack inversion can obtain the true P-wave data, and the resolution is much higher than the post-stack data since this method removes the effects of S-wave information. The study in this paper finds that feasible and reliable evaluations on the reservoir quality could be realized, based on which we can predict the high quality reservoirs. What's more, the mud-filled area at the top of Yingshan Formation is successfully identified, and could provide a unique and new insight for studying the development rules of carbonate reservoir.
This study assessed possible application of the friable sand model for shallow mechanically compacted overconsolidated sands from experientially derived velocity and porosity relations of six brine saturated sand aggregates and three published sand compaction datasets. The results showed that the friable sand model established with the assumption of the sediments normally consolidated can also be used for overconsolidated sands caused by stress released as long as the pore pressure is hydrostatic during unloading. We also found that the friable sand model used for overconsolidated sands not only describes change in depositional sorting but also variation in preconsolidation stress related to the amount of uplift. The study outcomes expand the rock physics diagnostic approach to predict seismic properties of shallow overconsolidated sands that have undergone complex burial history (i.e., loading-unloading-reloading) in uplifted basins like the Barents Sea.
The experience of integration of a broadband seismic dataset on CLOV project is presented. The fields to be developed are based in turbiditic reservoirs of Oligocene (Cravo, Lirio and Orquidea) and Miocene ages (Orquidea-Violeta) trapped within faulted structures in the passive margin transitional domain.
The survey acquired early 2012 with Geostreamer1 technology, mainly as 4D baseline, is replacing a reference 3DHR dataset in the geoscience field evaluation process. Hence in the optimization of the development wells and operational follow up, drawing benefits from the quality improvements associated with bandwidth broadening. Further S/N enhancement around the challenging areas results from more coherent low frequency content and more elaborated processing sequence.
As it proves to be valuable in terms of 3D seismic interpretation, too, the CLOV broadband dataset appears reliable as baseline dataset for the planed 4D reservoir surveillance seismic.
We develop a gradient computation for recovering anisotropic parameters associated with tilted-transverse isotropy (TTI) in full-waveform inversion (FWI). Our treatment is based on adjoint state theory and pseudo-analytic TTI wave propagation. The method provides reasonable estimates for anisotropic parameter updates and we observe that the migrated images for both the synthetic and real data show improvements with the inverted anisotropic parameters. However, due to the huge null space associated with the inversion we occasionally obtain parameters that are not in line with our understanding of the geology. To overcome this, we introduce a priori information within our inversion via dip annihilation filters. The introduction of this model styling goal within our inversion not only makes the results aesthetically pleasing but also in accordance with our physical intuition about the model. We illustrate our method with the help of a synthetic 3D example.
A series of geophysical surveys using ground-penetrating radar (GPR) and electromagnetic induction (EM) were conducted at Padre Island National Seashore (PAIS), Texas to characterize barrier island transgression at a range of spatial and temporal scales. The first study (June, 2010 - January, 2011) used GPR and vibracores to investigate variations within the internal structure of various-sized foredunes along a 2.5 km length of the non-driving (restricted vehicular access) section of PAIS. Three representative foredunes with different heights were selected within the storm impact regime proposed by Sallenger (2000). A small (1.79 m), intermediate (2.69 m) and large dune (3.77 m) represent inundation, overwash and collision regimes, respectively.
The second preliminary study (May, 2013) used the EM induction method to measure spatial variations in apparent conductivity across the intermediate dune site surveyed by GPR in 2010-2011. A 100 m long EM transect was conducted in close proximity to the 2010-2011 GPR survey area. Additionally, a new GPR survey was conducted along the same EM profile for comparison with the EM data as well as the GPR profile taken in 2010-2011. Currently, we plan to conduct several extensive GPR and EM surveys at regularly-spaced intervals (ㅇ km apart) within the driving (vehicular accessible) section of PAIS and across the Laguna Madre wind-tidal flat system.
We describe a modification of the design criteria normally adopted for Offset Vector Tiles, so that they can be used more effectively with converted-wave data. Our modification is based on the asymptotic conversion point correction and is a function of the Vp/Vs ratio.
There are a number of potential applications of this approach in converted-wave processing. Here, we demonstrate the improvement that results in converted-wave prestack time migration of single offset vector tiles using this PS Offset Vector Tile design.