A discontinuous Galerkin method for seismic wave propagation in coupled elastic and poroelastic media     

Yijie Zhang  Jinghuai Gao  Weimin Han  Yanbin He 《Geophysical Prospecting》2019,67(5):1392-1403

Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement.  相似文献   

The kinematics of S waves in acoustic orthorhombic media     

Shibo Xu  Alexey Stovas 《Geophysical Prospecting》2019,67(8):1983-2002

Orthorhombic models are often used in the seismic industry nowadays to describe azimuthal and polar anisotropy and reasonably realistic in capturing the features of the earth interior. It is challenging to handle so many model parameters in the seismic data processing. In order to reduce the number of the parameters for P wave, the acoustic orthorhombic medium is proposed by setting all on-axis S wave velocities to zero. However, due to the coupled behaviour for P and S waves in the orthorhombic model, the ‘S wave artefacts’ are still remained in the acoustic orthorhombic model, which kinematics needs to be defined and analysed. In this paper, we analyse the behaviour of S wave in acoustic orthorhombic media. By analysis of the slowness surface in acoustic orthorhombic media, we define the S waves (or S wave artefacts) that are more complicated in shape comparing to the one propagating in an acoustic transversely isotropic medium with a vertical symmetry axis. The kinematic properties of these waves are defined and analysed in both phase and group domain. The caustics, amplitude and the multi-layered case for S wave in acoustic orthorhombic model are also discussed. It is shown that there are two waves propagating in this acoustic orthorhombic medium. One of these waves is similar to the one propagating in acoustic vertical symmetry axis media, whereas another one has a very complicated shape consisting of two crossing surfaces.  相似文献   

A fast and robust two-point ray tracing method in layered media with constant or linearly varying layer velocity     

Xinding Fang  Xiaofei Chen 《Geophysical Prospecting》2019,67(7):1811-1824

A fast and robust method for two-point ray tracing in one-dimensional layered media is presented. This method is applicable to layered models with constant or linearly varying isotropic layer velocity. For given model properties and source and receiver positions, a ray path can be uniquely determined once its ray parameter (i.e. horizontal slowness) is known. The ray parameter can be obtained by numerically solving the nonlinear offset (i.e. source–receiver horizontal distance) equation using Newton's method, which generally works well at near and mid offsets. However, Newton's method becomes hard to converge at large offsets due to the oversensitivity of offset to ray parameter. Based on the analysis of the characteristic of the offset equation, a modified ray parameter is proposed and used to replace the generic ray parameter in numerical calculation. Numerical experiments show that the iteration process becomes stable and converges rapidly with the modified ray parameter. Moreover, a rational function that asymptotically approximates the shape of the offset equation is introduced for obtaining good initial estimates of the modified ray parameter. Numerical tests show that this method is robust in any situation, and an accurate ray parameter can be obtained within two or three iterations for a wide range of model velocity structure and source–receiver distance. Furthermore, the proposed two-point ray tracing method is easy to implement.  相似文献   

Rock elasticity as a function of the uniaxial stress: laboratory measurements and theoretical modelling of vertical transversely isotropic and orthorhombic shales     

V.A. Sviridov  S.I. Mayr  S.A. Shapiro 《Geophysical Prospecting》2019,67(7):1867-1881

We apply a rock-physics model that describes the relationship between the effective stress and rock elasticity. We experimentally obtain and analyse a data set containing one vertical transversely isotropic and one orthorhombic shale sample. The vertical transversely isotropic symmetry of the first sample is caused by the layered structure of the rock. The seismic orthorhombicity of the second sample could be explained after microscopic analysis of thin section, which demonstrates an imperfect disorder of inhomogeneities. Both samples were loaded uniaxially in a quasi-static regime. During the loading, we measured stress-dependent seismic velocities and sample deformations. For the analysis of the stress-dependent velocities and stiffnesses, we modelled the measured data set using a recent generalization of the porosity deformation approach. Comparison of the experimentally determined and numerically modelled data supports the applicability of the theory and helps in the interpretation of experimentally obtained data. In agreement with the theory, uniaxial stress increases the elliptic component of the seismic anisotropy and does not impact the anellipticity parameter. We demonstrate the distinct influence of the stiff and compliant porosities on the stress sensitivity of the elastic properties.  相似文献   

Forward plane-wave electromagnetic model in three dimensions using hybrid finite volume–integral equation scheme     

Musa A. Bello  Rongwen Guo  Jianxin Liu 《Geophysical Prospecting》2019,67(8):2213-2226

We present a concept of the hybrid finite volume–integral equation technique for solving Maxwell's equation in a quasi-static form. The divergence correction was incorporated to improve the convergence and stability of the governing linear system equations which pose a challenge on the discretization of the curl–curl Helmholtz equation. A staggered finite volume approach is applied for discretizing the system of equations on a structured mesh and solved in a secondary field technique. The bi-conjugate gradient stabilizer was utilized with block incomplete lower-upper factorization preconditioner to solve the system of equation. To obtain the electric and magnetic fields at the receivers, we use the integral Green tensor scheme. We verify the strength of our hybrid technique with benchmark models relative to other numerical algorithms. Importantly, from the tested models, our scheme was in close agreement with the semi-analytical solution. It also revealed that the use of a quasi-analytical boundary condition helps to minimize the runtime for the linear system equation. Furthermore, the integral Green tensor approach to compute at the receivers demonstrates better accuracy compared with the conventional interpolation method. This adopted technique can be applied efficiently to the inversion procedure.  相似文献   

Three-dimensional geophysical characterization of the La Rambla and Zafra de Záncara anticlines (Loranca Basin,Central Spain)     

C. Ayala  F.M. Rubio  C. Rey-Moral  M.I. Reguera  C. Biete 《Geophysical Prospecting》2019,67(3):580-594

The Zafra de Záncara anticline (also known as the El Hito anticline), located in the Loranca Cenozoic Basin (part of the Tagus Basin, Central Spain), had been studied by several oil companies during the late 1960s and early 1970s. In 2009, within the ‘Plan for selection and characterization of suitable structures of CO₂ geological storage’, this anticline was selected as a potential CO₂ storage site. A preliminary three-dimensional geological model, based on five geological cross sections that were constrained with the interpretation of the available seismic profiles (that are rather old and do not have very good quality), was created. With the aim of improving the geological knowledge of the Zafra de Záncara anticline and helping to investigate the suitability of a nearby anticline, namely La Rambla, as another structural closure that might make it a possible CO₂ storage site, a local gravity survey (1 station every km²) was carried out in the area, seven new geological cross sections, based on these existing seismic profiles and field geology, were build, and a new three-dimensional geological model that included both anticlines, improved through three-dimensional stochastic gravity inversion, was constructed. The densities needed for the geological formations of the model come from the analysis of rock samples, logging data from El Hito-1 drillhole and petrophysical information from Instituto Geológico y Minero de España database. The inversion has improved the knowledge about the geometry of the anticlines’ traps and seals as well as the geometry of the basement relief and the structural relationship between basement and cover.  相似文献   

Water retention effects on elastic properties of Opalinus shale     

Alexey Yurikov  Maxim Lebedev  Marina Pervukhina  Boris Gurevich 《Geophysical Prospecting》2019,67(4):984-996

Shales play an important role in many engineering applications such as nuclear waste, CO₂ storage and oil or gas production. Shales are often utilized as an impermeable seal or an unconventional reservoir. For both situations, shales are often studied using seismic waves. Elastic properties of shales strongly depend on their hydration, which can lead to substantial structural changes. Thus, in order to explore shaly formations with seismic methods, it is necessary to understand the dependency of shale elastic properties on variations in hydration. In this work, we investigate structural changes in Opalinus shale at different hydration states using laboratory measurements and X-ray micro-computed tomography. We show that the shale swells with hydration and shrinks with drying with no visible damage. The pore space of the shale deforms, exhibiting a reduction in the total porosity with drying and an increase in the total porosity with hydration. We study the elastic properties of the shale at different hydration states using ultrasonic velocities measurements. The elastic moduli of the shale show substantial changes with variations in hydration, which cannot be explained with a single driving mechanism. We suggest that changes of the elastic moduli with variations in hydration are driven by multiple competing factors: (1) variations in total porosity, (2) substitution of pore-filling fluid, (3) change in stiffness of contacts between clay particles and (4) chemical hardening/softening of clay particles. We qualitatively and quantitatively analyse and discuss the influence of each of these factors on the elastic moduli. We conclude that depending on the microstructure and composition of a particular shale, some of the factors dominate over the others, resulting in different dependencies of the elastic moduli on hydration.  相似文献   

Feasibility of utilizing wavelet phase to map the CO2 plume at the Ketzin pilot site,Germany          下载免费PDF全文

Fei Huang  Christopher Juhlin  Li Han  Daniel Sopher  Monika Ivandic  Ben Norden  Wubing Deng  Fengjiao Zhang  Thomas Kempka  Stefan Lüth 《Geophysical Prospecting》2017,65(2):523-543

Spectral decomposition is a powerful tool that can provide geological details dependent upon discrete frequencies. Complex spectral decomposition using inversion strategies differs from conventional spectral decomposition methods in that it produces not only frequency information but also wavelet phase information. This method was applied to a time‐lapse three‐dimensional seismic dataset in order to test the feasibility of using wavelet phase changes to detect and map injected carbon dioxide within the reservoir at the Ketzin carbon dioxide storage site, Germany. Simplified zero‐offset forward modelling was used to help verify the effectiveness of this technique and to better understand the wavelet phase response from the highly heterogeneous storage reservoir and carbon dioxide plume. Ambient noise and signal‐to‐noise ratios were calculated from the raw data to determine the extracted wavelet phase. Strong noise caused by rainfall and the assumed spatial distribution of sandstone channels in the reservoir could be correlated with phase anomalies. Qualitative and quantitative results indicate that the wavelet phase extracted by the complex spectral decomposition technique has great potential as a practical and feasible tool for carbon dioxide detection at the Ketzin pilot site.  相似文献   

Basement geology beneath the northeast Thelon Basin,Nunavut: insights from integrating new gravity,magnetic and geological data   总被引：1，自引：0，他引：1       下载免费PDF全文

V. Tschirhart  C.W. Jefferson  W.A. Morris 《Geophysical Prospecting》2017,65(2):617-636

Current models for unconformity‐associated uranium deposits predict fluid flow and ore deposition along reactivated faults in >1.76 Ga basement beneath Mesoproterozoic siliciclastic basins. In frontier regions such as the Thelon Basin in the Kivalliq region of Nunavut, little is known about the sub‐basin distribution of units and structures, making exploration targeting very tenuous. We constructed a geological map of the basement beneath the unconformity by extrapolating exposed features into the subsurface. The new map is constrained by detailed geological, geophysical, and rock property observations of outcrops adjacent to the basin and by aeromagnetic and gravity data over the geophysically transparent sedimentary basin. From rock property measurements, it is clear that the diverse magnetic and density characteristics of major rock packages provide quantitative three‐dimensional constraints. Gravity profiles forward modelled in four cross sections define broad synforms of the Amer Belt and Archean volcanic rocks that are consistent with the structural style outside the basin. Major lithotectonic entities beneath the unconformity include: supracrustal rocks of the Archean Woodburn Lake group and Marjorie Hills meta sedimentary gneiss and associated mixed granitoid and amphibolitic gneiss; the Amer Mylonite Zone and inferred mafic intrusions oriented parallel and sub‐parallel; other igneous intrusions of 2.6 Ga, 1.83 Ga, and 1.75 Ga vintage; and the <2.3 Ga to >1.84 Ga Amer Group. Four main brittle regional fault arrays (040°–060°, 075°–90°, 120°, and 150°) controlled development and preservation of the basin. The reactivated intersections of such faults along fertile basement units such as the Rumble assemblage, Marjorie Hills assemblage, Nueltin igneous rocks, and Pitz formation are the best targets for uranium exploration.  相似文献   

Modelling Chorabari Lake outburst flood,Kedarnath, India     

Mohammd RAFIQ  Shakil Ahmad ROMSHOO  Anoop Kumar MISHRA  Faizan JALAL 《山地科学学报》2019,(1):64-76

In this study, the Glacier Lake Outburst Flood(GLOF) that occurred over Kedarnath in June 2013 was modeled using integrated observations from the field and Remote Sensing(RS). The lake breach parameters such as area, depth, breach, and height have been estimated from the field observations and Remote Sensing(RS) data. A number of modelling approaches, including Snow Melt Runoff Model(SRM), Modified Single Flow model(MSF), Watershed Management System(WMS), Simplified Dam Breach Model(SMPDBK) and BREACH were used to model the GLOF. Estimations from SRM produced a runoff of about 22.7 m3 during 16–17, June 2013 over Chorabari Lake. Bathymetry data reported that the lake got filled to its maximum capacity(3822.7 m3) due to excess discharge. Hydrograph obtained from the BREACH model revealed a peak discharge of about 1699 m3/s during an intense water flow episode that lasted for 10–15 minutes on 17 th June 2013. Excess discharge from heavy rainfall and snowmelt into the lake increased its hydrostatic pressure and the lake breached cataclysmically.  相似文献