Reflection moveout inversion in azimuthally anisotropic media: accuracy, limitations and acquisition

Parameter estimation from the elliptical variations in the normal-moveout (NMO) velocity in azimuthally anisotropic media is sensitive to the angular separation between the survey lines in 2D, or equivalently, the source-to-receiver azimuth in 3D, and to the set of azimuths used in the inversion procedure. The accuracy in estimating the orientation of an NMO ellipse, in particular the parameter α, is also sensitive to the magnitude of anisotropy. On the other hand, the accuracy in estimating the semi-axes of the NMO-velocity ellipse is about the same for any magnitude of anisotropy.   To invert for the NMO ellipse parameters at least three NMO-velocity measurements along distinct azimuth directions are needed. In order to maximize the accuracy and stability in parameter estimation, it is best to have the azimuths for the three source-to-receiver directions 60° apart. Having more than three distinct source-to-receiver azimuths (e.g. full azimuthal coverage) provides a useful data redundancy that enhances the quality of the estimates.   In order to maximize quality in the inversion process, it is recommended to design the seismic data acquisition such that it contains small sectors (≤10°) with adequate fold and offset distribution.   Using three NMO-velocity measurements, 60° apart, an azimuthally anisotropic layer overlain by an azimuthally isotropic overburden (as might occur for fractured reservoirs) should have a relative thickness (in time) with respect to the total thickness at least equal to the ratio of the error in the NMO (stacking) velocity to the interval anisotropy of the fractured layer. Coverage along more than three azimuths, however, improves this limitation, which is imposed by Dix differentiation, by at most 50%, depending on the number of observations (NMO velocities) that enter the inversion procedure.     

Processes of initial collision and suturing between India and Asia

The initial collision between Indian and Asian continents marked the starting point for transformation of land-sea thermal contrast, uplift of the Tibet-Himalaya orogen, and climate change in Asia. In this paper, we review the published literatures from the past 30 years in order to draw consensus on the processes of initial collision and suturing that took place between the Indian and Asian plates. Following a comparison of the different methods that have been used to constrain the initial timing of collision, we propose that the tectono-sedimentary response in the peripheral foreland basin provides the most sensitive index of this event, and that paleomagnetism presents independent evidence as an alternative, reliable, and quantitative research method. In contrast to previous studies that have suggested collision between India and Asia started in Pakistan between ca. 55 Ma and 50 Ma and progressively closed eastwards, more recent researches have indicated that this major event first occurred in the center of the Yarlung Tsangpo suture zone (YTSZ) between ca. 65 Ma and 63 Ma and then spreading both eastwards and westwards. While continental collision is a complicated process, including the processes of deformation, sedimentation, metamorphism, and magmatism, different researchers have tended to define the nature of this event based on their own understanding, an intuitive bias that has meant that its initial timing has remained controversial for decades. Here, we recommend the use of reconstructions of each geological event within the orogenic evolution sequence as this will allow interpretation of collision timing on the basis of multidisciplinary methods.     

DHI evaluation by combining rock physics simulation and statistical techniques for fluid identification of Cambrian-to-Cretaceous clastic reservoirs in Pakistan

The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K _sat and K _sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V _P/V _S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K _sat–μ, EEI, V _P/V _S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration.     

Updated earthquake catalogue for seismic hazard analysis in Pakistan

A reliable and homogenized earthquake catalogue is essential for seismic hazard assessment in any area. This article describes the compilation and processing of an updated earthquake catalogue for Pakistan. The earthquake catalogue compiled in this study for the region (quadrangle bounded by the geographical limits 40–83° N and 20–40° E) includes 36,563 earthquake events, which are reported as 4.0–8.3 moment magnitude (M_W) and span from 25 AD to 2016. Relationships are developed between the moment magnitude and body, and surface wave magnitude scales to unify the catalogue in terms of magnitude M_W. The catalogue includes earthquakes from Pakistan and neighbouring countries to minimize the effects of geopolitical boundaries in seismic hazard assessment studies. Earthquakes reported by local and international agencies as well as individual catalogues are included. The proposed catalogue is further used to obtain magnitude of completeness after removal of dependent events by using four different algorithms. Finally, seismicity parameters of the seismic sources are reported, and recommendations are made for seismic hazard assessment studies in Pakistan.     

Seismic fragility functions for code compliant and non-compliant RC SMRF structures in Pakistan

Fragility functions are derived for low-rise code compliant & non-compliant special moment resisting frames (SMRFs). Non-compliant SMRFs those built in low strength concrete and lacking confining ties in joint panel zones, commonly found in developing countries. Shake table tests were performed on single-storey and two-storey 1:3 reduced scale representative frames to understand the damage mechanism and develop deformation-based damage scale. The non-compliant SMRF experienced column flexure cracking, longitudinal bar-slip in beam and observed with cover concrete spalling from the joint panels. The code compliant SMRF experienced flexure cracks in beam/column, and experienced joint cracking under extreme shaking. Numerical modeling technique is developed for inelastic modeling of reinforced concrete frame with beam bar-slip and joint damageability using SeismoStruct. Natural accelerograms were used to analyze the considered frames through incremental dynamic analyses in SeismoStruct. A probabilistic based approach was used to derive fragility functions for the considered frames. An example case study is presented for damageability evaluation of structures for earthquakes of various return periods (43, 72, 475, 2475 years).     

Towards Real‐Time Geodemographics: Clustering Algorithm Performance for Large Multidimensional Spatial Databases

Geodemographic classifications provide discrete indicators of the social, economic and demographic characteristics of people living within small geographic areas. They have hitherto been regarded as products, which are the final “best” outcome that can be achieved using available data and algorithms. However, reduction in computational cost, increased network bandwidths and increasingly accessible spatial data infrastructures have together created the potential for the creation of classifications in near real time within distributed online environments. Yet paramount to the creation of truly real time geodemographic classifications is the ability for software to process and efficiency cluster large multidimensional spatial databases within a timescale that is consistent with online user interaction. To this end, this article evaluates the computational efficiency of a number of clustering algorithms with a view to creating geodemographic classifications “on the fly” at a range of different geographic scales.     

An indicator based approach for assessing household’s perceptions of landslide risk in Murree hills of Pakistan

Natural Hazards - The Murree region in Pakistan is frequently affected by landslides and causes damages to lands, houses, life lines and even loss of livestock and human life. We have attempted in...     

The effectiveness of a pseudo-inverse extended Born operator to handle lateral heterogeneity for imaging and velocity analysis applications

Wave equation–based migration velocity analysis techniques aim to construct a kinematically accurate velocity model for imaging or as an initial model for full waveform inversion applications. The most popular wave equation–based migration velocity analysis method is differential semblance optimization, where the velocity model is iteratively updated by minimizing the unfocused energy in an extended image volume. However, differential semblance optimization suffers from artefacts, courtesy of the adjoint operator used in imaging, leading to poor convergence. Recent findings show that true amplitude imaging plays a significant role in enhancing the differential semblance optimization's gradient and reducing the artefacts. Here, we focus on a pseudo-inverse operator to the horizontally extended Born as a true amplitude imaging operator. For laterally inhomogeneous models, the operator required a derivative with respect to a vertical shift. Extending the image vertically to evaluate such a derivative is costly and impractical. The inverse operator can be simplified in laterally homogeneous models. We derive an extension of the approach to apply the full inverse formula and evaluate the derivative efficiently. We simplified the implementation by applying the derivative to the imaging condition and utilize the relationship between the source and receiver wavefields and the vertical shift. Specifically, we verify the effectiveness of the approach using the Marmousi model and show that the term required for the lateral inhomogeneity treatment has a relatively small impact on the results for many cases. We then apply the operator in differential semblance optimization and invert for an accurate macro-velocity model, which can serve as an initial velocity model for full waveform inversion.     

Mapping the change of coral reefs using remote sensing and in situ measurements: a case study in Pangkajene and Kepulauan Regency,Spermonde Archipelago,Indonesia

As elsewhere in Indonesia, local inhabitants in the Pangkajene and Kepulauan (PANGKEP) Regency, Spermonde Archipelago area and along the south-west coast of Sulawesi traditionally regard the coral reefs as their livelihood source. Since human activities as well as natural disturbances pose major threats to the coral reefs, these livelihoods may also be at risk. Currently, no comprehensive information on the status and condition of coral reefs in this area is available for this resource management. We determined the changes of coral reef habitat over a period of 20 years from 1994 to 2014 using a satellite Landsat multi-temporal image substantiated with in situ measurement data collected in 2014. The spectral value of coral reefs was extracted from multi-temporal Landsat imagery data, while the diffuse attenuation coefficient of water was obtained by using statistical analysis between the ratio of live coral cover and the spectral value of the visible bands. By using the unsupervised classification integrated with the data ground truth, it is stated that there has been a decline in live coral cover over a period of 20 years from 7716 ha in 1994 to 4236 ha in 2014, with a degradation rate of 174 ha/year. Based on the results, the ratio of the coral cover in the coral reef transects varied from the average of 24% for live corals to 96% for coral rubbles, implying the degraded status of coral reefs in the study area.