SH-SH wave inversion for S-wave velocity and density

SS-waves (SV-SV waves and SH-SH waves) are capable of inverting S-wave velocity ( VS) and density ( ρ) because they are sensitive to both parameters. SH-SH waves can be separated from multicomponent data sets more effectively than the SV-SV wave because the former is decoupled from the PP-wave in isotropic media. In addition, the SH-SH wave can be better modeled than the SV-SV wave in the case of strong velocity/impedance contrast because the SV-SV wave has multicritical angles, some of which can be quite small when velocity/ impedance contrast is strong. We derived an approximate equation of the SH-SH wave reflection coefficient as a function of VS and ρ in natural logarithm variables. The approximation has high accuracy, and it enables the inversion of VS and ρ in a direct manner. Both coefficients corresponding to VS and ρ are “model-parameter independent” and thus there is no need for prior estimate of any model parameter in inversion. Then, we developed an SH-SH wave inversion method, and demonstrated it by using synthetic data sets and a real SH-SH wave prestack data set from the west of China. We found that VS and ρ can be reliably estimated from the SH-SH wave of small angles.

SEISMIC RESONANCE: WAVELET-FREE REFLECTIVITY RETRIEVAL VIA MODIFIED CEPSTRAL DECOMPOSITION

Spectral decomposition is a proven tool in seismic interpretation, aiding interpreters to highlight channels, map temporal bed thickness and other geological discontinuities. Once seismic data is spectrally decomposed, notch patterns in the amplitude spectra are indicative of the reservoir layer’s thickness and/or its interval velocity. Additional cepstral decomposition will allow direct extraction of bed time-thickness or arrival time under particular reflectivity series setup. We build on these observations to establish a more generalized workflow for reflectivity retrieval method without the need to understand the details of the wavelet, provided the starting seismic is stably phased via phase correction during processing. We demonstrate reflector time and its ‘apparent strength’ can be identified in a transformed seismic resonance domain resulted from a modified cepstrum analysis. In this domain, each reflector can be characterized from obvious linear hot spots. The timing and strength of those linear hot spots will reveal reflector times and scaled reflectivity coefficients. This new method is subsequently applied for thickness prediction of a target reservoir in a complex geological setting, with large thickness variations and weak impedance contrast with underlying lithology previously complicating identification of base-reservoir. In a deep-water field blind test, the sand thicknesses evaluated from this method are found to be close to true vertical thickness found in wells.

Seismic microzonation of a region with complex surficial geology based on different site classification approaches

A seismic microzonation study was conducted to refine the seismic hazard model for the city of Saguenay, Canada. The Quaternary geology underlying Saguenay shows complex glacial and post-glacial stratigraphy with a number of buried valleys filled with fluvioglacial and glaciomarine sediments. High impedance contrast between rock formations and surficial sediments is prone to seismic amplification. To evaluate their applicability, advantages and limitations in capturing the geological specificity of the study area, four site classification methods were applied: the current National Building Code of Canada (NBCC) and Eurocode 8, both mainly based on the average shear-wave velocity for the surficial sediments (VS,avg) and for the top 30 m (VS,30); a method based on the fundamental site period (T0); and a hybrid method based on the combination of VS,30, T0 and VS,avg. The study specifically aimed to evaluate the importance of the site classification parameters on the resulting microzonation maps. VS,30 is capable to present the geological and geotechnical site conditions, however, the results may be further improved by considering Vs,avg in shallow and T0 in thick layers of soil sediments as secondary parameters. The T0 method gives also satisfactory results with T0 showing a better correlation to Vs,30 than to Vs,avg. The versatile hybrid method may be challenging to apply in certain cases with its nine different site categories and parameters.

Characterization of an Unstable Slope Using Geophysical, UAV, and Geological Techniques: Karakoram Himalaya, Northern Pakistan

Given active tectonism, rough terrain, and climate, the mountainous ranges in northern Pakistan are prone to geohazards, including earthquakes, unstable slopes, and landslides. The frequent landsliding in the region poses a risk to communities, economic activities, and transportation networks. In this context, the unstable slope above Mayun village calls for a multi-method approach for better assessment of the slope for planning interventions aimed at hazard mitigation. We conducted an integrated study including uncrewed aerial vehicle (UAV) and ground-penetrating radar (GPR) in coordination with geomorphic field observations to image the possible slip surfaces for a comprehensive understanding of a potential future rockslide with significant socioeconomic consequences. UAV-derived results helped delineate the overall extent of the unstable slope and its downslope area in a quick, remote, and safe way. GPR profiles have enabled the reconstruction of the bedrock’s morphology and its internal structure and the depth distribution of cracks running through the overburden and bedrock. The results provided insight into the stable and unstable compartments of the slope due to the thin cover of surficial deposits, high impedance contrast at the overburden-bedrock interface, lateral heterogeneities, and presence of open cracks, and almost detached blocks, respectively. These data on the dynamic properties of a landslide-prone slope could be used for the correct planning of civil infrastructure to minimize the potential risk of building damage in the seismically active Hunza valley.

Geophysical constraints on the properties of a subglacial lake in northwest Greenland

In this study, we report the results of an active-source seismology and ground-penetrating radar survey performed in northwestern Greenland at a site where the presence of a subglacial lake beneath the accumulation area has previously been proposed. Both seismic and radar results show a flat reflector approximately 830–845 m below the surface, with a seismic reflection coefficient of −0.43 ± 0.17, which is consistent with the acoustic impedance contrast between a layer of water and glacial ice. Additionally, in the seismic data we observe an intermittent lake bottom reflection arriving between 14–20 ms after the lake top reflection, corresponding to a lake depth of approximately 10–15 m. A strong coda following the lake top and lake bottom reflections is consistent with a package of lake bottom sediments although its thickness and material properties are uncertain. Finally, we use these results to conduct a first-order assessment of the lake origins using a one-dimensional thermal model and hydropotential modeling based on published surface and bed topography. Using these analyses, we narrow the lake origin hypotheses to either anomalously high geothermal flux or hypersalinity due to local ancient evaporite. Because the origins are still unclear, this site provides an intriguing opportunity for the first in situ sampling of a subglacial lake in Greenland, which could better constrain mechanisms of subglacial lake formation, evolution, and relative importance to glacial hydrology.

Modeling the effects of fracture infill on frequency-dependent anisotropy and AVO response of a fractured porous layer

In a fractured porous hydrocarbon reservoir, wave velocities and reflections depend on frequency and incident angle. A proper description of the frequency dependence of amplitude variations with offset (AVO) signatures should allow effects of fracture infills and attenuation and dispersion of fractured media. The novelty of this study lies in the introduction of an improved approach for the investigation of incident-angle and frequency variations-associated reflection responses. The improved AVO modeling method, using a frequency-domain propagator matrix method, is feasible to accurately consider velocity dispersion predicted from frequency-dependent elasticities from a rock physics modeling. And hence, the method is suitable for use in the case of an anisotropic medium with aligned fractures. Additionally, the proposed modeling approach allows the combined contributions of layer thickness, interbedded structure, impedance contrast and interferences to frequency-dependent reflection coefficients and, hence, yielding seismograms of a layered model with a dispersive and attenuative reservoir. Our numerical results show bulk modulus of fracture fluid significantly affects anisotropic attenuation, hence causing frequency-dependent reflection abnormalities. These implications indicate the study of amplitude versus angle and frequency (AVAF) variations provides insights for better interpretation of reflection anomalies and hydrocarbon identification in a layered reservoir with vertical transverse isotropy (VTI) dispersive media.

Compact Wideband Microstrip Patch Antenna Design for Breast Cancer Detection

The current breast cancer detection techniques are mostly invasive and suffer from high cost, high false rate and inefficacy in early detection. These limitations can be subdued by development of non-invasive microwave detection system whose performance is predominantly dependent on the antenna used in the system. The designing of a compact wideband antenna and matching its impedance with breast phantom is a challenging task. In this paper, we have designed a compact antenna matched with the breast phantom operating in wideband frequency from 1 to 6 GHz capable to detect the dielectric (or impedance) contrast of the benign and malignant tissue. The impedance of the antenna is matched to a cubically shaped breast phantom and a very small tumor (volume=1 cm3). The antenna is tuned to the possible range of electrical properties of breast phantom and tumour (permittivity ranging from 10 to 20 and conductivity from 1.5 to 2.5 S/m). The return loss (S11), E-field distribution and specific absorption rate (SAR) are simulated. The operating band of antenna placed near the phantom without tumor was found to be (1.11-5.47)GHz and with tumor inside phantom is (1.29-5.50)GHz. Results also show that the SAR of the antenna is within the safety limit.

Evaluation of site amplification factors for shallow rock sites of Islamabad, Pakistan.

Current seismic provisions in Building Code of Pakistan (BCP-2007), seismic site classifications, and corresponding site amplification factors were determined similar to Uniform Building Code (UBC-1997), which were based on the local site conditions of Western United States with low impedance contrast, deep sites, and high seismicity conditions. The use of these deep sites based amplification factors to the shallow sites may not be appropriate to capture shallow sites' response. In this study, we performed a non-linear site-specific response analysis suite to compute the amplification factors for the six (06) representative shallow bedrock sites of Islamabad, Pakistan. The computed amplification factors are compared with BCP-2007 code-based design estimates; it is found that BCP-2007 code-based guidelines underestimate amplification factors at a short period whereas overestimating at an extended period. The study's findings highlight the potential implications of current code-based amplification factors for Islamabad's shallow sites and suggest improving the current seismic guidelines.

Broadband Q-Factor Imaging for Geofluid Detection in the Gulf of Trieste (Northern Adriatic Sea)

Seismic surveys allow estimating lithological parameters, as P-wave velocity and anelastic absorption, which can detect the presence of fracture and fluids in the geological formations. Recently, a new method has been proposed for high-resolution imaging of anelastic absorption, which combines a macro-model from seismic tomography with a micro-model obtained by the pre-stack depth migration of a seismic attribute, i.e., the instantaneous frequency. As a result, we can get a broadband image that provides clues about the presence of saturating fluids. When the saturation changes sharply, as for gas reservoirs with an impermeable caprock, the acoustic impedance contrast produces “bright spots” because of the resulting high reflectivity at its top. When the fluid content changes smoothly, the anelastic absorption becomes a good detector, as fluid-filled formations absorb more seismic energy than hard rocks. We apply this method for imaging the anelastic absorption in a regional seismic survey acquired by OGS in the Gulf of Trieste (northern Adriatic Sea, Italy).

Mississippian Chat and Tripolite Zones in Osage County, Oklahoma: Paleokarst Interpretation based on 3D Seismic and Well Logs

Mississippian paleokarst Chat and tripolitic chert (tripolite) zones associated with the Mississippian Lime have been hydrocarbon exploration targets in Osage County for many decades. Chat is residual chert, either in place or transported, weathered out of chert-bearing Mississippian Limestone that was eroded at the Mississippian-Pennsylvanian unconformity. Thus, the formation of Chat is an epigene paleokarst process. Tripolite occurs as a highly porous, silica-rich interval within the Mississippi Lime. It is formed by in-place alteration of limestone by silica-rich surface waters or deep-seated hydrothermal fluids, making tripolite formation a mixed or hypogene paleokarst process. Here, we have studied Chat and tripolite by seismic analysis calibrated by well control with full-wave sonic and density log data. We identify that Chat and tripolite show clear separation from Mississippian Lime log-based acoustic impedance and VP/ VS, but there is no meaningful separation of Chat from tripolite, and they both exhibit total porosities greater than 20% with evidence of fracture porosity. We find that the sonic-based normal incidence wedge model for Chat bounded above by the Pennsylvanian Shale and below by the Mississippian Lime indicate that two seismic expressions are plausible: first, a strong negative amplitude when the Chat thickness is above the tuning thickness (for this survey it is 56 ft) and, second, a weak positive or negative amplitude associated with the small impedance contrast between Chat and overlying Pennsylvanian Shale. Our analysis suggests that the traditional Chat “strong response” and a new “dim-out” exploration strategy may be usefully applied in Osage County. We show that the tripolite response is consistently a negative amplitude event that strengthens with increasing tripolite thickness. We provide an interpretive framework for characterizing Chat and tripolite zones associated with the Mississippian Lime in the US Midcontinent, which may apply to regions around the world.