No Arabic abstract
The analysis of surface wave dispersion curves is a way to infer the vertical distribution of shear-wave velocity. The range of applicability is extremely wide going, for example, from seismological studies to geotechnical characterizations and exploration geophysics. However, the inversion of the dispersion curves is severely ill-posed and only limited efforts have been put into the development of effective regularization strategies. In particular, relatively simple smoothing regularization terms are commonly used, even when this is in contrast with the expected features of the investigated targets. To tackle this problem, stochastic approaches can be utilized, but they are too computationally expensive to be practical, at least, in the case of large surveys. Instead, within a deterministic framework, we evaluate the applicability of a regularizer capable of providing reconstructions characterized by tunable levels of sparsity. This adjustable stabilizer is based on the minimum support regularization, applied before on other kinds of geophysical measurements, but never on surface wave data. We demonstrate the effectiveness of this stabilizer on i) two benchmark - publicly available - datasets at crustal and near-surface scales, ii) an experimental dataset collected on a well-characterized site. In addition, we discuss a possible strategy for the estimation of the depth of investigation. This strategy relies on the integrated sensitivity kernel used for the inversion and calculated for each individual propagation mode. Moreover, we discuss the reliability, and possible caveats, of the direct interpretation of this particular estimation of the depth of investigation, especially in the presence of sharp boundary reconstructions.
Geologic shear fractures such as faults and slip surfaces involve marked friction along the discontinuities as they are subjected to significant confining pressures. This friction plays a critical role in the growth of these shear fractures, as revealed by the fracture mechanics theory of Palmer and Rice decades ago. In this paper, we develop a novel phase-field model of shear fracture in pressure-sensitive geomaterials, honoring the role of friction in the fracture propagation mechanism. Building on a recently proposed phase-field method for frictional interfaces, we formulate a set of governing equations for different contact conditions (or lack thereof) in which frictional energy dissipation emerges in the crack driving force during slip. We then derive the degradation function and the threshold fracture energy of the phase-field model such that the stress-strain behavior is insensitive to the length parameter for phase-field regularization. This derivation procedure extends a methodology used in recent phase-field models of cohesive tensile fracture to shear fracture in frictional materials in which peak and residual strengths coexist and evolve by confining pressure. The resulting phase-field formulation is demonstrably consistent with the theory of Palmer and Rice. Numerical examples showcase that the proposed phase-field model is a physically sound and numerically efficient method for simulating shear fracture processes in geologic materials, such as faulting and slip surface growth.
Time-lapse seismic monitoring of carbon storage and sequestration is often challenging because the time-lapse signature of the growth of CO2 plumes is weak in amplitude and therefore difficult to detect seismically. This situation is compounded by the fact that the surveys are often coarsely sampled and not replicated to reduce costs. As a result, images obtained for different vintages (baseline and monitor surveys) often contain artifacts that may be attributed wrongly to time-lapse changes. To address these issues, we propose to invert the baseline and monitor surveys jointly. By using the joint recovery model, we exploit information shared between multiple time-lapse surveys. Contrary to other time-lapse methods, our approach does not rely on replicating the surveys to detect time-lapse changes. To illustrate this advantage, we present a numerical sensitivity study where CO2 is injected in a realistic synthetic model. This model is representative of the geology in the southeast of the North Sea, an area currently considered for carbon sequestration. Our example demonstrates that the joint recovery model improves the quality of time-lapse images allowing us to monitor the CO2 plume seismically.
We report experimental observations of two canonical surface wave patterns --- ship waves and ring waves --- skewed by sub-surface shear, thus confirming effects predicted by recent theory. Observed ring waves on a still surface with sub-surface shear current are strikingly asymmetric, an effect of strongly anisotropic wave dispersion. Ship waves for motion across a sub--surface current on a still surface exhibit striking asymmetry about the ships line of motion, and large differences in wake angle and transverse wavelength for upstream vs downstream motion are demonstrated, all of which in good agreement with theoretical predictions. Neither of these phenomena can occur on a depth-uniform current. A quantitative comparison of measured vs predicted average phase shift for a ring wave is grossly mispredicted by no-shear theory, but in good agreement with predictions for the measured shear current. A clear difference in wave frequency within the ring wave packet is observed in the upstream vs downstream direction for all shear flows, while it conforms with theory for quiescent water for propagation normal to the shear current, as expected. Peak values of the measured 2-dimensional Fourier spectrum for ship waves are shown to agree well with the predicted criterion of stationary ship waves, with the exception of some cases where results are imperfect due to the limited wave-number resolution, transient effects and/or experimental noise. Experiments were performed on controlled shear currents created in two different ways, with a curved mesh, and beneath a blocked stagnant-surface flow. Velocity profiles were measured with particle image velocimetry, and surface waves with a synthetic schlieren method. Our observations lend strong empirical support to recent predictions that wave forces on vessels and structures can be greatly affected by shear in estuarine and tidal waters.
Seismic data quality is vital to geophysical applications, so methods of data recovery, including denoising and interpolation, are common initial steps in the seismic data processing flow. We present a method to perform simultaneous interpolation and denoising, which is based on double-sparsity dictionary learning. This extends previous work that was for denoising only. The original double sparsity dictionary learning algorithm is modified to track the traces with missing data by defining a masking operator that is integrated into the sparse representation of the dictionary. A weighted low-rank approximation algorithm is adopted to handle the dictionary updating as a sparse recovery optimization problem constrained by the masking operator. Compared to traditional sparse transforms with fixed dictionaries that lack the ability to adapt to complex data structures, the double-sparsity dictionary learning method learns the signal adaptively from selected patches of the corrupted seismic data while preserving compact forward and inverse transform operators. Numerical experiments on synthetic seismic data indicate that this new method preserves more subtle features in the dataset without introducing pseudo-Gibbs artifacts when compared to other directional multiscale transform methods such as curvelets.
We investigate the transport problem that a spinful matter wave is incident on a strong localized spin-orbit-coupled Bose-Einstein condensate in optical lattices, where the localization is admitted by atom interaction only existing at one particular site, and the spin-orbit coupling arouse spatial rotation of the spin texture. We find that tuning the spin orientation of the localized Bose-Einstein condensate can lead to spin-nonreciprocal / spin-reciprocal transport, meaning the transport properties are dependent on / independent of the spin orientation of incident waves. In the former case, we obtain the conditions to achieve transparency, beam-splitting, and blockade of the incident wave with a given spin orientation, and furthermore the ones to perfectly isolate incident waves of different spin orientation, while in the latter, we obtain the condition to maximize the conversion of different spin states. The result may be useful to develop a novel spinful matter wave valve that integrates spin switcher, beam-splitter, isolator, and converter. The method can also be applied to other real systems, e.g., realizing perfect isolation of spin states in magnetism, which is otherwise rather difficult.