# Journal of Geophysics and Engineering - Latest Papers

Latest articles for Journal of Geophysics and Engineering
Updated: 5 hours 29 min ago

### Using marine magnetic survey data to identify a gold ore-controlling fault: a case study in Sanshandao fault, eastern China

Sun, 02/25/2018 - 7:00pm
The Jiaodong Peninsula has the greatest concentration of gold ore in China and is characterized by altered tectonite-type gold ore deposits. This type of gold deposit is mainly formed in fracture zones and is strictly controlled by faults. Three major ore-controlling faults occur in the Jiaodong Peninsula—the Jiaojia, Zhaoping and Sanshandao faults; the former two are located on land and the latter is located near Sanshandao and its adjacent offshore area. The discovery of the world’s largest marine gold deposit in northeastern Sanshandao indicates that the shallow offshore area has great potential for gold prospecting. However, as two ends of the Sanshandao fault extend to the Bohai Sea, conventional geological survey methods cannot determine the distribution of the fault and this is constraining the discovery of new gold deposits. To explore the southwestward extension of the Sanshandao fault, we performed a 1:25 000 scale marine magnetic survey in this region and obtained hig...

### Spatio-temporal evolution of apparent resistivity during coal-seam hydraulic flushing

Sun, 02/25/2018 - 7:00pm
Hydraulic flushing in gas predrainage is widely used, but the hydraulic-flushing effect is evaluated in a traditional way, by determining the desorption volume, moisture content, gas drainage rate and other conventional indices. To verify the rationality and feasibility of the multielectrode resistivity method in the evaluation of coal-seam hydraulic flushing and to research the spatio-temporal evolution of apparent resistivity during hydraulic flushing, a field test was conducted in 17# coal seam at Nuodong Mine, Guizhou. During hydraulic flushing, four stages were defined according to the variation in coal rock resistivity with time, namely, the preparation stage, the sharply decreasing stage, the rapidly increasing stage and the steady stage. The apparent resistivity of the coal rock mass is affected mainly by its own degree of fragmentation and flushing volume. A more serious rupture and a greater flushing volume yield a smaller apparent resistivity during the sharply decrea...

### Numerical solution of fractured horizontal wells in shale gas reservoirs considering multiple transport mechanisms

Sun, 02/25/2018 - 7:00pm
The multiscale pore size and specific gas storage mechanism in organic-rich shale gas reservoirs make gas transport in such reservoirs complicated. Therefore, a model that fully incorporates all transport mechanisms and employs an accurate numerical method is urgently needed to simulate the gas production process. In this paper, a unified model of apparent permeability was first developed, which took into account multiple influential factors including slip flow, Knudsen diffusion (KD), surface diffusion, effects of the adsorbed layer, permeability stress sensitivity, and ad-/desorption phenomena. Subsequently, a comprehensive mathematical model, which included the model of apparent permeability, was derived to describe gas production behaviors. Thereafter, on the basis of unstructured perpendicular bisection grids and finite volume method, a fully implicit numerical simulator was developed using Matlab software. The validation and application of the new model were confirmed usin...

### Simultaneous multi-component seismic denoising and reconstruction via K-SVD

Thu, 02/22/2018 - 7:00pm
Data denoising and reconstruction play an increasingly significant role in seismic prospecting for their value in enhancing effective signals, dealing with surface obstacles and reducing acquisition costs. In this paper, we propose a novel method to denoise and reconstruct multicomponent seismic data simultaneously. This method lies within the framework of machine learning and the key points are defining a suitable weight function and a modified inner product operator. The purpose of these two processes are to perform missing data machine learning when the random noise deviation is unknown, and building a mathematical relationship for each component to incorporate all the information of multi-component data. Two examples, using synthetic and real multicomponent data, demonstrate that the new method is a feasible alternative for multi-component seismic data processing.

### Evaluation of subsoil competence for foundation studies at site III of the Delta State University, Nigeria

Tue, 02/20/2018 - 7:00pm
An integrated geophysical and geotechnical investigation has been carried out at site III of the Delta State University, Abraka, Nigeria. This took place in a bid to generate information on the competence of the soil in withstanding stress and strain emanating from overburden or pore pressure, swelling, cracking and other anthropogenic activity in relation to civil engineering and building structures. An electromagnetic method employing the very low frequency (VLF) technique, and electrical resistivity employing the Wenner and the vertical electrical sounding techniques were used for this study. Soil samples were also collected at depth for geotechnical analysis. Isoresistivity slices generated from the data of 33 VES stations at 1 m showed generally low resistivity values of subsurface earth materials, classified as clayey sand, sandy clay or clay, and ranging from 60–300 Ωm. However, at depths of 3 and 5 m, the result showed a generally high resistivity distribution with value...

### SNR enhancement for downhole microseismic data based on scale classification shearlet transform

Tue, 02/20/2018 - 7:00pm
Shearlet transform (ST) can be effective in 2D signal processing, due to its parabolic scaling, high directional sensitivity, and optimal sparsity. ST combined with thresholding has been successfully applied to suppress random noise. However, because of the low magnitude and high frequency of a downhole microseismic signal, the coefficient values of valid signals and noise are similar in the shearlet domain. As a result, it is difficult to use for denoising. In this paper, we present a scale classification ST to solve this problem. The ST is used to decompose noisy microseismic data into serval scales. By analyzing the spectrum and energy distribution of the shearlet coefficients of microseismic data, we divide the scales into two types: low-frequency scales which contain less useful signal and high-frequency scales which contain more useful signal. After classification, we use two different methods to deal with the coefficients on different scales. For the low-frequency scales,...

### A theoretical approach to quantify the effect of random cracks on rock deformation in uniaxial compression

Tue, 02/20/2018 - 7:00pm
Cracks have a significant effect on the uniaxial compression of rocks. Thus, a theoretically analytical approach was proposed to assess the effects of randomly distributed cracks on the effective Young’s modulus during the uniaxial compression of rocks. Each stage of the rock failure during uniaxial compression was analyzed and classified. The analytical approach for the effective Young’s modulus of a rock with only a single crack was derived while considering the three crack states under stress, namely, opening, closure-sliding, and closure-nonsliding. The rock was then assumed to have many cracks with randomly distributed direction, and the effect of crack shape and number during each stage of the uniaxial compression on the effective Young’s modulus was considered. Thus, the approach for the effective Young’s modulus was used to obtain the whole stress–strain process of uniaxial compression. Afterward, the proposed approach was employed to analyze the effects of related param...

### Improving fault image by determination of optimum seismic survey parameters using ray-based modeling

Tue, 02/20/2018 - 7:00pm
In complex structures such as faults, salt domes and reefs, specifying the survey parameters is more challenging and critical owing to the complicated wave field behavior involved in such structures. In the petroleum industry, detecting faults has become crucial for reservoir potential where faults can act as traps for hydrocarbon. In this regard, seismic survey modeling is employed to construct a model close to the real structure, and obtain very realistic synthetic seismic data. Seismic modeling software, the velocity model and parameters pre-determined by conventional methods enable a seismic survey designer to run a shot-by-shot virtual survey operation. A reliable velocity model of structures can be constructed by integrating the 2D seismic data, geological reports and the well information. The effects of various survey designs can be investigated by the analysis of illumination maps and flower plots. Also, seismic processing of the synthetic data output can describe the ta...

### Hybrid sparse blind deconvolution: an implementation of SOOT algorithm to real data

Tue, 02/20/2018 - 7:00pm
Getting information of seismic data depends on deconvolution as an important processing step; it provides the reflectivity series by signal compression. This compression can be obtained by removing the wavelet effects on the traces. The recently blind deconvolution has provided reliable performance for sparse signal recovery. In this study, two deconvolution methods have been implemented to the seismic data; the convolution of these methods provides a robust spiking deconvolution approach. This hybrid deconvolution is applied using the sparse deconvolution (MM algorithm) and the Smoothed-One-Over-Two algorithm (SOOT) in a chain. The MM algorithm is based on the minimization of the cost function defined by standards l1 and l2. After applying the two algorithms to the seismic data, the SOOT algorithm provided well-compressed data with a higher resolution than the MM algorithm. The SOOT algorithm requires initial values to be applied for real data, such as the wavelet coefficients ...

### Angle-domain common imaging gather extraction via Kirchhoff prestack depth migration based on a traveltime table in transversely isotropic media

Thu, 02/15/2018 - 7:00pm
Angle-domain common image-point gathers (ADCIGs) can alleviate the limitations of common image-point gathers in an offset domain, and have been widely used for velocity inversion and amplitude variation with angle (AVA) analysis. We propose an effective algorithm for generating ADCIGs in transversely isotropic (TI) media based on the gradient of traveltime by Kirchhoff pre-stack depth migration (KPSDM), as the dynamic programming method for computing the traveltime in TI media would not suffer from the limitation of shadow zones and traveltime interpolation. Meanwhile, we present a specific implementation strategy for ADCIG extraction via KPSDM. Three major steps are included in the presented strategy: (1) traveltime computation using a dynamic programming approach in TI media; (2) slowness vector calculation by the gradient of a traveltime table calculated previously; (3) construction of illumination vectors and subsurface angles in the migration process. Numerical examples are...

### Experimental and numerical study on frost heave of saturated rock under uniform freezing conditions

Thu, 02/15/2018 - 7:00pm
A series of freezing experiments are conducted on saturated sandstone and mortar specimens to investigate the frost heave of saturated rock under uniform freezing conditions. The experimental results show that the frost heave of saturated rock is isotropic under uniform freezing conditions. During the freezing process, three stages are observed in the curves of variation of total frost heaving strain versus time: the thermal contraction stage, the frost heaving stage and the steady stage. Moreover, the amount of final stable frost heave first increases and then decreases with decrease in freezing temperature, and the maximum final stable frost heave occurs at different freezing temperature in saturated sandstone and mortar. Furthermore, a coupled thermal–mechanical (TM) model of frost heave of saturated rock is proposed in which a constraint coefficient ##IMG## [http://ej.iop.org/images/1742-2140/15/2/593/jgeaa93acieqn1.gif] {$\zeta$} is used to consider t...

### The effect of grain size and cement content on index properties of weakly solidified artificial sandstones

Thu, 02/15/2018 - 7:00pm
The effects of textural characteristics, especially grain size, on index properties of weakly solidified artificial sandstones are studied. For this purpose, a relatively large number of laboratory tests were carried out on artificial sandstones that were produced in the laboratory. The prepared samples represent fifteen sandstone types consisting of five different median grain sizes and three different cement contents. Indices rock properties including effective porosity, bulk density, point load strength index, and Schmidt hammer values (SHVs) were determined. Experimental results showed that the grain size has significant effects on index properties of weakly solidified sandstones. The porosity of samples is inversely related to the grain size and decreases linearly as grain size increases. While a direct relationship was observed between grain size and dry bulk density, as bulk density increased with increasing median grain size. Furthermore, it was observed that the point l...

### Model identification and control of development of deeply buried paleokarst reservoir in the central Tarim Basin, northwest China

Thu, 02/15/2018 - 7:00pm
The paleokarst reservoirs of the Ordovician Yingshan formation, rich in oil and gas, are deeply buried in the central Tarim Basin, northwest China. Dozens of imaging well-logs in this region reveal five typical paleokarst features, including solution vugs, solution-enlarged fractures, filled caves, unfilled caves and collapsed caves, as well as two typical paleokarst structures located in different paleotopographic sites, including paleokarst vadose and phreatic zones. For seismic data, the large wave impedance contrast between the paleocave system and the surrounding rocks leads to a strong seismic reflection, which is highlighted as a bead-like ‘bright spot’ in a seismic section. By quantitatively estimating the seismic resolution limits of deep seismic reflections, a single paleocave cannot be identified from a seismic profile, and the bead-like reflection represents an entire paleocave complex. The spectral decomposition technique was employed to depict the planar shape and ...

### Seismic rock physics characterization of anisotropic shale—a Longmaxi Shale case study

Mon, 02/12/2018 - 7:00pm
An anisotropic rock physics model based on effective medium theories is constructed for shales and applied to the Longmaxi Shale formation. In the rock physics model, intrinsic VTI (transverse isotropy with a vertical asymmetry axis) anisotropy due to preferred orientations of clay particles is quantified by introducing the clay lamination ( CL ) index in the Backus averaging method, and additional anisotropy enhanced by bed-parallel fractures is described by the Chapman model for multi-scale pore spaces. Rock physics templates are constructed based on the model in order to better understand the implicit relations between reservoir properties, elastic and mechanical properties and seismic parameters. A model-based method is proposed for inversion of the CL and aspect ratio ( AR ) of horizontal fractures, as well as VTI anisotropy parameters for the shale formation. Results indicate that the obtained CL index shows a negative correlation with clay content,...

### Accuracy and sensitivity analysis on seismic anisotropy parameter estimation

Mon, 02/12/2018 - 7:00pm
There is significant uncertainty in measuring the Thomsen’s parameter δ in laboratory even though the dimensions and orientations of the rock samples are known. It is expected that more challenges will be encountered in the estimating of the seismic anisotropy parameters from field seismic data. Based on Monte Carlo simulation of vertical transversely isotropic layer cake model using the database of laboratory anisotropy measurement from the literature, we apply the commonly used quartic non-hyperbolic reflection moveout equation to estimate the seismic anisotropy parameters and test its accuracy and sensitivities to the source-receive offset, vertical interval velocity error and time picking error. The testing results show that the methodology works perfectly for noise-free synthetic data with short spread length. However, this method is extremely sensitive to the time picking error caused by mild random noises, and it requires the spread length to be greater than the de...

### Experimental analysis and application of the effect of stress on continental shale reservoir brittleness

Mon, 02/12/2018 - 7:00pm
Hydraulic fracturing is an effective measure of reservoir modification for the development of shale gas. The evaluation of rock brittleness can provide a basis for the optimization of fracturing. In this paper, the effect of stress on the brittleness of shale is systematically analyzed by designing triaxial mechanics tests. The strain analysis method was used to evaluate the shale brittleness. The research indicates that, with the increase of effective confining pressure, the value of the brittleness index ( B 1 ) decreases. There is a linear and positive correlation between the average reduction ratio of B 1 and the buried depth. The stress has a significant effect on the shale brittleness. Therefore, the rock brittleness can be overestimated without considering the influence of the buried depth or the stress of formation when using the mineral composition method. Being affected by the stress, when the brittle mineral content of the shale reser...

### An effective absorbing layer for the boundary condition in acoustic seismic wave simulation

Mon, 02/12/2018 - 7:00pm
Efficient numerical simulation of seismic wavefields generally involves truncating the Earth model in order to keep computing time and memory requirements down. Absorbing boundary conditions, therefore, are applied to remove the boundary reflections caused by this truncation, thereby allowing for accurate modeling of wavefields. In this paper, we derive an effective absorbing boundary condition for both acoustic and elastic wave simulation, through the simplification of the damping term of the split perfectly matched layer (SPML) boundary condition. This new boundary condition is accurate, cost-effective, and easily implemented, especially for high-performance computing. Stability analysis shows that this boundary condition is effectively as stable as normal (non-absorbing) wave equations for explicit time-stepping finite differences. We found that for full-waveform inversion (FWI), the strengths of the effective absorbing layer—a reduction of the computational and memory cost c...

### Comparative study of receiver-side ghost wavefield attenuation on different marine acquisition configurations

Mon, 02/12/2018 - 7:00pm
In marine seismic exploration, the ghost energies (down-going waves), which arise from the reflection at the surface, are often treated as unwanted signals for data processing. The ghost wave fields interfere with the desired primary signals, leads to frequency notches and attenuation of low frequencies, which in turn downgrade the resolution of the recorded seismic data. There are two main categories of methods to solve the ghost or the so-called notch problem: the non-conventional acquisition configuration-based technique and a deghosting algorithm-based solution. The variable-depth streamer (VDS) acquisition solution is one of the most representative methods in the first category, which has become a popular solution for marine seismic acquisition to obtain broad data bandwidth. However, this approach is not as economic as the conventional constant depth streamer (CDS) acquisition, due to the precise control of the towing streamer. In addition, there are large quantities of co...

### Sensitivity analyses of acoustic impedance inversion with full-waveform inversion

Mon, 02/12/2018 - 7:00pm
Acoustic impedance estimation has a significant importance to seismic exploration. In this paper, we use full-waveform inversion to recover the impedance from seismic data, and analyze the sensitivity of the acoustic impedance with respect to the source-receiver offset of seismic data and to the initial velocity model. We parameterize the acoustic wave equation with velocity and impedance, and demonstrate three key aspects of acoustic impedance inversion. First, short-offset data are most suitable for acoustic impedance inversion. Second, acoustic impedance inversion is more compatible with the data generated by density contrasts than velocity contrasts. Finally, acoustic impedance inversion requires the starting velocity model to be very accurate for achieving a high-quality inversion. Based upon these observations, we propose a workflow for acoustic impedance inversion as: (1) building a background velocity model with travel-time tomography or reflection waveform inversion; (2...

### Improved vertical displacements induced by a refined thermal expansion model and its quantitative analysis in GPS height time series

Mon, 02/12/2018 - 7:00pm
There are apparent seasonal variations in GPS height time series, and thermal expansion is considered to be one of the potential geophysical contributors. The displacements introduced by thermal expansion are usually derived without considering the annex height and underground part of the monument (e.g. located on roof or top of the buildings), which may bias the geophysical explanation of the seasonal oscillation. In this paper, the improved vertical displacements are derived by a refined thermal expansion model where the annex height and underground depth of the monument are taken into account, and then 560 IGS stations are adopted to validate the modeled thermal expansion (MTE) displacements. In order to evaluate the impact of thermal expansion on GPS heights, the MTE displacements of 80 IGS stations with less data discontinuities are selected to compare with their observed GPS vertical (OGV) displacements with the modeled surface loading (MSL) displacements removed in advanc...