Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPerturbations of the seismic reflectivity of a fluid-saturated depth-dependent poro-elastic medium
489
0
0.0
(
0
)
Authors
Louis De Barros
Ask ChatGPT about the research
No Arabic abstract
Analytical formulas are derived to compute the first-order effects produced by plane inhomogeneities on the point source seismic response of a fluid-filled stratified porous medium. The derivation is achieved by a perturbation analysis of the poro-elastic wave equations in the plane-wave domain using the Born approximation. This approach yields the Frechet derivatives of the P -- SV - and SH-wave responses in terms of the Greens functions of the unperturbed medium. The accuracy and stability of the derived operators are checked by comparing, in the time-distance domain, differential seismograms computed from these analytical expressions with complete solutions obtained by introducing discrete perturbations into the model properties. For vertical and horizontal point forces, it is found that the Frechet derivative approach is remarkably accurate for small and localized perturbations of the medium properties which are consistent with the Born approximation requirements. Furthermore, the first-order formulation appears to be stable at all source-receiver offsets. The porosity, consolidation parameter, solid density and mineral shear modulus emerge as the most sensitive parameters in forward and inverse modeling problems. Finally, the Amplitude-Versus-Angle response of a thin layer shows strong coupling effects between several model parameters.
rate research
Read More
We solve the problem of sparse signal deconvolution in the context of seismic reflectivity inversion, which pertains to high-resolution recovery of the subsurface reflection coefficients. Our formulation employs a nonuniform, non-convex synthesis sparse model comprising a combination of convex and non-convex regularizers, which results in accurate approximations of the l0 pseudo-norm. The resulting iterative algorithm requires the proximal average strategy. When unfolded, the iterations give rise to a learnable proximal average network architecture that can be optimized in a data-driven fashion. We demonstrate the efficacy of the proposed approach through numerical experiments on synthetic 1-D seismic traces and 2-D wedge models in comparison with the benchmark techniques. We also present validations considering the simulated Marmousi2 model as well as real 3-D seismic volume data acquired from the Penobscot 3D survey off the coast of Nova Scotia, Canada.
The coupled mechanics of fluid-filled granular media controls the behavior of many natural systems such as saturated soils, fault gouge, and landslides. The grain motion and the fluid pressure influence each other: It is well established that when the fluid pressure rises, the shear resistance of fluid-filled granular systems decreases, and as a result catastrophic events such as soil liquefaction, earthquakes, and accelerating landslides may be triggered. Alternatively, when the pore pressure drops, the shear resistance of these systems increases. Despite the great importance of the coupled mechanics of grains-fluid systems, the basic physics that controls this coupling is far from understood. We developed a new multi-scaled model based on the discrete element method, coupled with a continuum model of fluid pressure, to explore this dynamical system. The model was shown recently to capture essential feedbacks between porosity changes arising from rearrangement of grains, and local pressure variations due to changing pore configurations. We report here new results from numerical experiments of a continuously shearing layer of circular two-dimensional grains, trapped between two parallel rough boundaries. The experiments use a fixed confining stress on the boundary walls, and a constant velocity applied to one of the boundaries, as if this system was the interior of a sliding geological fault filled with fault gouge. In addition, we control the layer permeability and the drainage boundary conditions. This paper presents modeling results showing that the localization of shear (into a narrow shear band within the shearing layer) is strongly affected by the presence of fluids. While in dry granular layers there is no preferred position for the onset of localization, drained systems tend to localize shear on their boundary. We propose a scaling argument to describe the pressure deviations in a shear band, and use that to predict the allowable positions of shear localizations as a function of the fault and gouge properties.
Crack microgeometries pose a paramount influence on effective elastic characteristics and sonic responses. Geophysical exploration based on seismic methods are widely used to assess and understand the presence of fractures. Numerical simulation as a promising way for this issue, still faces some challenges. With the rapid development of computers and computational techniques, discrete-based numerical approaches with desirable properties have been increasingly developed, but have not yet extensively applied to seismic response simulation for complex fractured media. For this purpose, we apply the coupled LSM-DFN model (Liu and Fu, 2020b) to examining the validity in emulating elastic wave propagation and scattering in naturally-fractured media. By comparing to the theoretical values, the implement of the schema is validated with input parameters optimization. Moreover, dynamic elastic moduli from seismic responses are calculated and compared with static ones from quasi-static loading of uniaxial compression tests. Numerical results are consistent with the tendency of theoretical predictions and available experimental data. It shows the potential for reproducing the seismic responses in complex fractured media and quantitatively investigating the correlations and differences between static and dynamic elastic moduli.
We consider a long-wave transversely isotropic (TI) medium equivalent to a series of finely parallel-layered isotropic layers, obtained using the citet{Backus} average. In such a TI equivalent medium, we verify the citet{Berrymanetal} method of indicating fluids and the authors method citep{Adamus}, using anisotropy parameter $varphi$. Both methods are based on detecting variations of the Lame parameter, $lambda$, in a series of thin isotropic layers, and we treat these variations as potential change of the fluid content. To verify these methods, we use Monte Carlo (MC) simulations; for certain range of Lame parameters $lambda$ and $mu$---relevant to particular type of rocks---we generate numerous combinations of these parameters in thin layers and, after the averaging process, we obtain their TI media counterparts. Subsequently, for each of the aforementioned media, we compute $varphi$ and citet{Thomsen} parameters $epsilon$ and $delta$. We exhibit $varphi$, $epsilon$ and $delta$ in a form of cross-plots and distributions that are relevant to chosen range of $lambda$ and $mu$. We repeat that process for various ranges of Lame parameters. Additionally, to support the MC simulations, we consider several numerical examples of growing $lambda$, by using scale factors. As a result of the thorough analysis of the relations among $varphi$, $epsilon$ and $delta$, we find eleven fluid detectors that compose a new fluid detection method. Based on these detectors, we show the quantified pattern of indicating change of the fluid content.
Using the standard ETAS model of triggered seismicity, we present a rigorous theoretical analysis of the main statistical properties of temporal clusters, defined as the group of events triggered by a given main shock of fixed magnitude m that occurred at the origin of time, at times larger than some present time t. Using the technology of generating probability function (GPF), we derive the explicit expressions for the GPF of the number of future offsprings in a given temporal seismic cluster, defining, in particular, the statistics of the clusters duration and the clusters offsprings maximal magnitudes. We find the remarkable result that the magnitude difference between the largest and second largest event in the future temporal cluster is distributed according to the regular Gutenberg-Richer law that controls the unconditional distribution of earthquake magnitudes. For earthquakes obeying the Omori-Utsu law for the distribution of waiting times between triggering and triggered events, we show that the distribution of the durations of temporal clusters of events of magnitudes above some detection threshold u has a power law tail that is fatter in the non-critical regime $n<1$ than in the critical case n=1. This paradoxical behavior can be rationalised from the fact that generations of all orders cascade very fast in the critical regime and accelerate the temporal decay of the cluster dynamics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
