Identification of gas hydrate dissociation from wireline-log data in the Shenhu area, South China Sea

Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. B125-B134 ◽  
Author(s):  
Xiujuan Wang ◽  
Myung Lee ◽  
Shiguo Wu ◽  
Shengxiong Yang
Keyword(s):  
Wave Velocity ◽  
Gas Hydrate ◽  
Seismic Data ◽  
Pore Space ◽  
P Wave ◽  
Well Log ◽  
Hydrate Dissociation ◽  
Free Gas ◽  
P Wave Velocity

Wireline logs were acquired in eight wells during China’s first gas hydrate drilling expedition (GMGS-1) in April–June of 2007. Well logs obtained from site SH3 indicated gas hydrate was present in the depth range of 195–206 m below seafloor with a maximum pore-space gas hydrate saturation, calculated from pore water freshening, of about 26%. Assuming gas hydrate is uniformly distributed in the sediments, resistivity calculations using Archie’s equation yielded hydrate-saturation trends similar to those from chloride concentrations. However, the measured compressional (P-wave) velocities decreased sharply at the depth between 194 and 199 mbsf, dropping as low as [Formula: see text], indicating the presence of free gas in the pore space, possibly caused by the dissociation of gas hydrate during drilling. Because surface seismic data acquired prior to drilling were not influenced by the in situ gas hydrate dissociation, surface seismic data could be used to identify the cause of the low P-wave velocity observed in the well log. To determine whether the low well-log P-wave velocity was caused by in situ free gas or by gas hydrate dissociation, synthetic seismograms were generated using the measured well-log P-wave velocity along with velocities calculated assuming both gas hydrate and free gas in the pore space. Comparing the surface seismic data with various synthetic seismograms suggested that low P-wave velocities were likely caused by the dissociation of in situ gas hydrate during drilling.

Characterization of elastic properties of near-surface and subsurface deepwater hydrate-bearing sediments

Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. D169-D179 ◽  
Author(s):  
Zijian Zhang ◽  
De-hua Han ◽  
Daniel R. McConnell
Keyword(s):  
Wave Velocity ◽  
Elastic Constants ◽  
Gas Hydrate ◽  
Effective Medium Theory ◽  
Pore Space ◽  
Seismic Interpretation ◽  
P Wave ◽  
Near Surface ◽  
Free Gas ◽  
Hydrate Bearing Sediments

Hydrate-bearing sands and shallow nodular hydrate are potential energy resources and geohazards, and they both need to be better understood and identified. Therefore, it is useful to develop methodologies for modeling and simulating elastic constants of these hydrate-bearing sediments. A gas-hydrate rock-physics model based on the effective medium theory was successfully applied to dry rock, water-saturated rock, and hydrate-bearing rock. The model was used to investigate the seismic interpretation capability of hydrate-bearing sediments in the Gulf of Mexico by computing elastic constants, also known as seismic attributes, in terms of seismic interpretation, including the normal incident reflectivity (NI), Poisson’s ratio (PR), P-wave velocity ([Formula: see text]), S-wave velocity ([Formula: see text]), and density. The study of the model was concerned with the formation of gas hydrate, and, therefore, hydrate-bearing sediments were divided into hydrate-bearing sands, hydrate-bearing sands with free gas in the pore space, and shallow nodular hydrate. Although relations of hydrate saturation versus [Formula: see text] and [Formula: see text] are different between structures I and II gas hydrates, highly concentrated hydrate-bearing sands may be interpreted on poststack seismic amplitude sections because of the high NI present. The computations of elastic constant implied that hydrate-bearing sands with free gas could be detected with the crossplot of NI and PR from prestack amplitude analysis, and density may be a good hydrate indicator for shallow nodular hydrate, if it can be accurately estimated by seismic methods.

scholarly journals Acoustic Velocity Log Numerical Simulation and Saturation Estimation of Gas Hydrate Reservoir in Shenhu Area, South China Sea

2013 ◽  
Vol 2013 ◽  
pp. 1-13 ◽  
Author(s):  
Kun Xiao ◽  
Changchun Zou ◽  
Biao Xiang ◽  
Jieqiong Liu
Keyword(s):  
Numerical Simulation ◽  
South China Sea ◽  
Wave Velocity ◽  
South China ◽  
Gas Hydrate ◽  
P Wave ◽  
Shenhu Area ◽  
Free Gas ◽  
P Wave Velocity ◽  
Hydrate Saturation

Gas hydrate model and free gas model are established, and two-phase theory (TPT) for numerical simulation of elastic wave velocity is adopted to investigate the unconsolidated deep-water sedimentary strata in Shenhu area, South China Sea. The relationships between compression wave (P wave) velocity and gas hydrate saturation, free gas saturation, and sediment porosity at site SH2 are studied, respectively, and gas hydrate saturation of research area is estimated by gas hydrate model. In depth of 50 to 245 m below seafloor (mbsf), as sediment porosity decreases, P wave velocity increases gradually; as gas hydrate saturation increases, P wave velocity increases gradually; as free gas saturation increases, P wave velocity decreases. This rule is almost consistent with the previous research result. In depth of 195 to 220 mbsf, the actual measurement of P wave velocity increases significantly relative to the P wave velocity of saturated water modeling, and this layer is determined to be rich in gas hydrate. The average value of gas hydrate saturation estimated from the TPT model is 23.2%, and the maximum saturation is 31.5%, which is basically in accordance with simplified three-phase equation (STPE), effective medium theory (EMT), resistivity log (Rt), and chloride anomaly method.

Using Seismic Data to Predict Shale Pore Pressure and Overpressure Sweet Spots: A Case Study from the Lower Silurian Longmaxi Formation in Sichuan Basin, China

2021 ◽  
Author(s):  
Sheng Chen ◽  
Qingcai Zeng ◽  
Xiujiao Wang ◽  
Qing Yang ◽  
Chunmeng Dai ◽  
...  
Keyword(s):  
Prediction Model ◽  
Wave Velocity ◽  
Shale Gas ◽  
Seismic Data ◽  
P Wave ◽  
Formation Pressure ◽  
S Wave ◽  
P Wave Velocity ◽  
New Formation ◽  
Sweet Spots

Abstract Practices of marine shale gas exploration and development in south China have proved that formation overpressure is the main controlling factor of shale gas enrichment and an indicator of good preservation condition. Accurate prediction of formation pressure before drilling is necessary for drilling safety and important for sweet spots predicting and horizontal wells deploying. However, the existing prediction methods of formation pore pressures all have defects, the prediction accuracy unsatisfactory for shale gas development. By means of rock mechanics analysis and related formulas, we derived a formula for calculating formation pore pressures. Through regional rock physical analysis, we determined and optimized the relevant parameters in the formula, and established a new formation pressure prediction model considering P-wave velocity, S-wave velocity and density. Based on regional exploration wells and 3D seismic data, we carried out pre-stack seismic inversion to obtain high-precision P-wave velocity, S-wave velocity and density data volumes. We utilized the new formation pressure prediction model to predict the pressure and the spatial distribution of overpressure sweet spots. Then, we applied the measured pressure data of three new wells to verify the predicted formation pressure by seismic data. The result shows that the new method has a higher accuracy. This method is qualified for safe drilling and prediction of overpressure sweet spots for shale gas development, so it is worthy of promotion.

scholarly journals Simulation of the effect of stress-induced anisotropy on borehole compressional wave propagation

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. D205-D216 ◽  
Author(s):  
Xinding Fang ◽  
Michael C. Fehler ◽  
Arthur Cheng
Keyword(s):  
Wave Propagation ◽  
Stress Concentration ◽  
Elastic Properties ◽  
Wave Velocity ◽  
Wave Amplitude ◽  
P Wave ◽  
Azimuthal Variation ◽  
P Wave Velocity ◽  
Velocity Region

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.

Elastic reflectivity vectors and the geometry of intercept-gradient crossplots

2019 ◽  
Vol 38 (10) ◽  
pp. 762-769
Author(s):  
Patrick Connolly
Keyword(s):  
Elastic Property ◽  
Elastic Properties ◽  
Wave Velocity ◽  
Seismic Data ◽  
Measurement Errors ◽  
Rock Physics ◽  
P Wave ◽  
S Wave ◽  
P Wave Velocity ◽  
Well Log Analysis

Reflectivities of elastic properties can be expressed as a sum of the reflectivities of P-wave velocity, S-wave velocity, and density, as can the amplitude-variation-with-offset (AVO) parameters, intercept, gradient, and curvature. This common format allows elastic property reflectivities to be expressed as a sum of AVO parameters. Most AVO studies are conducted using a two-term approximation, so it is helpful to reduce the three-term expressions for elastic reflectivities to two by assuming a relationship between P-wave velocity and density. Reduced to two AVO components, elastic property reflectivities can be represented as vectors on intercept-gradient crossplots. Normalizing the lengths of the vectors allows them to serve as basis vectors such that the position of any point in intercept-gradient space can be inferred directly from changes in elastic properties. This provides a direct link between properties commonly used in rock physics and attributes that can be measured from seismic data. The theory is best exploited by constructing new seismic data sets from combinations of intercept and gradient data at various projection angles. Elastic property reflectivity theory can be transferred to the impedance domain to aid in the analysis of well data to help inform the choice of projection angles. Because of the effects of gradient measurement errors, seismic projection angles are unlikely to be the same as theoretical angles or angles derived from well-log analysis, so seismic data will need to be scanned through a range of angles to find the optimum.

scholarly journals Semiglobal viscoacoustic full-waveform inversion

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R271-R293 ◽  
Author(s):  
Nuno V. da Silva ◽  
Gang Yao ◽  
Michael Warner
Keyword(s):  
Physical Properties ◽  
Wave Velocity ◽  
Seismic Data ◽  
Waveform Inversion ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Intrinsic Attenuation ◽  
Full Waveform ◽  
P Wave Velocity

Full-waveform inversion deals with estimating physical properties of the earth’s subsurface by matching simulated to recorded seismic data. Intrinsic attenuation in the medium leads to the dispersion of propagating waves and the absorption of energy — media with this type of rheology are not perfectly elastic. Accounting for that effect is necessary to simulate wave propagation in realistic geologic media, leading to the need to estimate intrinsic attenuation from the seismic data. That increases the complexity of the constitutive laws leading to additional issues related to the ill-posed nature of the inverse problem. In particular, the joint estimation of several physical properties increases the null space of the parameter space, leading to a larger domain of ambiguity and increasing the number of different models that can equally well explain the data. We have evaluated a method for the joint inversion of velocity and intrinsic attenuation using semiglobal inversion; this combines quantum particle-swarm optimization for the estimation of the intrinsic attenuation with nested gradient-descent iterations for the estimation of the P-wave velocity. This approach takes advantage of the fact that some physical properties, and in particular the intrinsic attenuation, can be represented using a reduced basis, substantially decreasing the dimension of the search space. We determine the feasibility of the method and its robustness to ambiguity with 2D synthetic examples. The 3D inversion of a field data set for a geologic medium with transversely isotropic anisotropy in velocity indicates the feasibility of the method for inverting large-scale real seismic data and improving the data fitting. The principal benefits of the semiglobal multiparameter inversion are the recovery of the intrinsic attenuation from the data and the recovery of the true undispersed infinite-frequency P-wave velocity, while mitigating ambiguity between the estimated parameters.

Iterative P-wave velocity inversion using Markov chain Monte Carlo method

2020 ◽  
Author(s):  
Hyunggu Jun ◽  
Hyeong-Tae Jou ◽  
Han-Joon Kim ◽  
Sang Hoon Lee
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
Velocity Structure ◽  
Low Frequency ◽  
P Wave ◽  
Velocity Analysis ◽  
Seismic Section ◽  
Traveltime Tomography ◽  
P Wave Velocity ◽  
Frequency Components

<p>Imaging the subsurface structure through seismic data needs various information and one of the most important information is the subsurface P-wave velocity. The P-wave velocity structure mainly influences on the location of the reflectors during the subsurface imaging, thus many algorithms has been developed to invert the accurate P-wave velocity such as conventional velocity analysis, traveltime tomography, migration velocity analysis (MVA) and full waveform inversion (FWI). Among those methods, conventional velocity analysis and MVA can be widely applied to the seismic data but generate the velocity with low resolution. On the other hands, the traveltime tomography and FWI can invert relatively accurate velocity structure, but they essentially need long offset seismic data containing sufficiently low frequency components. Recently, the stochastic method such as Markov chain Monte Carlo (McMC) inversion was applied to invert the accurate P-wave velocity with the seismic data without long offset or low frequency components. This method uses global optimization instead of local optimization and poststack seismic data instead of prestack seismic data. Therefore, it can avoid the problem of the local minima and limitation of the offset. However, the accuracy of the poststack seismic section directly affects the McMC inversion result. In this study, we tried to overcome the dependency of the McMC inversion on the poststack seismic section and iterative workflow was applied to the McMC inversion to invert the accurate P-wave velocity from the simple background velocity and inaccurate poststack seismic section. The numerical test showed that the suggested method could successfully invert the subsurface P-wave velocity.</p>

scholarly journals Seismic P  Wave Velocity Model From 3‐D Surface and Borehole Seismic Data at the Alpine Fault DFDP‐2 Drill Site (Whataroa, New Zealand)

2021 ◽  
Author(s):  
V Lay ◽  
S Buske ◽  
SB Bodenburg ◽  
John Townend ◽  
R Kellett ◽  
...  
Keyword(s):  
New Zealand ◽  
Wave Velocity ◽  
Seismic Data ◽  
Velocity Model ◽  
P Wave ◽  
Alpine Fault ◽  
P Wave Velocity ◽  
Drill Site ◽  
Borehole Seismic

No description supplied

Sign in / Sign up

Export Citation Format

Share Document