Surveys in Geophysics最新文献

Sparse representation and inversion have been widely used in the acquisition and processing of geophysical data. In particular, the low-rank representation of seismic signals shows that they can be determined by a few elementary modes with predominantly large singular values. We review global and local low-rank representation for seismic reflectivity models and then apply it to least-squares migration (LSM) in acoustic and viscoacoustic media. In the global singular value decomposition (SVD), the elementary modes determined by singular vectors represent horizontal and vertical stratigraphic segments sorted from low to high wavenumbers, and the corresponding singular values reflect the contribution of these basic modes to form a broadband reflectivity model. In contrast, local SVD for grouped patch matrices can capture nonlocal similarity and thus accurately represent the reflectivity model with fewer ranks than the global SVD method. Taking advantage of this favorable sparsity, we introduce a local low-rank regularization into LSM to estimate subsurface reflectivity models. A two-step algorithm is developed to solve this low-rank constrained inverse problem: the first step is for least-squares data fitting and the second is for weighted nuclear-norm minimization. Numerical experiments for synthetic and field data demonstrate that the low-rank constraint outperforms conventional shaping and total-variation regularizations, and can produce high-quality reflectivity images for complicated structures and low signal-to-noise data.

The rheological models of Lomnitz and Jeffreys have been widely used in earthquake seismology (to simulate a nearly constant Q medium) and to describe the creep and relaxation behavior of rocks as a function of time. Other similar models, such as those of Becker, Scott Blair and Kolsky, show similar properties, particularly the Scott Blair model describes a perfectly constant Q as a function of frequency. We first give a historical overview of the main scientists and the development and versions of the various models and priorities of discovery. Then, we clarify the relationship between the different versions of these models in terms of mathematical expressions of the complex modulus and calculate the phase velocity and quality factor Q as a function of frequency, illustrating the various special cases. In addition, we give useful hints for the numerical calculation of these moduli, which include special cases of the hypergeometric function.

The accurate extraction of useful signals from the measurement data is one of the important parts and challenges of the understanding of subsurface information. The desired signal is usually hidden in the background noise, and the amplitude is weak due to the particular geological environment of the subsurface or the consistency of the measuring instrument. By extending the theory of seismic interferometry to include the effects of adjacent virtual channels and by combining super- and reverse-virtual interferometry, we obtain a hybrid virtual interferometry technique based on stacking neighboring virtual traces for wave reconstruction. We have verified the effectiveness of the processing method in suppressing noise interference and extracting useful signals using synthetic data tests. The method is applied to the processing and interpretation of acoustic measurements acquired in a cased borehole of a coal formation and an open hole of an igneous formation, where the processed waveforms are finely reconstructed, and the estimated slowness results are in good agreement with other measurements, thus providing an effective tool for data analysis.

The Rodl-Kaplice-Blanice fault system (RKB) of Variscan shear origin, repeatedly active since the Late Paleozoic to the Recent, is expressed by a number of lithological contacts, distinct geophysical gradients and many landforms. A general trend of the RKB as well as linear configuration of its internal architecture is fairly similar to those of topical near Rhine Graben and Alpine-Carpathian transition area as the two other consistent recently reactivated large-scale tectonic structures in the extended (thinned) crust of central Europe. In middle part of the RKB, the occurring linear topographic and geological features parallel to the main RKB sections point to the existence of a wide tectonic zone in the crust following the fault system. Our multidisciplinary study includes a summary of corresponding basic geological data, overview of seismic, regional geophysical and geomorphological conditions, primary model of recent kinematic activity in the RKB area derived from the space (Global Navigation Satellite System—GNSS) monitoring and terrestrial (repeated high precision levelling) geodetic data and comparison of these various information.

The obtained knowledge indicates that the RKB is active up to ~ 1.0 mm horizontally and > 0.5 mm vertically. The fault system area in the Bohemian Massif can be subdivided into the three parts of diverse tectonic structure and block kinematics. Sinistral horizontal movements are highest near the southern surface sections (Rodl-Kaplice, Rudolfov and Drahotěšice faults), whereas noticeable vertical differentiation is going on mainly along the Blanice and Kouřim faults in the north where the RKB activity is gradually decreasing towards the extensive Elbe shear zone with transverse movements. The middle part of the RKB is dislocated by a large active transverse tectonic structure of the South Bohemian Basins (SBB) with variable horizontal velocity vectors of surface GNSS stations. Most of the weak regional earthquakes have been recorded west of the RKB. Besides faults of the SBB, these were mainly associated with the RKB-subparallel Lhenice fault. Based on the earthquake distribution and foci depths, the latter fault can have similar structural position as the RKB related to lower part of the Variscan level in the ~ 10–12 km depth.

For polyhedral mass bodies having arbitrary shape and density distribution of polynomial type we present a tensorial approach to derive analytical expressions of the gravitational potential and gravity vector. They are evaluated at an arbitrary point by means of formulas, referred to a Cartesian reference frame having an arbitrary origin, that are shown to be singularity-free whatever is the position of the observation point with respect to the body. The solution is expressed as a sum of algebraic quantities depending solely upon the 3D coordinates of the polyhedron vertices and the coefficients of the polynomial density function. Hence, no recursive expression needs to be invoked as in the recent contribution by Ren et al. (Surv Geophys 41:695–722, 2020). Moreover, the tensorial formalism developed in the paper allows one to obtain more concise, coordinate-free expressions that can also be extended to address polynomial functions of greater order. The analytical expressions of the gravitational potential and gravity vector are numerically validated and compared with alternative methods retrieved from the literature.

The occurrence of subsurface karst caves can cause the development of superficial depressions which, in turn, may pose a construction hazard. Identifying such a substratum requires integrated non-invasive measurement methods. The main objective of the study was to demonstrate the effectiveness of the non-invasive ERT, TLS, and GPR survey techniques in identifying the karst floor and determining the direction of discontinuities around the cave. The paper analyzes the limitations of the methods used in the study of heterogeneous media. These limitations are related to the methodology and measurement conditions, data processing, and interpretation in the context of the resolution and depth range. The study was conducted using the example of the Jaskinia pod Świecami cave, formed in the Sarmatianal calcarenites in Poland. The research confirmed its complex karst-anthropogenic genesis. The cave was formed as a result of the infiltration of rainwater and the dissolution of limestone by groundwater, while the paleokarst forms that are characteristic of it and of the surrounding caves and occur in their vicinity, i.e., narrow ridges called "karst candles", were formed as a result of water circulation during the local permafrost degradation in the middle Pleistocene. However, these forms were modified in the Upper Pleistocene and Holocene, as indicated by ERT images.

Compared with surface wave corresponding to the normal mode, which is widely studied, there is less research on guided-P wave corresponding to the leaking mode. Guided-P wave carries the dispersion information that can be used to construct the subsurface velocity structures. In this paper, to simultaneously estimate P-wave velocity (({{v}}_{{P}})) and S-wave velocity (({{v}}_{{S}})) structures, an integrated inversion method of guided-P and surface wave dispersion curves is proposed. Through the calculation of Jacobian matrix, the sensitivity of dispersion curves is quantitatively analyzed. It shows that the dispersion curves of guided-P and surface waves are, respectively, sensitive to the ({{v}}_{{P}}) and ({{v}}_{{S}}). Synthetic model tests demonstrate the proposed integrated inversion method can estimate the ({{v}}_{{P}}) and ({{v}}_{{S}}) models accurately and effectively identify low-velocity interlayers. The integrated inversion method is also applied to the field seismic data acquired for oil and gas prospecting. The pseudo-2D ({{v}}_{{P}}), ({{v}}_{{S}}) and Poisson’s ratio inversion results are of significance for near-surface geological interpretation. The comparison with the result of first-arrival traveltime tomography further demonstrates the accuracy and practicality of the proposed integrated inversion method. Not only in the field of exploration seismic, the guided-P wave dispersion information can also be extracted from the earthquake seismic, engineering seismic and ambient noise. The proposed inversion method can exploit previously neglected guided-P wave to characterize the subsurface ({{v}}_{{P}}) structures, showing broad and promising application prospects. This compensates for the inherent defect that the surface wave dispersion curve is mainly sensitive to the ({{v}}_{{S}}) structure.