Grid dispersion in generating finite-differences synthetic seismograms

Ｇｒｉｄｄｉｓｐｅｒｓｉｏｎｉｎｇｅｎｅｒａｔｉｎｇｆｉｎｉｔｅ－ｄｉｆｆｅｒｅｎｃｅｓｓｙｎｔｈｅｔｉｃｓｅｉｓｍｏｇｒａｍｓＡｂｄｏｌｒａｈｉｍＪａｖａｈｅｒｉａｎ（ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ，ｔｈｅＵｎｉｖｅｒｓｉｔｙｏｆＴｅｈ...     

Nonstationary deconvolution using maximum kurtosis optimization

Deconvolution is an essential step for high-resolution imaging in seismic data processing. The frequency and phase of the seismic wavelet change through time during wave propagation as a consequence of seismic absorption. Therefore, wavelet estimation is the most vital step of deconvolution, which plays the main role in seismic processing and inversion. Gabor deconvolution is an effective method to eliminate attenuation effects. Since Gabor transform does not prepare the information about the phase, minimum-phase assumption is usually supposed to estimate the phase of the wavelet. This manner does not return the optimum response where the source wavelet would be dominantly a mixed phase. We used the kurtosis maximization algorithm to estimate the phase of the wavelet. First, we removed the attenuation effect in the Gabor domain and computed the amplitude spectrum of the source wavelet; then, we rotated the seismic trace with a constant phase to reach the maximum kurtosis. This procedure was repeated in moving windows to obtain the time-varying phase changes. After that, the propagating wavelet was generated to solve the inversion problem of the convolutional model. We showed that the assumption of minimum phase does not reflect a suitable response in the case of mixed-phase wavelets. Application of this algorithm on synthetic and real data shows that subtle reflectivity information could be recovered and vertical seismic resolution is significantly improved.     

The optimum source－time function for generating finite－difference synthetic seismograms

Ｔｈｅｏｐｔｉｍｕｍｓｏｕｒｃｅ－ｔｉｍｅｆｕｎｃｔｉｏｎｆｏｒｇｅｎｅｒａｔｉｎｇｆｉｎｉｔｅ－ｄｉｆｆｅｒｅｎｃｅｓｙｎｔｈｅｔｉｃｓｅｉｓｍｏｇｒａｍｓＡｂｄｏｌｒａｈｉｍＪＡＶＡＨＥＲＩＡＮ（ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ，ｔｈｅ...     

Development of flood mitigation strategies toward sustainable development

Natural Hazards - Decreasing flood damages in a basin and achieving sustainable development reveal the necessity of integrated management of flood. Integrated flood management requires the...     

Ground-roll attenuation using modified common-offset–common-reflection-surface stacking

We modified the common-offset–common-reflection-surface (COCRS) method to attenuate ground roll, the coherent noise typically generated by a low-velocity, low-frequency, and high-amplitude Rayleigh wave. The COCRS operator is based on hyperbolas, thus it fits events with hyperbolic traveltimes such as reflection events in prestack data. Conversely, ground roll is linear in the common-midpoint (CMP) and common-shot gathers and can be distinguished and attenuated by the COCRS operator. Thus, we search for the dip and curvature of the reflections in the common-shot gathers prior to the common-offset section. Because it is desirable to minimize the damage to the reflection amplitudes, we only stack the multicoverage data in the ground-roll areas. Searching the CS gathers before the CO section is another modification of the conventional COCRS stacking. We tested the proposed method using synthetic and real data sets from western Iran. The results of the ground-roll attenuation with the proposed method were compared with results of the f–k filtering and conventional COCRS stacking after f–k filtering. The results show that the proposed method attenuates the aliased and nonaliased ground roll better than the f–k filtering and conventional CRS stacking. However, the computation time was higher than other common methods such as f–k filtering.     

Seismic channel edge detection using 3D shearlets—a study on synthetic and real channelised 3D seismic data

Automatic feature detection from seismic data is a demanding task in today's interpretation workstations. Channels are among important stratigraphic features in seismic data both due to their reservoir capability or drilling hazard potential. Shearlet transform as a multi‐scale and multi‐directional transformation is capable of detecting anisotropic singularities in two and higher dimensional data. Channels occur as edges in seismic data, which can be detected based on maximizing the shearlet coefficients through all sub‐volumes at the finest scale of decomposition. The detected edges may require further refinement through the application of a thinning methodology. In this study, a three‐dimensional, pyramid‐adapted, compactly supported shearlet transform was applied to synthetic and real channelised, three‐dimensional post‐stack seismic data in order to decompose the data into different scales and directions for the purpose of channel boundary detection. In order to be able to compare the edge detection results based on three‐dimensional shearlet transform with some famous gradient‐based edge detectors, such as Sobel and Canny, a thresholding scheme is necessary. In both synthetic and real data examples, the three‐dimensional shearlet edge detection algorithm outperformed Sobel and Canny operators even in the presence of Gaussian random noise.     

Automatic Method for Identifying Photospheric Bright Points and Granules Observed by Sunrise

In this study, we propose methods for the automatic detection of photospheric features (bright points and granules) from ultra-violet (UV) radiation, using a feature-based classifier. The methods use quiet-Sun observations at 214 nm and 525 nm images taken by Sunrise on 9 June 2009. The function of region growing and mean shift procedure are applied to segment the bright points (BPs) and granules, respectively. Zernike moments of each region are computed. The Zernike moments of BPs, granules, and other features are distinctive enough to be separated using a support vector machine (SVM) classifier. The size distribution of BPs can be fitted with a power-law slope ?1.5. The peak value of granule sizes is found to be about 0.5 arcsec². The mean value of the filling factor of BPs is 0.01, and for granules it is 0.51. There is a critical scale for granules so that small granules with sizes smaller than 2.5 arcsec² cover a wide range of brightness, while the brightness of large granules approaches unity. The mean value of BP brightness fluctuations is estimated to be 1.2, while for granules it is 0.22. Mean values of the horizontal velocities of an individual BP and an individual BP within the network were found to be 1.6 km?s^?1 and 0.9 km?s^?1, respectively. We conclude that the effect of individual BPs in releasing energy to the photosphere and maybe the upper layers is stronger than what the individual BPs release into the network.     

Seismic time–frequency analysis as a robust method to estimate the fluid saturation: A case study of carbonates reservoir,Iran

The propagation of seismic waves through a saturated reservoir compresses the fluid in the pore spaces. During this transition, parts of seismic energy would be attenuated because of intrinsic absorption. Rock physics models make the bridge between the seismic properties and petrophysical reality in the earth. Attenuation is one of the significant seismic attributes used to describe the fluid behaviour in the reservoirs. We examined the core samples using ultrasonic experiments at the reservoir conditions. Given the rock properties of the carbonate reservoir and experiment results, the patchy saturation mechanism was solved for substituted fluid using the theory of modulus frequency. The extracted relationship between the seismic attenuation and water saturation was used in time–frequency analysis. We performed the peak frequency method to estimate the Q factor in the Gabor domain and determined the water saturation based on the computed rock physics model. The results showed how the probable fault in the reservoir has stopped the fluid movement in the reservoir and caused touching the water‐bearing zone through drilling.     

Study of blasting vibrations in Sarcheshmeh copper mine

Introduction In spite of development of mechanized methods of ground excavation, drilling and blasting is still extensively employed because of its low capital investment and simplicity. Its extensive use is not even limited by extension of mines close to residential areas and vital establishments. If it is not used in a controlled way, blasting operation can cause instability, failure of mine slopes and damps and damage to the nearby structures. The main objective here could be to reduce the …